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launcher,  the  freed  weight  can  then  be  reallocated  to  the 
Cube  Satellites.  This  will  increase  the  maximum  allowed 
weight  per  satellite. 

The  current  design  has  six  pounds  of  available  free 
weight.  During  the  course  of  the  analysis,  it  is  found  that 
removal  of  unused  portions  of  the  four  side  walls  of  the 
launcher  structure  will  reduce  the  weight  without 
undermining  the  strength.  Additionally,  when  looking  at  the 
use  of  different  metal  options  for  the  launcher,  it  is  found 
that  the  structure  needs  to  maintain  the  use  of  an  aluminum 
alloy.  Maintaining  space-grade  aluminum  in  conjunction  with 
the  removal  of  unused  sections  of  the  side  walls  frees  an 
additional  four  pounds.  This  totals  ten  pounds  that  can  be 
allocated  back  into  the  Cube  Satellites. 


v 


THIS  PAGE  INTENTIONALLY  LEFT  BLANK 


vi 


TABLE  OF  CONTENTS 


I.  HISTORY  AND  PURPOSE . 1 

A.  SATELLITES . 1 

B.  STUDENT  INVOLVEMENT . 4 

C.  NPSCUL . 6 

II.  THE  DYNAMIC  STRUCTURE . 11 

A.  LAUNCH  EFFECTS . 11 

B.  MODAL  SHAPES . 11 

C.  VIBRATIONS  TESTING . 14 

D.  OBJECTIVE . 17 

III.  SIMPLIFICATION  OF  THE  STRUCTURE . 19 

A.  SINGLE  WALL . 19 

B.  COMPARISON  TECHNIQUE . 20 

C.  FLAT  PLATE  ANALYSIS . 20 

IV.  MASS  REDUCTION  OF  THE  ALUMINUM  WALL . 23 

A.  REMOVAL  OF  WALL  VOLUME . 23 

B.  TRUSS  ANALYSIS . 25 

C.  CALCULATIONS . 27 

D.  OPTIMAL  CONFIGURATION . 35 

V.  CHANGE  OF  MATERIAL . 39 

A.  CURRENT  MATERIAL . 39 

B.  OTHER  MATERIALS . 39 

C.  MINITAB . 40 

VI.  COMPOUND  LOAD  ANALYSIS . 45 

A.  ADDITIONAL  ANALYSIS . 45 

B.  MATLAB  RESULTS . 47 

VII.  CONCLUSION . 49 

A.  CONFIGURATION . 49 

B.  ALUMINUM  7075  AND  COMBINED  LOADING . 49 

C.  FUTURE  WORK . 50 

APPENDIX  A . 51 

APPENDIX  B . 123 

APPENDIX  C . 127 

APPENDIX  D . 129 

APPENDIX  E . 131 

APPENDIX  F . 165 

APPENDIX  G . 167 


Vll 


LIST  OF  REFERENCES . 169 

INITIAL  DISTRIBUTION  LIST . 173 


viii 


LIST  OF  FIGURES 


Figure  1.  JCSAT-12  being  built.  From  [2] . 2 

Figure  2.  Size  of  a  Cube  Satellite.  From  [5] . 4 

Figure  3.  Cal  Poly  Mk  II  P-POD.  From  [7] . 5 

Figure  4.  Image  of  NPSCuL . 7 

Figure  5.  ESPA  ring  on  the  EELV.  From  [13] . 8 

Figure  6.  ULA  ABC.  From  [14] . 8 

Figure  7.  Mode  1  analyses.  From  [17] . 12 

Figure  8.  Mode  2  analyses.  From  [17] . 12 

Figure  9.  Mode  3  analyses.  From  [17] . 13 

Figure  10.  Mode  6  analyses.  From  [17] . 13 

Figure  11.  First  natural  frequency,  42.96  Hz,  structural 

deflection.  From  [10] . 15 

Figure  12.  "Rocking"  natural  frequency,  57.83  Hz, 

structural  deflection.  From  [10] . 16 

Figure  13.  "Rocking"  natural  frequency,  59.32  Hz, 

structural  deflection.  From  [10] . 16 

Figure  14.  "Po-go"  natural  frequency,  100.50  Hz, 

structural  deflection.  From  [10] . 17 

Figure  15.  Two-dimensional  view  of  side  wall,  including 

the  modifiable  area  (MA) . 19 

Figure  16.  Two-dimensional  view  of  side  wall  with 

sequencer  bolt  holes . 22 

Figure  17.  Visual  display  of  full  wall  versus  removing  the 

MA . 24 

Figure  18.  Truss  variations  A,  B  and  C . 26 

Figure  19.  Truss  variations  D,  E,  and  F . 26 

Figure  20.  Truss  variations  G  and  H . 27 

Figure  21.  Example  of  what  different  beam  widths  look  like 

in  the  MA . 2  8 

Figure  22.  Reduction  of  stiffness  versus  mass  for  0.25- 

inch  beam  width . 31 

Figure  23.  Reduction  of  stiffness  versus  mass  for  0.5-inch 

beam  width . 31 

Figure  24.  Reduction  of  stiffness  versus  mass  for  1-inch 

beam  width . 32 

Figure  25.  Total  freed  mass  from  all  truss  types  at  0.25 

inches . 33 

Figure  26.  Total  freed  mass  from  all  truss  types  at  0.5 

inches . 33 

Figure  27.  Total  freed  mass  from  all  truss  types  at  1.0 

inches . 34 

Figure  28.  Comparison  of  four  favored  solutions  for  0.25- 

inch  beam  width . 36 


IX 


Figure  29.  Comparison  of  four  favored  solutions  for  0.5- 

inch  beam  width . 36 

Figure  30.  Comparison  of  four  favored  solutions  for  1-inch 

beam  width . 37 

Figure  31.  Parameters  used  for  MINITAB  input . 41 

Figure  32.  MINITAB  results  of  holding  the  beam  width  at 

0.25  inches . 42 

Figure  33.  MINITAB  results  of  holding  the  beam  width  at 

0.5  inches . 43 

Figure  34 .  MINITAB  results  of  holding  the  beam  width  at  1 

inch . 43 

Figure  35.  Diagrams  of  compound  forces  analyzed  (Fl-Flb) . . .45 

Figure  36.  Diagrams  of  compound  forces  analyzed  (Flc-Fle) .. 46 

Figure  37.  Diagrams  of  compound  forces  analyzed  (Flf-Flh) ..46 

Figure  38.  Diagrams  of  compound  forces  analyzed  (Fli-Flk) . . 47 

Figure  39.  Stress  and  deflections  for  all  setups . 48 

Figure  40.  Percent  reduction  in  stiffness  variations  for 

compound  loadings . 4  8 


x 


LIST  OF  TABLES 


Table  1.  Current  NPSCuL  structure  data.  From  [11]  and 

[12] . 9 

Table  2.  Design  constraints  for  NPSCuL  system.  From 

[11]  and  [15] . 10 

Table  3.  Comparison  of  I-DEAS  modeling  and  test  results.. 15 

Table  4.  Aluminum  7075-T7351.  From  [18] . 21 

Table  5.  Comparison  of  removing  all  of  the  modifiable 

area . 23 

Table  6.  Buckling  failure  of  full  MA  removal . 25 

Table  7.  Percent  reduction  in  mass  and  stiffness 

calculations . 30 

Table  8.  Buckling  failure  for  each  configuration . 35 

Table  9.  Maximum  and  minimum  metallic  values  for 

modulus,  density  and  beam  width . 40 


xi 


THIS  PAGE  INTENTIONALLY  LEFT  BLANK 


Xll 


LIST  OF  ACRONYMS  AND  ABBREVIATIONS 


ADAMS at 

CubeSat 
Cal  Poly 
COTS 
FEM 
MA 

MAT LAB 
NPS 

NPSCuL 

NASA 

P-POD 

U 

VMA 


AS&T  Technology  Demonstration  And  Maturation 
Satellite 

Cube  Satellite 

California  Polytechnic  State  University 

Commercial  off-the-shelf 

Finite  Element  Model 

Modifiable  area 

Analysis  software 

Naval  Postgraduate  School 

Naval  Postgraduate  School  Cube  Satellite 
Launcher 

National  Aeronautics  and  Space  Administration 
Poly  Picosatellite  Orbital  Deployer 
Unit  of  10  cm3  picosatellite 
Volume  of  Modifiable  Area 


Xlll 


THIS  PAGE  INTENTIONALLY  LEFT  BLANK 


xiv 


I. 


HISTORY  AND  PURPOSE 


A.  SATELLITES 

The  space  age  kicked  off  in  1957  with  the  Russian 
launch  of  Sputnik.  Sputnik-1  was  immediately  followed  one 
month  later  with  Sputnik  II,  which  carried  a  dog  named 
Laika.  The  United  States  countered  with  the  launch  of 
Explorer  1  in  1958,  which  went  on  to  discover  the  Van  Allen 
radiation  belts.  These  feats  led  to  the  launch  of  Echo  1,  an 
aluminum-coated  plastic  sphere  that  passively  reflected 
voice  and  picture  signals.  Next  was  the  first  true 
communications  satellite,  Telstar.  Telstar  actively  relayed 
communications  between  locations  that  were  at  great 
distances  from  one  another  [1]  .  Over  the  years,  the  use  of 
satellites  has  evolved,  and  now  they  function  in  many 
different  roles.  Satellites  are  used  for  navigation,  as 
relay  stations  for  radio  and  television,  communications, 
weather,  observations  of  the  Earth  and  space,  and  for  wars. 
The  use  of  these  satellites  will  not  diminish;  rather,  as 
time  progresses,  there  will  be  new  inventive  things  that  can 
be  done  with  a  satellite.  Additionally,  as  resources  become 
scarcer,  the  satellite  will  need  to  get  smaller  and  cheaper. 

The  Galileo,  Cassini  and  Hubble  spacecrafts  had  life 
cycle  costs  over  one  billion  dollars.  The  Clementine  space 
craft  cost  $85  million  and  the  Lunar  Prospector  cost  $63 
million.  After  the  Cold  War  ended,  the  exclusivity  of  space 
to  a  select  few  countries  started  to  change.  There  became  a 
commercialization  of  space,  from  the  previously 
military/government  specific  one.  There  was  an  increase  in 
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the  availability  for  launch  services  and  a  dedicated  launch 
vehicle  was  no  longer  required.  It  became  cost  effective  to 
get  as  many  different  companies,  or  entities,  on  one  launch 
vehicle  as  possible  to  offset  the  costs.  Large  satellites, 
such  as  in  Figure  1,  can  cost  from  $5,000  to  $20,000  per 
kilogram  to  launch,  and  from  $30,000  to  $160,000  per 
kilogram  to  build  [1] .  One  driving  force  in  the  design  cost 
for  satellites  are  the  use  of  space-qualified  components. 


Figure  1.  JCSAT-12  being  built.  From  [2] 

Almost  all  companies,  including  the  U.S.  government, 
require  the  components  of  a  space  craft  or  satellite  to 
contain  space-qualified  components.  These  components  must 
meet  a  military  or  NASA  (National  Aeronautics  and  Space 
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Administration)  standard  for  the  material,  fabrication, 
storage  and  transport.  These  requirements  significantly 
increase  the  price  of  a  satellite.  Additionally,  anywhere 
from  two  to  four  models  must  be  made;  one  for  testing,  one 
for  flight,  a  possible  flight  spare  and  a  possible 
engineering  model  [1] .  Costs  need  to  be  cut  wherever 
possible.  Decreasing  the  size  requires  more  advanced 
technology  that  is  smaller  than  its  predecessors.  Satellites 
need  to  be  lighter  and  more  efficient.  The  use  of  commercial 
off-the-shelf  (COTS)  products  dramatically  reduces  these 
costs.  Multiple  spares  can  be  purchased  for  a  fraction  of 
the  price  of  designing  and  fabricating  required  circuitry. 
These  steps  have  led  to  an  increase  in  the  use  of  satellites 
by  universities  for  education  and  research  [1]. 

Satellites  are  classified  by  weight  into  size 
categories.  Cube  Satellites,  Figure  2,  which  are  known  as 
CubeSats,  can  be  nano-satellites  (weighing  from  1-10  kgs)  or 
pico-satellites  (weighing  less  than  1  kg)  and  are  generally 
10  centimeters  cubed.  They  are  orders  of  magnitude  smaller 
than  their  full  satellite  cousins.  The  idea  of  a  CubeSat  is 
not  only  important  to  the  space  industry,  but  to  education 
as  well.  CubeSats  are  becoming  important  and  useful  to  high 
schools,  colleges,  and  universities.  Students  are  designing, 
building  and  launching  these  small  objects  through  their  own 
programs .  The  knowledge  they  gain  in  these  programs  prepares 
them  for  joining  the  industry's  satellite  community.  Their 
small  size  allows  CubeSats  to  be  placed  as  a  secondary 
payload  on  an  already  scheduled  launch  vehicle,  vice  having 
the  dedicated  launch  required  for  larger  satellites.  For  the 
industry,  a  CubeSat  can  be  launched  on  a  budget  of  only 

$100,000  and  the  launch  can  take  place  years  before  a  full- 
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sized  satellite  could  [3] .  For  student  programs,  there  is  a 
discount  to  stimulate  research  and  development  for  new 
innovations,  where  a  CubeSat  can  be  launched  for  as  little 
as  $40,000  [4]  .  This  price  allows  new  and  radical 
technologies  to  be  tested  in  space  that  would  not  normally 
take  place  due  to  the  large  monetary  investment  that  the 
industry  has  on  a  given  launch.  Industry  needs  systems  that 
are  guaranteed  to  work  and  do  not  have  the  money  to  spend  on 
a  "what  if."  These  CubeSats  allow  for  a  much  lower  cost  of 
space  testing  for  these  advancements .  Launch  opportunities 
still  remain  an  issue.  Two  schools,  California  Polytechnic 
State  University  (Cal  Poly)  in  San  Luis  Obispo,  CA,  and  the 
Naval  Postgraduate  School  (NPS)  in  Monterey,  CA,  have  been 
instrumental  in  assisting  in  this  area. 


Figure  2.  Size  of  a  Cube  Satellite.  From  [5] 

B .  STUDENT  INVOLVEMENT 

Cal  Poly  has  standardized  the  industry  by  creating  the 
Poly  Picosatellite  Orbital  Deployer,  or,  as  it  is  commonly 
known,  the  P-POD.  The  P-POD  (Figure  3)  can  hold  up  to  three 
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units  (3-Us)  of  CubeSats .  A  single  unit  of  CubeSat  for  the 
P-POD  is  defined  as  a  10  cm  x  10  cm  x  10  cm  satellite.  As 
such,  a  two  unit  (2-U)  satellite  would  be  10  cm  x  10  cm  x  20 
cm.  Each  P-POD  contains  and  shelters  the  satellite (s) 
throughout  the  launch  process  and  once  orbital  altitude  is 
reached,  deploys  the  cargo.  The  P-POD  system  is  self 
contained  and  standardized.  CubeSats  going  into  the  P-PODs 
have  maximum  weight  limits  and  center  of  gravity 
constraints.  Once  CubeSats  are  built,  they  are  sent  to  Cal 
Poly,  where  they  are  incorporated  into  the  P-POD  and  sent  to 
the  launch  site  for  transportation  into  space  via  a 
previously  arranged  launch  vehicle  [6]. 


Figure  3.  Cal  Poly  Mk  II  P-POD.  From  [7] 

A  common  problem  in  satellite  launches  is  missed  launch 
opportunities.  Parts  not  being  delivered  on  time,  lack  of 
adequate  access  to  clean  rooms,  failed  inspections/tests, 
and  countless  other  reasons  also  can  cause  the  delivery  date 
of  satellites  to  be  missed.  Current  rocket  launch  procedure, 
the  method  of  CubeSat  delivery  to  space,  requires  a  hard  and 
fast  deadline  for  delivery  of  all  cargo.  This  deadline  is 
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essential  to  calculations  of  fuel  and  launch  criteria  for 
the  rocket.  If  a  CubeSat/P-POD  fails  to  be  delivered  on 
time,  weight  equivalence  is  used  and  it  is  no  longer 
possible  to  get  onto  that  flight.  The  flight  opportunity  is 
missed  and  must  be  rescheduled.  Flight  opportunities  are 
rare  and  that  is  why  the  P-POD  is  crucial.  This  allows  three 
CubeSats  to  be  launched  instead  of  one. 

C .  NPSCUL 

Currently,  there  are  about  40  to  50  university  programs 
developing  CubeSats  [4] .  All  of  these  programs  would  like  to 
test  out  their  systems  in  space  but,  even  with  Cal  Poly's 
three-unit  launcher,  there  are  still  not  enough 
opportunities.  As  such,  NPS  has  stepped  in  to  assist  with 
the  creation  of  the  NPS  CubeSat  Launcher,  or  NPSCuL  as  seen 
in  Figure  4.  NPSCuL,  in  short,  is  a  unit  that  contains  eight 
P-PODs,  increasing  a  rocket's  delivery  potential  from  3-Us 
of  CubeSats  to  24-Us.  Therefore,  when  one  launch  opportunity 
comes  along,  multiple  educational  programs  can  get  on  board. 
This  can  also  lower  the  individual  cost  per  launch. 
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Figure  4.  Image  of  NPSCuL. 


NPSCuL  has  already  been  designed,  built  and 
qualification  tested  to  ULA  ABC  launch  environment  [8]  . 
NPSCuL  was  designed  specifically  for  ESPA  (EELV  or  Evolved 
Expendable  Launch  Vehicle  secondary  payload  adaptor) .  The 
ESPA  ring  (Figure  5),  allows  for  six  payload  attachments. 
NPSCuL  can  be  connected  to  ESPA  compatible  fixtures  on  the 
Atlas,  Athena,  Delta,  Falcon  9,  Falcon  le,  Minotaur  IV, 
Taurus  and  European  Space  Agency's  VEGA  rockets  [9] .  On  the 
Atlas  V,  one  possible  union  of  the  launcher  and  rocket  is 
through  the  United  Launch  Alliance  (ULA)  Aft  Bulkhead 
Carrier  (ABC)  [10]  (Figure  6) .  The  NPSCuL  structure  consists 
of  a  5-sided  box,  8  P-Pods,  a  sequencer,  and  fasteners.  All 
specifications  for  the  current  design  of  NPSCuL,  as  tested, 
are  shown  in  Table  1 . 
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Figure  5.  ESPA  ring  on  the  EELV.  From  [13] 


Centaur  Aft 


Existing 

Mounting 

Brackets 


Aft  Bulkhead 
Carrier  (ABC) 


Secondary 
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ULA  ABC.  From  [14] 
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Table  1.  Current  NPSCuL  structure  data.  From  [11]  and  [12] 


NPSCuL 

lbs 

8  P-Pods 

116.84 

Sequencer  plus  fasteners 

8.80 

Box  plus  fasteners 

42 . 15 

Harness 

2 . 90 

Total  weight  of  NPSCuL 

170.63 

Max  allowable  weight  of  structure 

177 . 00 

Current  usable  weight 

6.37 

Even  though  there  is  a  current  design  for  the  NPSCuL,  it 
could  use  some  improvement.  While  technology  is  getting 
smaller  and  better,  it  also  allows  satellite  designers  to  put 
more  equipment  into  the  CubeSat  structure.  With  the 
advancement  in  technology  making  parts  smaller,  developers 
want  to  fit  more  into  the  CubeSat.  As  such,  this  increases 
the  demand  of  the  previously  allotted  1.33-kg  weight 
allowance  [6]  .  By  making  the  structure  lighter,  this  will 
allow  the  freed  weight  to  be  applied  back  into  the  CubeSats. 
The  design  constraints  are  shown  in  Table  2 .  The  weight 
constraint  for  launch  still  remains  with  a  maximum  weight  of 
NPSCuL,  including  the  CubeSats,  not  exceeding  177  lbs  or 
being  lighter  than  169  lbs  [11] . 
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Table  2  . 


Design  constraints  for  NPSCuL  system. 
From  [11]  and  [15] 


Max  weight  of  structure 

177 .00 

lbs 

Minimum  weight  of  structure 

169.00 

lbs 

Weight  of  a  side  of  the  box 

7.44 

lbs 

Expected  g-forces  in  X  direction  [12] 

5 

g '  s 

Expected  g-forces  in  Y  direction  [12] 

5 

g '  s 

Expected  g-forces  in  Z  direction  [12] 

7 

g '  s 

Density  of  A1  7075-T7351  [11] 

0 . 101 

lbs/ inA3 

Modulus  of  elasticity  [11] 

10.3E6 

psi 

Yield  Stress  [11] 

57 

ksi 
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II.  THE  DYNAMIC  STRUCTURE 


A.  LAUNCH  EFFECTS 

NPSCuL  must  get  to  space  to  be  effective.  The  current 
method  is  to  go  via  rocket  and  launching  on  a  rocket  is  not 
a  smooth  ride.  There  are  strong  forces,  up  to  seven  times 
the  gravitational  constant,  and  random  vibrations.  The 
natural  frequencies  of  the  structure  are  of  the  utmost 
interest.  At  these  natural  frequencies,  the  CubeSat  launcher 
can  see  a  large  amplification  of  forces,  which  in  turn  will 
be  transmitted  to  the  CubeSats.  Even  though  the  launcher 
structure  can  easily  withstand  these  forces,  the  safety  of 
the  CubeSats  inside  the  launcher  becomes  a  question.  The 
satellites,  though  built  tough  for  the  space  environment, 
still  have  sensitive  sensors  that  can  be  damaged.  The 
current  design  for  NPSCuL  has  already  been  analyzed  using  a 
Finite  Element  Model  (FEM)  in  NX  I-DEAS  6;  a  modeling 
software.  The  first  six  modes  of  the  structure  are  selected 
for  reference  in  this  research. 

B .  MODAL  SHAPES 

The  modal  shapes  refer  to  the  movement  of  the  structure 
at  the  system's  specific  natural  frequencies.  Of  NPSCuL' s 
lowest  six  modes,  this  analysis  will  only  show  modes  1,  2,  3 
and  6.  Mode  1  is  the  lowest  natural  frequency,  and  modes  2, 
3,  and  6  have  significant  center  of  gravity  movement  and 
structural  displacements.  Figures  7  through  10  show  the 
analytical  display  of  the  analysis  of  these  frequencies  [17] . 
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Figure  7.  Mode  1  analyses.  From  [17] 


Figure  8.  Mode  2  analyses.  From  [17] 
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Figure  9.  Mode  3  analyses.  From  [17] 


Figure  10.  Mode  6  analyses.  From  [17] 
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Mode  1  is  considered  "corner  flexing"  [17]  and  is  the 
lowest  natural  frequency  for  the  structure.  The  requirements 
for  the  structure  is  to  have  a  first  natural  frequency  no 
lower  than  35  Hz  [15]  .  The  first  natural  frequency  was 
analytically  evaluated  at  41.39  Hz.  Figure  7  shows  the 
rendering  of  movement  for  mode  1.  Mode  2  and  mode  3,  53.75 

and  54.73  Hz,  respectively,  are  modes  that  are  a  "rocking 
motion"  in  the  X-Y  plane,  as  seen  in  Figure  8  and  the  Y-Z 
plane,  as  seen  in  Figure  9.  Mode  6  is  a  "po-go"  motion  of 
the  entire  structure  up  and  down  onto  the  adaptor  ring.  This 
is  shown  in  Figure  10  [17]  . 

C .  VIBRATIONS  TESTING 

The  next  step  is  to  take  the  results  determined  by 
analytical  methods  and  verify  them  with  a  vibrations  test 
(Table  3)  .  The  engineering  model  of  NPSCuL  was  loaded  with 
mass  models  for  the  P-PODS  and,  with  the  assistance  of  tri- 
axial  accelerometers,  the  actual  natural  frequencies  were 
found.  Table  4  shows  the  results  of  this  test  [10]  and  Table  5 
shows  this  data  compared  to  what  was  expected  from  the 
analytical  results.  Figures  11-14  show  the  structural  movement 
of  the  first  and  primary  modes  (modes  1,  2,  3  and  6) . 
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Table  3.  Comparison  of  I-DEAS  modeling  and  test  results. 


Mode 

Analytical 

Frequency 

Tested 

Frequency 

Percent 

Error 

1 

41.39 

42 .96 

3.65 

2 

53.75 

57 . 83 

7 .06 

3 

54 .73 

59 . 32 

7 .74 

4 

84 . 68 

83 . 93 

0 .89 

5 

87 . 62 

85.25 

2 .78 

6 

95 .06 

100 . 5 

5 .41 

OCrORMATXOH:  6  6:  OCMQIMCO_100408_500r_5l4APE/42 . 05880 
MODE:  6  FREQ:  42.0580  damp:  0.5224406 

Acceleration  -  MAG  MIH:  O.OOE^OO  MAX:  4  3?E*00 
FRAME  OF  REF:  PART 


■~y 


Figure  11.  First  natural  frequency,  42.96  Hz,  structural 

deflection.  From  [10] 
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DEFORMATION:  4-4:  CCMElNED_100408_SDOF_$HAPE/57.  83442 
MODE:  4  FREQ:  57.83441  DAMP  1.987138 

Acceleration  -  MAG  MIN:  O.00E«O0  MAX:  2.52E*01 
FRAME  OF  REF  PART 


Figure  12.  "Rocking"  natural  frequency,  57.83  Hz,  structural 

deflection.  From  [10] 


DEFORMATION:  7-7: OCME I NE D_1 00408_SDOF_S NAPE /5 9 . 32238 
MODE:  7  FREQ:  59.32238  DAMP :  1.93505 

Acceleration  -  MAG  MIN:  0. C0e*00  MAX:  1  29e*01 
FRAME  OF  REF:  PART 


Figure  13. 


"Rocking"  natural  frequency,  59.32  Hz,  structural 
deflection.  From  [10] 
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OE  FORMATION :  2-2:  CCME I NE  0_1 00408.5  OOF_S  HA  PE  /1 00 . 4971 
MODE:  2  FREQ:  100.497  DAMP:  0.2834248 

Acceleration  -  MAG  MIN:  0.0OE*0O  MAX:  ?.13E*00 
FRAME  OF  KEF:  PART 


Figure  14.  "Po-go"  natural  frequency,  100.50  Hz,  structural 

deflection.  From  [10] 

D .  OBJECTIVE 

The  goal  is  to  reduce  the  mass  of  the  current 
structure.  The  NPSCuL  is  designed  to  create  a  way  to  launch 
more  CubeSats  per  launch  opportunity.  The  design  shape  was 
determined  by  fitting  the  structure  into  the  given  volume 
and  mass  constraints.  The  following  analysis  is  the  first 
step  in  optimizing  the  mass  of  the  structure  to  see  where 
weight  can  be  removed  without  unduly  reducing  its  strength. 
Once  a  reduced  weight  is  achieved,  further  analysis  to 
include  dynamic  response  can  be  made  with  the  new  design. 

There  are  two  possible  directions  for  weight  removal: 
modify  the  current  structure,  or  completely  redesign  it.  The 
first  option  of  modifying  the  structure  is  the  cheapest  and 
fastest  way  as  the  time,  materials,  and  manpower  do  not  have 
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to  be  spent  to  start  over.  An  analysis  needs  to  be  done  on 
the  weight  of  the  system  to  see  if  weight  can  be  removed. 
This  will  be  the  initial  iteration  at  looking  at 
modifications  to  the  structure.  Once  weight  can  be  removed, 
it  can  either  be  integrated  back  into  the  CubeSats'  weight 
allowance  or  remain  off.  Additional  alternatives  can  include 
external  stiffeners,  a  damping  system,  or  modifications  to 
the  adapter  ring.  All  of  these  options  eventually  go  toward 
reducing  the  total  mass  of  the  launcher  or,  in  the  case  of 
the  alternatives,  will  reduce  the  forces  seen  in  the 
CubeSats  during  launch. 

The  second  option  would  be  to  change  the  structure  from 
its  current  box  configuration  and  to  start  with  an  entirely 
new  design.  Using  a  different  shape  or  reducing  the  number 
of  CubeSats  the  structure  can  hold  are  all  possible 
alternatives,  but  are  beyond  the  scope  of  this  analysis. 
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III.  SIMPLIFICATION  OF  THE  STRUCTURE 


A.  SINGLE  WALL 


The  first  simplification  needed  is  to  take  the  four- 
walled  box  and  look  at  one  single  wall.  As  seen  in  the 
Figure  15,  a  single  wall  has  two  P-PODS  attached  onto  the 
plate  on  the  right-hand  side.  This  area  is  not  going  to  be 
modified  since  nearly  a  quarter  of  the  weight  of  the 
structure  is  concentrated  in  this  region  and  there  is  a 
large  presence  of  fasteners.  The  dots  represent  the 
locations  of  the  fasteners  for  the  P-PODS,  and  the  area  to 
be  modified  is  the  6.365-inch  by  12.0-inch  section  on  the 
left.  For  this  analysis,  the  wall  depth  remains  at  a  uniform 
0.25  inches  throughout  the  wall. 


Figure  15.  Two-dimensional  view  of  side  wall,  including  the 

modifiable  area  (MA) . 
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B.  COMPARISON  TECHNIQUE 


Modifying  the  wall  volume  of  the  current  structure, 
while  maintaining  all  remaining  dimensions,  allows  an 
analytical  comparison  to  the  original  test  data.  The 
original  structure,  of  a  solid  A1  7075  wall  plate  with 
dimensions  of  18  inches  by  13  inches  by  0.25  inches,  is  used 
as  the  baseline  when  comparing  the  total  mass  of  a  full 
wall.  When  comparing  solutions  for  the  modifiable  area  (MA) , 
the  baseline  is  an  A1  7075  wall  plate  with  dimensions  of 
6.365  inches  by  12  inches  by  0.25  inches.  Analyzing  the 
desired  alternatives  would  show  an  increase  or  decrease  in 
mass  and/or  an  increase  or  decrease  in  stiffness.  The 
following  analysis  determines  the  baseline  of  the  MA  by 
using  A1  7075. 

C.  FLAT  PLATE  ANALYSIS 

To  analyze  the  stress  and  deflection  seen  in  the  MA, 
the  gravitational  forces  must  be  turned  into  a  distributed 
load  across  the  top  of  the  wall.  Once  the  equivalent  force 
is  calculated,  the  stress  and  deflection  can  be  determined 
and  from  there  the  stiffness  can  be  calculated.  The  maximum 
weight  of  the  structure  is  177  lbs.  During  launch,  the 
structure  can  see  up  to  5  gs  of  force  in  the  X-direction,  5 
gs  of  force  in  the  Y-direction  and  7  gs  of  force  in  the  Z- 
direction.  Using  the  static-equivalent  to  these  dynamic 
forces,  the  following  calculations  show  why  1200  lbs  is 
used.  The  6.365  inch  by  12  inch  MA  is  approximately  one- 
third  of  the  total  volume  of  one  side  wall.  The  factor  of 
safety  used  for  this  calculation  is  two. 
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1  gravitational  unit  =  9.8067  N/kg  =  1.0000  lbf/lbm 


Combined  load  =  a/(52+52+72)  =  9.95  g' s  =  9.95  lbf/lbm 
9.95  lbf/lbm  *  177  lbm  =  1761.15  lbf 
1761.15  lbf  *  (1/3)  =  587.05  lbf 

587.05  lbf  *  2  S  1200  lbf 

Now,  the  full  unmodified  wall  can  be  analyzed.  The 
current  structure  is  made  from  Aluminum  7075-T7351  and  a 
full  plate  analysis  yields  the  baseline  data  for  comparison 
to  volume  removal  in  the  modifiable  area.  The  material 
properties  are  found  in  Table  4 . 

Table  4.  Aluminum  7075-T7351.  From  [18] 


Modulus  of  Elasticity  (E) 

O 
\ — 1 

X 

00 

o 
\ — 1 

psi 

Density 

0.101 

lbm/ in3 

Yield  stress  (gy) 

57 

ksi 

Ultimate  stress  (a0) 

68 

ksi 

The  mass  of  a  single  full  wall  is  equal  to 

p*t*l*h  =  0.101  lbm/in3  *  0.25  in  *  18  in  *  13  in  =  5.909  lbm 
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The  mass  of  the  modifiable  section  of  a  single  wall  equals 
p*t*l*h  =  0.101  lbm/ in3  *  0.25  in  *  6.365  in  *  12  in 

=  1.929  lbm 

Maximum  stress  is  [16] 

a  =  F/A  =  1200  lbf  /  (0.25  in  *  6.365  in)  =  754.12  psi 


The  stiffness  in  the  area  to  be  modified  [16] 
k  =  [ AE ] / 1  =  [(0.25  in  *  6.365  in)  *  10.3xl06  psi]  /  12  in  = 

1366xl03  lbf/in 


These  values  will  be  essential  in  further  calculations  and 
show  that  the  modifiable  area  is  approximately  one-third  of 
the  total  mass  of  one  side  of  a  wall.  It  is  also  important 
to  keep  in  mind  that  one  of  the  four  side  walls  contains  the 
sequencer  and  will  not  be  modified.  This  wall  is  shown  in 
Figure  16. 


</> 


-C 


© 

eg 


8.115  Inches 


k - N 

6.365  Inches 


Figure  16.  Two-dimensional  view  of  side  wall  with  sequencer 

bolt  holes. 
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IV.  MASS  REDUCTION  OF  THE  ALUMINUM  WALL 


A.  REMOVAL  OF  WALL  VOLUME 

Due  to  the  presence  of  corner  flexing  in  the  current 
structure  during  modal  movement  [10]  and  [17],  simply 
removing  the  6.365”  by  12”  section  is  unwise.  The  complete 
removal  would  further  undermine  the  corner  regions,  adding 
to  the  flexing  already  seen.  There  need  to  be  diagonal 
stiffeners  in  this  region.  The  stiffness  of  this  modifiable 
area  started  as  1366xl03  lbf/in;  however,  once  the  center 
section  is  removed,  it  becomes  approximately 

k  =  [ AE ] / 1  =  [(0.25  in  *  1  in)  *  10.3xl06  psi]  /  12  in  = 

214 . 6xl03  lbf/in 

This  reduces  the  original  stiffness  by  nearly  85%.  Table  5 
shows  the  comparison  of  a  full  wall  of  Al-7075  compared  to 
removing  the  entire  modifiable  area.  Figure  17  shows  this 
data  graphically. 

Table  5.  Comparison  of  removing  all  of  the  modifiable  area. 


Mass 

Percent  of 

mass 

Stiffness 

Percent  of 

stiffness 

Full  Wall 

5.909 

100 

1366 

100 

Removal  of  MA 

3.980 

67.4 

214.5 

15.7 
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Figure  17.  Visual  display  of  full  wall  versus  removing  the 

MA. 


This  drop  in  stiffness  is  dramatic.  The  verification  to 
ensure  that  this  stiffness  drop  does  not  undermine  the 
system  is  to  check  for  buckling  of  the  thinnest  members. 
Checking  for  buckling  along  with  verifying  that  the  yield 
stress  of  the  material  is  not  violated  and  the  deflections 
are  minimal  determines  the  suitability  of  the  alternative. 
The  sides  of  the  MA  are  0.5  inches  wide;  the  wall  has  a 
thickness  of  0.25  inches,  and  the  member  seeing  the  largest 
force  is  12  inches  long.  Using  the  buckling  criteria  [19] 
for  a  beam,  the  critical  load  that  can  be  carried  before 
failure  is  calculated.  This  value  can  then  be  compared  to 
the  maximum  force  seen  in  the  member.  The  results  of  this 
analysis,  as  shown  in  Table  6,  explain  that  the  removal  of 
the  entire  MA  results  in  failure. 

Per  =  [n2*E*I]  /L2 

=  [w*t3]/12 
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Where  I 


Table  6.  Buckling  failure  of  full  MA  removal. 


MA 

F 

(lbs) 

PCr 

(lbs) 

Status 

LO 

o 

600 

459.6 

B .  TRUSS  ANALYSIS 

The  possible  solution  to  optimizing  the  mass  is  for 
material  from  the  MA  to  be  removed  while  keeping  portions  as 
support  members.  To  do  this,  a  Finite  Element  Methods  (FEM) 
technique  is  utilized.  FEM  is  a  way  to  take  a  complicated 
problem  and  simplify  it  by  putting  it  into  a  computer 
program  to  be  analyzed.  Three-dimensional  problems  can  be 
simplified  as  two-dimensional  ones  with  the  use  of  beam 
elements.  Using  the  basic  MATLAB  code  found  in  [20]  for  a 
truss  system,  the  modifiable  area  was  turned  into  a  truss 
problem.  FEAPLYC2 . m,  FEASMBL1 . m,  FEELDOF . m  and  FETRUSS2.m 
were  all  used  as  the  basic  code  and  then  a  function  file  was 
created  to  run  these  programs  with  various  truss  designs. 
Only  the  code  in  the  function  file  was  modified  for  each 
variation  in  the  design.  The  analysis  was  conducted  on  eight 
different  truss  designs,  as  shown  in  Figures  18-20. 
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Figure  18 . 


Truss  variations  A,  B  and  C. 
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Figure  20.  Truss  variations  G  and  H. 


C .  CALCULATIONS 

Not  only  are  all  eight  configurations  analyzed,  but 
additionally  each  configuration  is  evaluated  using  three 
different  truss  beam  widths;  1  inch,  0.5  inches  and  0.25 
inches.  For  all  cases,  each  beam  remained  at  0.25  inches  of 
thickness  and  the  beam  width  refers  to  the  width  of  the 
face.  An  example  of  how  the  three  different  beam  widths 
would  look  on  configuration  A  is  shown  in  Figure  21.  The 
function  codes  for  each  variation  can  be  found  in  Appendix 
A.  Each  file  is  designated  first  by  the  letter  assigned  to 
the  configuration  as  shown  in  Figures  18-20,  and  secondly  by 
the  size  of  the  beam  width.  FEMA1  denotes  the  FEM  analysis 
of  a  1-inch  beam  width,  FEMA5  denotes  the  0.5-inch  beam 
width  and  FEMA25  denotes  the  0.25-inch  beam  width. 
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1"  width 


0.5"  width 


0.25"  width 


Figure  21.  Example  of  what  different  beam  widths  look  like  in 

the  MA. 


The  output  created  by  MATLAB  contains  the  maximum 
stress  in  each  beam  element  and  also  each  node's  deflections 
in  the  X  and  Y  directions.  Elements  and  nodes  are  numbered 
in  the  clockwise  direction.  Each  node  has  two  degrees  of 
freedom,  the  X  and  Y  directions.  As  such,  the  first  node 
would  have  the  designation  of  the  X  direction' s  deflection 
denoted  as  one,  its  Y  direction's  deflection  denoted  as  two 
and  node  two  would  have  X  and  Y  directions  listed  as  three 
and  four,  respectively.  These  results  for  all  24 
permutations  are  found  in  Appendix  B.  From  the  stress, 
deflection,  and  volume  of  removed  material,  each  variant's 
percent  of  total  wall  mass  and  stiffness  can  be  calculated. 
The  following  data  is  used  to  create  the  output  found  in 
Table  7 . 


o  =  F/A 


k  =  F/5 
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Cross-sectional  area  of  1"  width  beam  =  0.25  in2 
Cross-sectional  area  of  0.5"  width  beam  =  0.125  in2 
Cross-sectional  area  of  0.25"  width  beam  =  0.0625  in2 
Volume  of  Modifiable  Area  (VMA)  =  0.25  *  6.365  *  12 

=  19.095  in3 

Volume  of  Voids  (VV)  =  Volume  of  material  no  longer  present 

in  MA 

Table  7  shows  the  calculations  for  the  stiffness  and  mass 
loss  for  each  situation.  The  percent  decrease  in  mass  is 
plotted  against  the  percent  of  stiffness  from  the  original 
solid  wall  baseline.  Figures  22-24  show  the  results  for  beam 
width  thicknesses  of  0.25,  0.5  and  1  inch,  respectively.  The 
most  advantageous  system  would  have  the  largest  percent  of 
mass  reduction  while  having  the  largest  percent  of 
stiffness.  As  an  example,  in  Figure  22,  configuration  A  has 
the  highest  reduction  in  mass  but  the  lowest  amount  of 
stiffness.  Conversely,  configuration  C  has  the  largest 
stiffness  but  is  only  average  in  mass  reduction. 
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Table  7 


Percent  reduction  in  mass  and  stiffness  calculations 


Bum 

Width 

IWl(ln) 

8  on) 

o(U) 

o(pd) 

A 

(h*2) 

F=a*A 

(lbs) 

tef/5 

(bVln) 

% 

Stiffness 
of  full  wall 

w 

On*3) 

9fi 

redux 
of  MA 

Mare  of 

fill  wall 
(lbs) 

%  Reduction 

of  Massof 
full  wall 

Man  freed 
whole 

struct  ire  (lbs) 

A-l" 

1 

-00053 

2.4000 
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Figure  22.  Reduction  of  stiffness  versus  mass  for  0.25-inch 

beam  width. 
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Figure  23.  Reduction  of  stiffness  versus  mass  for  0.5-inch 

beam  width. 
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Figure  24.  Reduction  of  stiffness  versus  mass  for  1-inch  beam 

width . 


A  definitive  optimal  configuration  is  not  easily 
discerned.  As  mass  is  the  most  important  parameter  for  this 
analysis,  looking  at  each  configuration's  contribution 
toward  total  weight  removal  for  the  full  structure  is 
assessed.  The  computations  in  Table  7  are  used  to  achieve 
Figures  25-27.  The  best  four  choices  for  all  three  beam 
width  sizes  are  configurations  A,  B,  C  and  F.  The  next  step 
is  to  ensure  that  none  of  the  configurations  buckle  in  all 
three  width  sizes.  This  would  eliminate  a  configuration  from 
the  analysis.  Table  8  shows  the  calcuations  and  verification 
to  whether  buckling  occurs.  As  can  be  seen,  there  are  three 
cases  where  buckling  would  take  place.  Threfore, 
configuration  A  can  be  used  only  if  1-inch  beam  widths  are 
used,  and  configuration  F  can  be  used  if  0.5-  or  1-inch 
widths  are  used. 
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Figure  25.  Total  freed  mass  from  all  truss  types  at  0 

inches . 
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Figure  26.  Total  freed  mass  from  all  truss  types  at  0 

inches . 
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Figure  27 . 
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Table  8.  Buckling  failure  for  each  configuration. 


Beam 

Width  - 
W  (in) 

F=s*A 

(lbs) 

Per  (lbs) 

Status 

A  -  1" 

1 

600.00 

919.21 

Safe 

A -0.5" 

0.5 

600.00 

459.60 

A  -  0.25" 

0.25 

600.00 

229.80 

B- 1" 

1 

600.00 

3676.83 

Safe 

B  -  0.5" 

0.5 

600.00 

1838.41 

Safe 

B  -  0.25" 

0.25 

600.00 

919.21 

Safe 

C- 1" 

1 

600.00 

3676.83 

Safe 

C  -  0.5" 

0.5 

600.00 

1838.41 

Safe 

C  -  0.25" 

0.25 

600.00 

919.21 

Safe 

D- 1" 

1 

362.38 

3676.83 

Safe 

D  -  0.5" 

0.5 

362.38 

1838.41 

Safe 

D  -  0.25" 

0.25 

362.38 

919.21 

Safe 

E- 1" 

1 

600.00 

3676.83 

Safe 

E  -  0.5" 

0.5 

600.00 

1838.41 

Safe 

E  -  0.25" 

0.25 

600.00 

919.21 

Safe 

F- 1" 

1 

362.38 

919.21 

Safe 

F  -  0.5" 

0.5 

362.38 

459.60 

Safe 

F  -  0.25" 

0.25 

362.38 

229.80 

G- 1" 

1 

600.00 

3676.83 

Safe 

G  -  0.5" 

0.5 

600.00 

1838.41 

Safe 

G  -  0.25" 

0.25 

600.00 

919.21 

Safe 

H- 1" 

1 

499.60 

3676.83 

Safe 

H-  0.5" 

0.5 

499.60 

1838.41 

Safe 

H  -  0.25" 

0.25 

499.60 

919.21 

Safe 

D .  OPTIMAL  CONFIGURATION 

Now  that  the  best  configurations  based  on  weight  have 
been  established,  it  is  worthwhile  to  compare  the  remaining 
four  designs.  Figures  28-30  show  the  four  configurations 
enlarged.  Configuration  A  has  the  highest  reduction  of  mass 
but  the  lowest  stiffness.  Configuration  C  has  the  highest 
stiffness  but  is  tied  with  configuration  B  for  the  lowest 
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mass  reduction.  That  leaves  configurations  F  as  a 
tradeoff  between  mass  and  stiffness. 
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Comparison  of  four  favored  solutions  for  0.5-inch 
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V. 


CHANGE  OF  MATERIAL 


A.  CURRENT  MATERIAL 

The  current  structure  is  made  from  space-grade 
Aluminum,  A1  7075.  This  metal  is  used  due  to  its  light 
weight  and  high  resistance  to  stress  corrosion  cracking.  The 
material  is  finished  with  a  chemical  conversion  coating  for 
corrosion  protection,  and  to  lower  the  electrical  resistance 
[15]  .  A  different  approach  to  lightening  the  structure  is 
to  look  at  changing  the  metal.  An  analysis  is  to  be 
conducted  looking  at  not  only  lighter  metals  but  also 
heavier  ones  as  compared  to  A1  7075.  In  order  to  properly 
analyze  a  heavier  metal,  the  use  of  the  MA  will  be  required. 
Lighter  materials  could  incorporate  the  use  of  the  MA  but  do 
not  have  to.  Firstly,  there  are  six  pounds  of  free  weight 
currently  available  for  use  in  the  NPSCuL  configuration  when 
using  the  space-grade  Aluminum.  If  a  heavier  metal  is  used 
that  exceeds  this  weight  constraint,  voids  will  have  to  be 
made  to  maintain  the  system  within  parameters .  From  Chapter 
IV,  it  is  seen  that  for  Aluminum,  the  "F"  configuration  is 
the  most  efficient.  However,  with  a  change  in  material,  that 
configuration  may  change. 

B .  OTHER  MATERIALS 

Using  a  list  of  general,  known  and  available  metals 
[19],  a  simultaneous  assessment  can  be  made  of  different 
metals  with  different  beam  widths.  Appendix  C  shows  the  data 
used  for  the  material  analyzed.  Table  9  shows  the  range  of 
variables  input  into  a  computing  program  called  MINITAB. 
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Table  9.  Maximum  and  minimum  metallic  values  for  modulus , 

density  and  beam  width. 


Minimum 

Maximum 

Modulus  of  Elasticity  (*106  psi) 

6.5  (Magnesium) 

58  (Tungsten) 

Density  (lb/in3) 

0 . 0639  (Magnesium) 

0.697  (Gold) 

Beam  width  (in) 

0.25" 

1.0" 

C .  MINITAB 

MINITAB  is  a  systems  engineering  analysis  tool  that 
takes  a  given  set  of  inputs  and  returns  desired  outputs.  For 
this  analysis,  the  inputs  are  to  reduce  the  area  of  material 
in  the  modifiable  section,  change  the  material  to  something 
other  than  Aluminum,  or  a  combination  of  the  two.  The 
outputs  of  the  MINITAB  program  are  the  new  mass  of  the 
changed  wall  and  the  corresponding  stiffness.  These  are  the 
same  outputs  seen  from  MATLAB.  However,  an  improvement  of 
MINITAB  over  MATLAB  is  its  ability  to  analyze  three  or  more 
parameters  simultaneously.  In  essence,  a  three-dimensional 
graph  can  be  analyzed  by  holding  one  of  the  variables  and 
looking  at  how  the  other  two  variables  respond  to  the  held 
value.  This  is  the  Response  Surface  Model  (RSM) .  In  the  RSM, 
the  feasibility  solution  space  is  identified.  The  structure 
is  not  to  exceed  177  lbs.  The  output  graphs  will  show 
regions  of  acceptable  material  Modulus  of  Elasticity  and 
density.  From  here,  a  material  only  needs  to  be  selected 
that  fits  these  criteria  to  satisfy  the  constraints. 

Using  the  information  from  Table  9,  a  list  of  20  random 
combinations  is  created.  This  data  is  then  input  into  MATLAB 
using  configuration  F  program  code.  The  maximum  stress  and 
displacements  are  found  and  in  conjunction  with  the 
calculated  mass,  the  input  is  ready  for  MINITAB.  This  data 
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is  shown  in  Appendix  D.  Once  the  mass  and  stiffness  data  is 
entered,  the  program  generates  the  region  of  acceptance  in 
the  solution.  Figure  31  shows  the  constraint  of  the  data 
entered  into  MINITAB. 


Figure  31.  Parameters  used  for  MINITAB  input. 

The  outputs  of  looking  at  the  desired  0.25-,  0.5-  and  1- 
inch  beam  width  elements  are  shown  in  Figures  32-34, 
respectively.  In  Figure  32,  where  the  beam  width  is  0.25 
inches,  the  acceptability  region  is  located  in  the  white 
section  in  the  upper  left-hand  corner.  Flagged  on  the  graph 
are  approximately  where  the  metals  Titanium  (density  of  0.163 
lb/in3).  Aluminum  (density  of  0.101  lb/inA3)  and  Magnesium 
(density  of  0.0639)  would  be  located.  These  are  the  three 
lowest  density  materials  looked  at  and,  for  0.25-inch 
thickness  of  beams,  do  not  satisfy  the  mass  and  or  stiffness 
constraints.  Figure  33  has  a  larger  feasibility  region,  but 
the  materials  still  do  not  satisfy  the  constraints.  In  Figure 
34,  the  feasibility  region  shows  that  the  material  needs  to 

have  a  density  of  less  than  0.108  lbs/in3  and  a  Modulus  of 
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Elasticity  between  9.25  x  106  psi  and  40.83  x  106  psi. 
Aluminum  and  Aluminum  only,  becomes  an  acceptable  solution. 
This  hereby  reconfirms  the  design  already  established  in 
previous  chapters  that  Aluminum  7075  should  be  used  in 
Configuration  F  with  a  beam  width  of  1  inch.  The  units  of 
M_Elasticity  are  *106  psi  and  the  units  for  density  are 
lb/in3. 


Contour  Plot  of  Stiffness,  Mass 
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Figure  32.  MINITAB  results  of  holding  the  beam  width  at  0.25 

inches . 
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Contour  Plot  of  Stiffness,  Mass 
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Figure  33.  MINITAB  results  of  holding  the  beam  width 

at  0.5  inches. 
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Figure  34 .  MINITAB  results  of  holding  the  beam  width 

at  1  inch. 


43 


THIS  PAGE  INTENTIONALLY  LEFT  BLANK 


44 


VI.  COMPOUND  LOAD  ANALYSIS 


A.  ADDITIONAL  ANALYSIS 

Now  that  the  optimum  material  is  determined  to  be 
space-grade  Aluminum  and  the  optimum  configuration  is  the 
pattern  designated  by  F,  the  next  step  can  be  taken.  A 
thorough  analysis  needs  to  be  done  on  Configuration  F.  The 
first  analysis  was  conducted  with  the  gravitational  forces 
being  applied  in  the  negative  Y  direction.  However,  for  this 
two-dimensional  analysis,  the  structure  also  needs  to  be 
analyzed  for  forces  in  the  positive  and  negative  X 
directions,  the  positive  Y  direction,  and  a  combination  of 
the  two.  Figures  35-38  show  all  force  combinations  analyzed. 


6001b  6001b  6001b  6001b  6001b  6001b 


Figure  35.  Diagrams  of  compound  forces  analyzed  (Fl-Flb) . 
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Figure  36.  Diagrams  of  compound  forces  analyzed  (Flc-Fle) 
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Figure  37.  Diagrams  of  compound  forces  analyzed  (Flf-Flh) 
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Figure  38.  Diagrams  of  compound  forces  analyzed  (Fli-Flk) . 

B .  MATLAB  RESULTS 

The  compound  loads  from  Figures  35-38  are  put  into 
MATLAB  to  return  the  maximum  stress  and  deflections  for  each 
design.  Appendix  E  contains  the  MATLAB  code  used  for 
analysis,  and  Appendix  F  shows  the  stress  and  displacements 
for  each  element  and  node.  The  maximum  values  were  then  used 
in  Appendix  G  to  solve  for  how  each  combined  load  affects 
the  reduction  in  stiffness.  Figure  39  shows  each  setup's 
corresponding  stress  and  deflection.  From  this  chart,  it  can 
be  seen  that  the  largest  deflection  and  stress  is  seen  in 
setups  a,  b,  d  and  e.  Figure  40  shows  how  each  setup 
compares  in  stiffness  percentages.  The  ultimate  stress  of 
Aluminum  7075  is  68  ksi,  and  the  yield  stress  is  57  ksi.  The 
maximum  setup  stress  is  6.2  ksi,  which  is  vastly  smaller 
than  the  limits  of  the  material. 
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Figure  39.  Stress  and  deflections  for  all  setups. 
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Figure  40.  Percent  reduction  in  stiffness  variations  for 

compound  loadings . 
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VII .  CONCLUSION 


A .  CONFIGURATION 

Mass  is  the  leading  cause  of  high  cost  prices  for  the 
construction  and  launch  of  satellites.  Something  with  a  lot 
of  mass  costs  more  in  materials  and  takes  more  time  to 
integrate  all  the  systems.  This  integration  also  contributes 
a  large  portion  to  the  full  cost.  As  a  result ,  the  goal  for 
any  satellite  is  to  have  minimal  mass.  The  NPSCuL  structure 
carries  up  to  24  CubeSats  into  space.  With  the  decrease  in 
size  of  components  for  these  satellites,  the  CubeSat 
designers  are  able  to  put  more  and  more  mass  into  the  10  cm3 
volume.  These  designers  constantly  desire  more  and  more  mass 
to  be  allowed.  Optimizing  the  weight  of  NPSCuL  can  either 
allow  more  mass  to  be  put  into  the  CubeSats,  or  can  go 
toward  attempting  to  lower  the  magnification  of  forces  seen 
at  resonance  frequencies  during  launch  of  the  rocket  that 
carries  the  structure  to  space.  Looking  at  ways  to  optimize 
mass  utilization  of  modifications  to  the  modifiable  area  is 
ascertained.  Looking  at  various  truss  like  structures 
contained  inside  the  MA  to  remove  as  much  mass  as  possible, 
while  maintaining  some  stiffness,  was  completed.  Using  an 
"X"  shape  of  crossbeams  in  the  desired  area  proved  to  be  the 
best  design. 

B.  ALUMINUM  7075  AND  COMBINED  LOADING 

The  NPSCuL  structure  does  not  have  to  operate  in  space 
for  the  designed  life  that  a  satellite  does.  It  is  merely 
the  transportation  of  the  satellites  on  their  way  to  space. 
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As  such,  the  next  step  in  optimization  was  to  evaluate 
whether  Aluminum  7075  (space-grade  Aluminum)  was  the  best 
choice  for  materials.  Analyzing  the  range  of  metals  from 
Magnesium  to  Tungsten  to  Gold,  a  simultaneous  analysis  was 
conducted  using  the  "X"-shaped  configuration.  The  percent 
reduction  in  stiffness  and  the  percent  reduction  in  mass 
were  evaluated  at  different  crossbeam  widths  for  all  the 
metals  using  an  analysis  program  called  MINITAB.  The 
constraint  is  that  the  structure  cannot  exceed  177  lbs.  It 
was  found  that  the  crossbeam  width  must  be  1  inch  and  the 
material  must  be  an  Aluminum  alloy  to  meet  the  requirements. 

The  initial  analysis  used  only  forces  in  one  direction 
for  calculations  and  comparisons.  The  next  step  to  evaluate 
was  if  the  compound  loading,  which  the  structure  will  see, 
renders  the  design  a  failure.  The  yield  stress  of  Aluminum 
7075  is  57  ksi.  A  combination  of  directions  and  forces  were 
applied  and  the  maximum  stress  seen  in  any  of  the  members 
was  6.2  ksi . 

C .  FUTURE  WORK 

This  optimization  of  mass  looked  at  the  two-dimensional 
response  of  the  separate  walls.  Each  wall  only  had  six  ksi 
of  a  stress  when  analyzed  with  compound  loads  and  the  yield 
of  the  material  is  57.  The  ultimate  of  the  material  is  68 
ksi  and  stresses  seen  were  ten  times  less  than  that.  The 
next  step  in  the  process  will  be  to  take  the  modifications 
in  the  modifiable  area  and  construct  a  model  in  IDEAS.  From 
there,  a  dynamic  response  can  be  determined;  if  it  is  a 
plausible  solution,  the  structure  should  be  built  and 
vibration  tested.  The  analysis  of  nonmetal  options  for  the 

structure  should  also  be  undertaken. 
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APPENDIX  A 


FEMA1 

nel=4 ; 
nnel=2 ; 
ndof =2 ; 
nnode=4 ; 

sdof=nnode*ndof ; 

gcoord ( 1 , 1 ) =0 ; gcoord ( 1 , 2 ) =0 ; 
gcoord (2 , 1 ) =0 ; gcoord (2 , 2 ) =12 ; 
gcoord (3, 1) =  6 . 3 65; gcoord (3, 2) =12; 
gcoord (4,1)=6.365; gcoord ( 4 , 2 ) =0 ; 

elprop (1,1) =10300000; 
elprop ( 1 ,  2 ) = . 25 ; 
elprop (2 , 1)=10300000; 
elprop (2 , 2 ) = . 25 ; 
elprop (3, 1) =10300000; 
elprop (3 , 2 ) = . 25 ; 
elprop (4, 1)=10300000; 
elprop (4  r  2 )  =  . 2  5 ; 

nodes ( 1 , 1 ) =1 ; nodes ( 1 f  2 ) =2 ; 
nodes (2,1) =2 ; nodes (2, 2) =4; 
nodes ( 3 , 1 ) =2 ; nodes ( 3 , 2 ) =3 ; 
nodes (4,1) =3 ; nodes (4,2) =4 ; 

bcdof (1) =1; 
bcval ( 1 ) =0 ; 
bcdof (2) =2; 
bcval ( 2 ) =0 ; 
bcdof (3) =7; 
bcval ( 3 ) =0 ; 
bcdof (4) =8; 
bcval ( 4 ) =0 ; 

ff=zeros (sdof , 1) ; 
kk=zeros (sdof , sdof ) ; 
index=zeros (nnel*ndof ,  1)  ; 
elf orce=zeros (nnel*ndof ,  1)  ; 
eldisp=zeros (nnel*ndof , nnel^ndof ) ; 
stress=zeros (nel, 1) ; 

ff (4) =-600; 
ff (6)=-600; 

for  iel=l:nel 

nd ( 1 ) =nodes ( iel ,  1 )  ; 
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nd (2 ) =nodes ( iel ,  2 )  ; 


xl=gcoord (nd ( 1 ) ,  1 ) ; yl=gcoord (nd ( 1 ) ,  2 ) 
x2=gcoord (nd (2 ) ,  1 ) ; y2=gcoord (nd (2 ) ,2) 

leng=sqrt ( (x2-xl) A2+ (y2-yl) A2) ; 

if  (x2-xl ) ==0 ; 

beta=2*atan (1) ; 

else 

beta=atan ( (y2-yl) / (x2-xl) ) ; 

end 

el=elprop  (iel, 1) ; 
area=elprop  (iel, 2 ) ; 

[ index] =FEELDOF (nd, nnel , ndof) ; 

k=FETRUSS2 (el, leng, area, 0 , beta, 1 ) ; 

kk=FEASMBLl (kk, k, index)  ; 


end 

[ kk, f f ] =FEAPLYC2 ( kk, f f , bcdof , bcval ) ; 

disp=kk\f f ; 

for  iel=l:nel 

nd ( 1 ) =nodes  ( iel , 1 ) ; 
nd (2 ) =nodes  ( iel , 2 ) ; 

xl=gcoord (nd ( 1 ) , 1 ) ; yl=gcoord (nd  ( 1 ) , 2 ) 
x2=gcoord (nd ( 2 ) , 1 ) ; y2=gcoord (nd ( 2 )  ,  2 ) 

leng=sqrt ( (x2-xl) A2+ (y2-yl) A2) ; 

if  (x2-xl ) ==0 ; 

beta=2^atan (1) ; 

else 

beta=atan ( (y2-yl) / (x2-xl) ) ; 

end 

el=elprop  (iel, 1) ; 
area=elprop  (iel, 2 ) ; 
index=FEELDOF (nd, nnel, ndof) ; 

k=FETRUSS2 (el, leng, area, 0 , beta, 1 ) ; 

for  i=l :  (nnel^ndof) 

eldisp (i) =disp (index (i) ) ; 

end 
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elf orce=k*eldisp; 

stress (iel) =sqrt (el force (1) A2+elforce (2) A2) /area 

if ( (x2-xl) *elforce  (3) )<0; 

stress (iel) ^-stress  (iel) ; 

end 


end 

num=l : 1 : sdof ; 
displ=[num’  disp] 

numm=l : 1 : nel ; 
stresses^ [numm 1  stress] 
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FEMA5 


nel=4 ; 
nnel=2 ; 
ndof=2 ; 
nnode=4 ; 

sdof =nnode*ndof ; 

gcoord ( 1 , 1 ) =0 ; gcoord ( 1 , 2 ) =0 ; 
gcoord (2 , 1) =0 ; gcoord (2 , 2 ) =12 ; 
gcoord (3,1)=6.365; gcoord (3 , 2 ) =12 ; 
gcoord (4,1)=6.365; gcoord ( 4  ,  2 ) =0 ; 

elprop (1, 1) =10300000; 
elprop (1,2)=. 125; 
elprop (2 , 1)=10300000; 
elprop (2  ,  2 )  =  . 12  5 ; 
elprop (3, 1)=10300000; 
elprop ( 3 , 2 ) = . 125 ; 
elprop (4, 1) =10300000; 
elprop ( 4 , 2 ) = . 125 ; 

nodes ( 1 , 1 ) =1 ; nodes ( 1 , 2 ) =2 ; 
nodes (2 , 1) =2 ; nodes (2 , 2 ) =4 ; 
nodes ( 3 , 1 ) =2 ; nodes ( 3 , 2 ) =3 ; 
nodes ( 4 , 1 ) =3 ; nodes ( 4 , 2 ) =4 ; 

bcdof (1)=1; 
bcval ( 1 ) =0 ; 
bcdof (2) =2; 
bcval (2 ) =0 ; 
bcdof (3) =7; 
bcval ( 3 ) =0 ; 
bcdof (4) =8; 
bcval ( 4 ) =0 ; 

ff=zeros (sdof , 1) ; 
kk=zeros (sdof ,  sdof) ; 
index=zeros (nnel*ndof ,  1)  ; 
elf orce=zeros (nnel*ndof ,  1)  ; 
eldisp=zeros (nnel*ndof , nnel*ndof ) ; 
stress=zeros (nel, 1) ; 

ff  (4)=-600; 
ff (6) =-600; 

for  iel=l:nel 

nd (1) =nodes  (iel, 1) ; 
nd (2 ) =nodes (iel, 2 ) ; 

xl=gcoord (nd ( 1 ) , 1 ) ; yl=gcoord  (nd  ( 1 ) ,  2 ) 
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x2=gcoord (nd (2 ) ,  1 ) ; y2=gcoord (nd (2 )  ,  2 )  ; 


leng=sqrt ( (x2-xl) A2+ (y2-yl) A2) ; 

if  (x2-xl ) ==0 ; 

beta=2*atan (1) ; 

else 

beta=atan ( (y2-yl) / (x2-xl) ) ; 

end 

el=elprop  (iel, 1) ; 
area=elprop  (iel, 2 ) ; 

[ index] =FEELDOF (nd, nnel , ndof) ; 

k=FETRUSS2 (el, leng, area, 0 , beta, 1 ) ; 

kk=FEASMBLl (kk, k, index)  ; 


end 

[ kk, f f ] =FEAPLYC2 ( kk, f f , bcdof , bcval ) ; 

disp=kk\f f ; 

for  iel=l : nel 

nd (1) =nodes (iel, 1) ; 
nd (2 ) =nodes (iel, 2 ) ; 

xl=gcoord (nd ( 1 ) , 1 ) ; yl=gcoord  (nd  ( 1 ) , 2 ) ; 
x2=gcoord (nd ( 2 )  ,  1 ) ; y2=gcoord (nd ( 2 )  ,  2 )  ; 

leng=sqrt ( (x2-xl) A2+ (y2-yl) A2) / 

if  (x2-xl ) ==0 ; 

beta=2*atan (1) ; 

else 

beta=atan ( (y2-yl) / (x2-xl) ) ; 

end 

el=elprop  (iel, 1) ; 
area=elprop  (iel, 2 ) ; 
index=FEELDOF (nd, nnel , ndof ) ; 

k=FETRUSS2 (el, leng, area, 0 , beta, 1 ) ; 

for  1=1 :  (nnel*ndof) 

eldisp (i) =disp (index (i) ) ; 

end 

elf orce=k*eldisp; 

stress (iel) =sqrt (elforce (1) A2+elforce (2) A2) /area 
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if ( (x2-xl) *elforce  (3) )<0; 

stress (iel) =- stress  (iel) 

end 


end 

num=l : 1 : sdof ; 
displ=[num'  disp] 

numm=l : 1 : nel / 
stresses^ [numm 1  stress] 


FEMA25 


nel=4 ; 
nnel=2 ; 
ndof=2 ; 
nnode=4 ; 

sdof =nnode*ndof ; 

gcoord ( 1 , 1 ) =0 ; gcoord ( 1 , 2 ) =0 ; 
gcoord (2 , 1) =0 ; gcoord (2 , 2 ) =12 ; 
gcoord (3,1)=6.365; gcoord (3 , 2 ) =12 ; 
gcoord (4,1)=6.365; gcoord ( 4  ,  2 ) =0 ; 

elprop (1, 1) =10300000; 
elprop (1,2)=. 0625; 
elprop (2 , 1)=10300000; 
elprop (2,2)=. 0625; 
elprop (3, 1)=10300000; 
elprop (3,2)=. 0625; 
elprop (4, 1) =10300000; 
elprop (4,2)=. 0625; 

nodes ( 1 , 1 ) =1 ; nodes (1,2) =2; 
nodes (2 , 1) =2 ; nodes (2 , 2 ) =4 ; 
nodes ( 3 , 1 ) =2 ; nodes ( 3 , 2 ) =3 ; 
nodes ( 4 , 1 ) =3 ; nodes ( 4 , 2 ) =4 ; 

bcdof ( 1 ) =1 ; 
bcval ( 1 ) =0 ; 
bcdof (2) =2; 
bcval (2 ) =0 ; 
bcdof (3) =7; 
bcval ( 3 ) =0 ; 
bcdof (4) =8; 
bcval ( 4 ) =0 ; 

ff=zeros (sdof , 1) ; 
kk=zeros (sdof , sdof ) ; 
index=zeros (nnel*ndof ,  1)  ; 
elf orce=zeros (nnel*ndof ,  1)  ; 
eldisp=zeros (nnel*ndof , nnel*ndof ) ; 
stress=zeros (nel, 1) ; 

ff  (4)=-600; 
ff (6) =-600; 

for  iel=l:nel 

nd (1) =nodes  (iel, 1) ; 
nd (2 ) =nodes (iel, 2 ) ; 

xl=gcoord (nd ( 1 ) , 1 ) ; yl=gcoord  (nd  ( 1 ) ,2) 
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x2=gcoord (nd (2 ) ,  1 ) ; y2=gcoord (nd (2 )  ,  2 )  ; 


leng=sqrt ( (x2-xl) A2+ (y2-yl) A2) ; 

if  (x2-xl ) ==0 ; 

beta=2*atan (1) ; 

else 

beta=atan ( (y2-yl) / (x2-xl) ) ; 

end 

el=elprop  (iel, 1) ; 
area=elprop  (iel, 2 ) ; 

[ index] =FEELDOF (nd, nnel , ndof) ; 

k=FETRUSS2 (el, leng, area, 0 , beta, 1 ) ; 

kk=FEASMBLl (kk, k, index)  ; 


end 

[ kk, f f ] =FEAPLYC2 ( kk, f f , bcdof , bcval ) ; 

disp=kk\f f ; 

for  iel=l : nel 

nd (1) =nodes (iel, 1) ; 
nd (2 ) =nodes (iel, 2 ) ; 

xl=gcoord (nd ( 1 ) , 1 ) ; yl=gcoord  (nd  ( 1 ) , 2 ) ; 
x2=gcoord (nd ( 2 )  ,  1 ) ; y2=gcoord (nd ( 2 )  ,  2 )  ; 

leng=sqrt ( (x2-xl) A2+ (y2-yl) A2) / 

if  (x2-xl ) ==0 ; 

beta=2*atan (1) ; 

else 

beta=atan ( (y2-yl) / (x2-xl) ) ; 

end 

el=elprop  (iel, 1) ; 
area=elprop  (iel, 2 ) ; 
index=FEELDOF (nd, nnel , ndof ) ; 

k=FETRUSS2 (el, leng, area, 0 , beta, 1 ) ; 

for  1=1 :  (nnel*ndof) 

eldisp (i) =disp (index (i) ) ; 

end 

elf orce=k*eldisp; 

stress (iel) =sqrt (elforce (1) A2+elforce (2) A2) /area 
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if ( (x2-xl) *elforce  (3) )<0; 

stress (iel) =- stress  (iel) 

end 


end 

num=l : 1 : sdof ; 
displ=[num'  disp] 

numm=l : 1 : nel / 
stresses^ [numm 1  stress] 


FEMB1 


nel=8 ; 
nnel=2 ; 
ndof=2 ; 
nnode=6 ; 

sdof =nnode*ndof ; 

gcoord ( 1 , 1 ) =0 ; gcoord ( 1 , 2 ) =0 ; 
gcoord  (2, ‘1)  =0;  gcoord  (2 , 2)  =  6; 
gcoord (3,1) =0; gcoord (3, 2) =  12; 
gcoord (4,1)=6.365; gcoord (4 , 2 ) =12 ; 
gcoord (5,1)=6.365; gcoord (5, 2) =6; 
gcoord (6, 1 ) =  6 . 3  6  5 ; gcoord ( 6, 2 ) =0 ; 

elprop ( 1 , 1 ) =10300000 ; 
elprop ( 1 , 2 )  =  . 25 ; 
elprop (2, 1)=10300000; 
elprop (2 , 2 ) = . 25 ; 
elprop (3, 1) =10300000 ; 
elprop  (3, 2)  =  . 25; 
elprop (4, 1)=10300000; 
elprop  (4 , 2 )  =  . 25 ; 
elprop (5, 1) =10300000; 
elprop (5, 2) =.25; 
elprop  (6,  1)=10300000; 
elprop ( 6, 2 ) = . 25 ; 
elprop (7, 1) =10300000; 
elprop  ( 7 , 2 )  =  . 2  5 ; 
elprop (8, 1) =10300000; 
elprop ( 8 , 2 ) = . 25 ; 

nodes ( 1 , 1 ) =1 ; nodes ( 1 , 2 ) =2 ; 
nodes (2,1) =2; nodes (2, 2) =6; 
nodes (3,1) =2; nodes (3, 2) =5; 
nodes (4,1) =2 ; nodes (4,2) =3 ; 
nodes (5, 1) =3; nodes (5, 2) =5; 
nodes (6, 1) =3; nodes ( 6, 2 ) =4 ; 
nodes (7, 1) =4; nodes (7, 2) =5; 
nodes (8,1) =5 ; nodes (8, 2) =6; 

bcdof ( 1 ) =1 ; 
bcval ( 1 ) =0 ; 
bcdof (2) =2; 
bcval (2 ) =0 ; 
bcdof (3) =41; 
bcval (3) =0 ; 
bcdof (4)=12; 
bcval ( 4 ) =0 ; 

ff=zeros (sdof , 1) ; 
kk=zeros (sdof ,  sdof )  ; 
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index=zeros (nnel*ndof ,  1)  ; 
elf orce=zeros (nnel*ndof ,  1)  ; 
eldisp=zeros (nnel*ndof , nnel*ndof ) ; 
stress=zeros (nel, 1) ; 

ff (6) =-600; 
ff  (8)—  600; 

for  iel=l : nel 

nd (1) =nodes (iel,  1)  ; 
nd (2 ) =nodes (iel, 2)  ; 

xl=gcoord (nd ( 1 ) ,1) ; yl=gcoord (nd ( 1 )  ,  2 ) 
x2=gcoord (nd (2 ) ,1) ; y2=gcoord (nd (2 ) f2) 

leng=sqrt ( (x2-xl) A2+ (y2-yl) A2) ; 

if  (x2-xl ) ==0 ; 

beta=2*atan (1) ; 

else 

beta=atan ( (y2-yl) / (x2-xl) ) ; 

end 

el=elprop  (iel, 1) ; 
area=elprop  (iel , 2 ) ; 

[ index] =FEELDOF (nd, nnel , ndof) ; 

k=FETRUSS2 (el, leng, area, 0, beta, 1) ; 

kk=FEASMBLl (kk, k, index) ; 


end 

[ kk, f f ] =FEAPLYC2 ( kk, f f , bcdof , bcval ) ; 

disp=kk\f f ; 

for  iel**!  :  nel 

nd ( 1 ) =nodes ( iel ,  1 )  ; 
nd (2 ) =nodes (iel,  2 )  ; 

xl=gcoord (nd ( 1 ) , 1 ) ; yl=gcoord (nd ( 1 ) , 2 ) 
x2=gcoord (nd (2 ) , 1 ) ; y2=gcoord (nd (2 )  ,  2 ) 

leng=sqrt ( (x2-xl) A2+ (y2-yl) A2) ; 

if  (x2-xl ) ==0 ; 

beta=2*atan (1) ; 

else 

beta=atan ( (y2-yl) / (x2-xl) ) ; 
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end 


el=elprop  (iel, 1) ; 
area=elprop  (iel, 2 ) ; 
index=FEELDOF (nd, nnel , ndof ) ; 

k=FETRUSS2 (el , leng, area, 0 , beta, 1 ) ; 

for  i=l :  (nnel*ndof) 

eldisp (i) =disp (index (i) ) ; 

end 

elf orce=k* eldisp; 

stress (iel) =sqrt (el force (1) A2+elforce (2) A2) /area 

if ( (x2-xl) *el force  (3) )  <0 ; 

stress (iel) ^-stress  (iel) ; 

end 


end 

num=l : 1 : sdof ; 
displ=[num'  disp] 

numm=l : 1 : nel / 
stresses^ [numm 1  stress] 
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FEMB5 


nel=8 ; 
nnel=2 ; 
ndof=2 ; 
nnode=6 ; 

sdof =nnode*ndof ; 

gcoord ( 1 , 1 ) =0 ; gcoord ( 1 , 2 ) =0 ; 
gcoord  (2, ‘1)  =0;  gcoord  (2 , 2)  =  6; 
gcoord (3,1) =0; gcoord (3, 2) =  12; 
gcoord (4,1)=6.365; gcoord (4 , 2 ) =12 ; 
gcoord (5,1)=6.365; gcoord (5, 2) =6; 
gcoord (6, 1 ) =  6 . 3  6  5 ; gcoord ( 6, 2 ) =0 ; 

elprop (1, 1)=10300000; 
elprop ( 1 , 2 ) = . 125 ; 
elprop (2, 1)=10300000; 
elprop (2 , 2 ) = . 125 ; 
elprop (3, 1) =10300000 ; 
elprop  (3, 2)  =  . 12  5; 
elprop (4, 1)=10300000; 
elprop (4 , 2 ) = . 125 ; 
elprop (5, 1) =10300000; 
elprop (5,2)=. 125; 
elprop  (6,  1)=10300000; 
elprop ( 6, 2 ) = . 125 ; 
elprop (7, 1) =10300000; 
elprop  ( 7 , 2 )  =  . 12  5 ; 
elprop (8, 1) =10300000; 
elprop  ( 8 , 2 )  =  . 12  5 ; 

nodes ( 1 , 1 ) =1 ; nodes ( 1 , 2 ) =2 ; 
nodes (2,1) =2; nodes (2, 2) =6; 
nodes (3,1) =2; nodes (3, 2) =5; 
nodes (4,1) =2 ; nodes (4,2) =3 ; 
nodes (5, 1) =3; nodes (5, 2) =5; 
nodes (6, 1) =3; nodes ( 6, 2 ) =4 ; 
nodes (7, 1) =4; nodes (7, 2) =5; 
nodes (8,1) =5 ; nodes (8, 2) =6; 

bcdof ( 1 ) =1 ; 
bcval ( 1 ) =0 ; 
bcdof (2) =2; 
bcval (2 ) =0 ; 
bcdof (3) =41; 
bcval (3) =0 ; 
bcdof (4)=12; 
bcval ( 4 ) =0 ; 

ff=zeros (sdof , 1) ; 
kk=zeros (sdof ,  sdof )  ; 
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index=zeros (nnel*ndof ,  1)  ; 
elf orce=zeros (nnel*ndof ,  1)  ; 
eldisp=zeros (nnel*ndof , nnel*ndof ) ; 
stress=zeros (nel, 1) ; 

ff (6) =-600; 
ff  (8)—  600; 

for  iel=l : nel 

nd (1) =nodes (iel,  1)  ; 
nd (2 ) =nodes (iel, 2)  ; 

xl=gcoord (nd ( 1 ) ,1) ; yl=gcoord (nd ( 1 )  ,  2 ) 
x2=gcoord (nd (2 ) ,1) ; y2=gcoord (nd (2 ) f2) 

leng=sqrt ( (x2-xl) A2+ (y2-yl) A2) ; 

if  (x2-xl ) ==0 ; 

beta=2*atan (1) ; 

else 

beta=atan ( (y2-yl) / (x2-xl) ) ; 

end 

el=elprop  (iel, 1) ; 
area=elprop  (iel , 2 ) ; 

[ index] =FEELDOF (nd, nnel , ndof) ; 

k=FETRUSS2 (el, leng, area, 0, beta, 1) ; 

kk=FEASMBLl (kk, k, index) ; 


end 

[ kk, f f ] =FEAPLYC2 ( kk, f f , bcdof , bcval ) ; 

disp=kk\f f ; 

for  iel**!  :  nel 

nd ( 1 ) =nodes ( iel ,  1 )  ; 
nd (2 ) =nodes (iel,  2 )  ; 

xl=gcoord (nd ( 1 ) , 1 ) ; yl=gcoord (nd ( 1 ) , 2 ) 
x2=gcoord (nd (2 ) , 1 ) ; y2=gcoord (nd (2 )  ,  2 ) 

leng=sqrt ( (x2-xl) A2+ (y2-yl) A2) ; 

if  (x2-xl ) ==0 ; 

beta=2*atan (1) ; 

else 

beta=atan ( (y2-yl) / (x2-xl) ) ; 
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end 


el=elprop  (iel, 1) ; 
area=elprop  (iel, 2 ) ; 
index=FEELDOF (nd, nnel , ndof ) ; 

k=FETRUSS2 (el , leng, area, 0 , beta, 1 ) ; 

for  i=l :  (nnel*ndof) 

eldisp (i) =disp (index (i) ) ; 

end 

elf orce=k* eldisp; 

stress (iel) =sqrt (el force (1) A2+elforce (2) A2) /area 

if ( (x2-xl) *el force  (3) )  <0 ; 

stress (iel) ^-stress  (iel) ; 

end 


end 

num=l : 1 : sdof ; 
displ=[num'  disp] 

numm=l : 1 : nel / 
stresses^ [numm 1  stress] 
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FEMB25 


nel=8 ; 
nnel=2 ; 
ndof=2 ; 
nnode=6 ; 

sdof =nnode*ndof ; 

gcoord ( 1 , 1 ) =0 ; gcoord ( 1 , 2 ) =0 ; 
gcoord  (2, ‘1)  =0;  gcoord  (2 , 2)  =  6; 
gcoord (3,1) =0; gcoord (3, 2) =  12; 
gcoord (4,1)=6.365; gcoord (4 , 2 ) =12 ; 
gcoord (5,1)=6.365; gcoord (5, 2) =6; 
gcoord (6, 1 ) =  6 . 3  6  5 ; gcoord ( 6, 2 ) =0 ; 

elprop (1, 1)=10300000; 
elprop (1,2)=. 0625; 
elprop (2, 1)=10300000; 
elprop (2,2)=.0625; 
elprop (3, 1) =10300000 ; 
elprop  (3,2)  =  .0625; 
elprop (4, 1)=10300000; 
elprop  (4,2)  =  .0625; 
elprop (5, 1) =10300000; 
elprop  (5,2)  =  . 0625 ; 
elprop  (6,  1)=10300000; 
elprop (6,2)=.0625; 
elprop (7, 1) =10300000; 
elprop  (7, 2)  =  . 0625; 
elprop (8, 1) =10300000; 
elprop  (8, 2)  =  . 0625; 

nodes (1, 1) =1; nodes (1,2) =2 ; 
nodes (2,1) =2; nodes (2, 2) =6; 
nodes (3,1) =2; nodes (3, 2) =5; 
nodes (4,1) =2 ; nodes (4,2) =3 ; 
nodes (5, 1) =3; nodes (5, 2) =5; 
nodes (6, 1) =3; nodes ( 6, 2 ) =4 ; 
nodes (7, 1) =4; nodes (7, 2) =5; 
nodes (8,1) =5 ; nodes (8, 2) =6; 

bcdof ( 1 ) =1 ; 
bcval ( 1 ) =0 ; 
bcdof (2) =2; 
bcval (2 ) =0 ; 
bcdof (3) =41; 
bcval (3) =0 ; 
bcdof (4)=12; 
bcval ( 4 ) =0 ; 

ff=zeros (sdof , 1) ; 
kk=zeros (sdof ,  sdof )  ; 
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index=zeros (nnel*ndof ,  1)  ; 
elf orce=zeros (nnel*ndof ,  1)  ; 
eldisp=zeros (nnel*ndof , nnel*ndof ) ; 
stress=zeros (nel, 1) ; 

ff (6) =-600; 
ff  (8)—  600; 

for  iel=l : nel 

nd (1) =nodes (iel,  1)  ; 
nd (2 ) =nodes (iel, 2)  ; 

xl=gcoord (nd ( 1 ) ,1) ; yl=gcoord (nd ( 1 )  ,  2 ) 
x2=gcoord (nd (2 ) ,1) ; y2=gcoord (nd (2 ) f2) 

leng=sqrt ( (x2-xl) A2+ (y2-yl) A2) ; 

if  (x2-xl ) ==0 ; 

beta=2*atan (1) ; 

else 

beta=atan ( (y2-yl) / (x2-xl) ) ; 

end 

el=elprop  (iel, 1) ; 
area=elprop  (iel , 2 ) ; 

[ index] =FEELDOF (nd, nnel , ndof) ; 

k=FETRUSS2 (el, leng, area, 0, beta, 1) ; 

kk=FEASMBLl (kk, k, index) ; 


end 

[ kk, f f ] =FEAPLYC2 ( kk, f f , bcdof , bcval ) ; 

disp=kk\f f ; 

for  iel**!  :  nel 

nd ( 1 ) =nodes ( iel ,  1 )  ; 
nd (2 ) =nodes (iel,  2 )  ; 

xl=gcoord (nd ( 1 ) , 1 ) ; yl=gcoord (nd ( 1 ) , 2 ) 
x2=gcoord (nd (2 ) , 1 ) ; y2=gcoord (nd (2 )  ,  2 ) 

leng=sqrt ( (x2-xl) A2+ (y2-yl) A2) ; 

if  (x2-xl ) ==0 ; 

beta=2*atan (1) ; 

else 

beta=atan ( (y2-yl) / (x2-xl) ) ; 
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end 


el=elprop  (iel, 1) ; 
area=elprop  (iel, 2 ) ; 
index=FEELDOF (nd, nnel , ndof ) ; 

k=FETRUSS2 (el , leng, area, 0 , beta, 1 ) ; 

for  i=l :  (nnel*ndof) 

eldisp (i) =disp (index (i) ) ; 

end 

elf orce=k* eldisp; 

stress (iel) =sqrt (el force (1) A2+elforce (2) A2) /area 

if ( (x2-xl) *el force  (3) )  <0 ; 

stress (iel) ^-stress  (iel) ; 

end 


end 

num=l : 1 : sdof ; 
displ=[num'  disp] 

numm=l : 1 : nel / 
stresses^ [numm 1  stress] 
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FEMC1 


nel=8 ; 
nnel=2 ; 
ndof=2 ; 
nnode=6 ; 

sdof =nnode*ndof ; 

gcoord ( 1 , 1 ) =0 ; gcoord ( 1 , 2 ) =0 ; 
gcoord  (2, ‘1)  =0;  gcoord  (2 , 2)  =  6; 
gcoord (3,1) =0; gcoord (3, 2) =  12; 
gcoord (4,1)=6.365; gcoord (4 , 2 ) =12 ; 
gcoord (5,1)=6.365; gcoord (5, 2) =6; 
gcoord (6, 1 ) =  6 . 3  6  5 ; gcoord ( 6, 2 ) =0 ; 

elprop ( 1 , 1 ) =10300000 ; 
elprop ( 1 , 2 )  =  . 25 ; 
elprop (2, 1)=10300000; 
elprop (2 , 2 ) = . 25 ; 
elprop (3, 1) =10300000 ; 
elprop  (3, 2)  =  . 25; 
elprop (4, 1)=10300000; 
elprop  (4 , 2 )  =  . 25 ; 
elprop (5, 1) =10300000; 
elprop (5, 2) =.25; 
elprop  (6,  1)=10300000; 
elprop ( 6, 2 ) = . 25 ; 
elprop (7, 1) =10300000; 
elprop  ( 7 , 2 )  =  . 2  5 ; 
elprop (8, 1) =10300000; 
elprop ( 8 , 2 ) = . 25 ; 

nodes ( 1 , 1 ) =1 ; nodes ( 1 , 2 ) =5 ; 
nodes (2 , 1 ) =1 ; nodes (2 , 2 ) =2 ; 
nodes (3,1) =2; nodes (3, 2) =5; 
nodes (4,1) =2 ; nodes (4,2) =3 ; 
nodes (5, 1) =3; nodes (5, 2) =5; 
nodes (6, 1) =3; nodes ( 6, 2 ) =4 ; 
nodes (7, 1) =4; nodes (7, 2) =5; 
nodes (8,1) =5 ; nodes (8, 2) =6; 

bcdof ( 1 ) =1 ; 
bcval ( 1 ) =0 ; 
bcdof (2) =2; 
bcval (2 ) =0 ; 
bcdof (3) =41; 
bcval (3) =0 ; 
bcdof (4)=12; 
bcval ( 4 ) =0 ; 

ff=zeros (sdof , 1) ; 
kk=zeros (sdof ,  sdof )  ; 
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index=zeros (nnel*ndof ,  1)  ; 
elf orce=zeros (nnel*ndof ,  1)  ; 
eldisp=zeros (nnel*ndof , nnel*ndof ) ; 
stress=zeros (nel, 1) ; 

ff (6) =-600; 
ff  (8)—  600; 

for  iel=l : nel 

nd (1) =nodes (iel,  1)  ; 
nd (2 ) =nodes (iel, 2)  ; 

xl=gcoord (nd ( 1 ) ,1) ; yl=gcoord (nd ( 1 )  ,  2 ) 
x2=gcoord (nd (2 ) ,1) ; y2=gcoord (nd (2 ) f2) 

leng=sqrt ( (x2-xl) A2+ (y2-yl) A2) ; 

if  (x2-xl ) ==0 ; 

beta=2*atan (1) ; 

else 

beta=atan ( (y2-yl) / (x2-xl) ) ; 

end 

el=elprop  (iel, 1) ; 
area=elprop  (iel , 2 ) ; 

[ index] =FEELDOF (nd, nnel , ndof) ; 

k=FETRUSS2 (el, leng, area, 0, beta, 1) ; 

kk=FEASMBLl (kk, k, index) ; 


end 

[ kk, f f ] =FEAPLYC2 ( kk, f f , bcdof , bcval ) ; 

disp=kk\f f ; 

for  iel**!  :  nel 

nd ( 1 ) =nodes ( iel ,  1 )  ; 
nd (2 ) =nodes (iel,  2 )  ; 

xl=gcoord (nd ( 1 ) , 1 ) ; yl=gcoord (nd ( 1 ) , 2 ) 
x2=gcoord (nd (2 ) , 1 ) ; y2=gcoord (nd (2 )  ,  2 ) 

leng=sqrt ( (x2-xl) A2+ (y2-yl) A2) ; 

if  (x2-xl ) ==0 ; 

beta=2*atan (1) ; 

else 

beta=atan ( (y2-yl) / (x2-xl) ) ; 
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end 


el=elprop  (iel, 1) ; 
area=elprop  (iel, 2 ) ; 
index=FEELDOF (nd, nnel , ndof ) ; 

k=FETRUSS2 (el , leng, area, 0 , beta, 1 ) ; 

for  i=l :  (nnel*ndof) 

eldisp (i) =disp (index (i) ) ; 

end 

elf orce=k* eldisp; 

stress (iel) =sqrt (el force (1) A2+elforce (2) A2) /area 

if ( (x2-xl) *el force  (3) )  <0 ; 

stress (iel) ^-stress  (iel) ; 

end 


end 

num=l : 1 : sdof ; 
displ=[num'  disp] 

numm=l : 1 : nel / 
stresses^ [numm 1  stress] 
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FEMC5 


nel=8 ; 
nnel=2 ; 
ndof=2 ; 
nnode=6 ; 

sdof =nnode*ndof ; 

gcoord ( 1 , 1 ) =0 ; gcoord ( 1 , 2 ) =0 ; 
gcoord  (2, ‘1)  =0;  gcoord  (2 , 2)  =  6; 
gcoord (3,1) =0; gcoord (3, 2) =  12; 
gcoord (4,1)=6.365; gcoord (4 , 2 ) =12 ; 
gcoord (5,1)=6.365; gcoord (5, 2) =6; 
gcoord (6, 1 ) =  6 . 3  6  5 ; gcoord ( 6, 2 ) =0 ; 

elprop (1, 1)=10300000; 
elprop ( 1 , 2 ) = . 125 ; 
elprop (2, 1)=10300000; 
elprop (2 , 2 ) = . 125 ; 
elprop (3, 1) =10300000 ; 
elprop  (3, 2)  =  . 12  5; 
elprop (4, 1)=10300000; 
elprop (4 , 2 ) = . 125 ; 
elprop (5, 1) =10300000; 
elprop (5,2)=. 125; 
elprop  (6,  1)=10300000; 
elprop ( 6, 2 ) = . 125 ; 
elprop (7, 1) =10300000; 
elprop  ( 7 , 2 )  =  . 12  5 ; 
elprop (8, 1) =10300000; 
elprop  ( 8 , 2 )  =  . 12  5 ; 

nodes ( 1 , 1 ) =1 ; nodes ( 1 , 2 ) =5 ; 
nodes (2,1) Ml ; nodes (2,2) =2 ; 
nodes (3,1) =2; nodes (3, 2) =5; 
nodes (4,1) =2 ; nodes (4,2) =3 ; 
nodes (5, 1) =3; nodes (5, 2) =5; 
nodes (6, 1) =3; nodes ( 6, 2 ) =4 ; 
nodes (7, 1) =4; nodes (7, 2) =5; 
nodes (8,1) =5 ; nodes (8, 2) =6; 

bcdof ( 1 ) =1 ; 
bcval ( 1 ) =0 ; 
bcdof (2) =2; 
bcval (2 ) =0 ; 
bcdof (3) =41; 
bcval (3) =0 ; 
bcdof (4)=12; 
bcval ( 4 ) =0 ; 

ff=zeros (sdof , 1) ; 
kk=zeros (sdof ,  sdof )  ; 
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index=zeros (nnel*ndof ,  1)  ; 
elf orce=zeros (nnel*ndof ,  1)  ; 
eldisp=zeros (nnel*ndof , nnel*ndof ) ; 
stress=zeros (nel, 1) ; 

ff (6) =-600; 
ff  (8)—  600; 

for  iel=l : nel 

nd (1) =nodes (iel,  1)  ; 
nd (2 ) =nodes (iel, 2)  ; 

xl=gcoord (nd ( 1 ) ,1) ; yl=gcoord (nd ( 1 )  ,  2 ) 
x2=gcoord (nd (2 ) ,1) ; y2=gcoord (nd (2 ) f2) 

leng=sqrt ( (x2-xl) A2+ (y2-yl) A2) ; 

if  (x2-xl ) ==0 ; 

beta=2*atan (1) ; 

else 

beta=atan ( (y2-yl) / (x2-xl) ) ; 

end 

el=elprop  (iel, 1) ; 
area=elprop  (iel , 2 ) ; 

[ index] =FEELDOF (nd, nnel , ndof) ; 

k=FETRUSS2 (el, leng, area, 0, beta, 1) ; 

kk=FEASMBLl (kk, k, index) ; 


end 

[ kk, f f ] =FEAPLYC2 ( kk, f f , bcdof , bcval ) ; 

disp=kk\f f ; 

for  iel**!  :  nel 

nd ( 1 ) =nodes ( iel ,  1 )  ; 
nd (2 ) =nodes (iel,  2 )  ; 

xl=gcoord (nd ( 1 ) , 1 ) ; yl=gcoord (nd ( 1 ) , 2 ) 
x2=gcoord (nd (2 ) , 1 ) ; y2=gcoord (nd (2 )  ,  2 ) 

leng=sqrt ( (x2-xl) A2+ (y2-yl) A2) ; 

if  (x2-xl ) ==0 ; 

beta=2*atan (1) ; 

else 

beta=atan ( (y2-yl) / (x2-xl) ) ; 
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end 


el=elprop  (iel, 1) ; 
area=elprop  (iel, 2 ) ; 
index=FEELDOF (nd, nnel , ndof ) ; 

k=FETRUSS2 (el , leng, area, 0 , beta, 1 ) ; 

for  i=l :  (nnel*ndof) 

eldisp (i) =disp (index (i) ) ; 

end 

elf orce=k* eldisp; 

stress (iel) =sqrt (el force (1) A2+elforce (2) A2) /area 

if ( (x2-xl) *el force  (3) )  <0 ; 

stress (iel) ^-stress  (iel) ; 

end 


end 

num=l : 1 : sdof ; 
displ=[num'  disp] 

numm=l : 1 : nel / 
stresses^ [numm 1  stress] 
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FEMC25 


nel=8 ; 
nnel=2 ; 
ndof=2 ; 
nnode=6 ; 

sdof =nnode*ndof ; 

gcoord ( 1 , 1 ) =0 ; gcoord ( 1 , 2 ) =0 ; 
gcoord  (2, ‘1)  =0;  gcoord  (2 , 2)  =  6; 
gcoord (3,1) =0; gcoord (3, 2) =  12; 
gcoord (4,1)=6.365; gcoord (4 , 2 ) =12 ; 
gcoord (5,1)=6.365; gcoord (5, 2) =6; 
gcoord (6, 1 ) =  6 . 3  6  5 ; gcoord ( 6, 2 ) =0 ; 

elprop (1, 1)=10300000; 
elprop (1,2)=. 0625; 
elprop (2, 1)=10300000; 
elprop (2,2)=.0625; 
elprop (3, 1) =10300000 ; 
elprop  (3,2)  =  .0625; 
elprop (4, 1)=10300000; 
elprop  (4,2)  =  .0625; 
elprop (5, 1) =10300000; 
elprop  (5,2)  =  . 0625 ; 
elprop  (6,  1)=10300000; 
elprop (6,2)=.0625; 
elprop (7, 1) =10300000; 
elprop  (7, 2)  =  . 0625; 
elprop (8, 1) =10300000; 
elprop  (8, 2)  =  . 0625; 

nodes (1, 1) =1; nodes (1,2) =5 ; 
nodes (2, 1) =1; nodes (2,2) =2 ; 
nodes (3,1) =2; nodes (3, 2) =5; 
nodes (4,1) =2 ; nodes (4,2) =3 ; 
nodes (5, 1) =3; nodes (5, 2) =5; 
nodes (6, 1) =3; nodes ( 6, 2 ) =4 ; 
nodes (7, 1) =4; nodes (7, 2) =5; 
nodes (8,1) =5 ; nodes (8, 2) =6; 

bcdof ( 1 ) =1 ; 
bcval ( 1 ) =0 ; 
bcdof (2) =2; 
bcval (2 ) =0 ; 
bcdof (3) =41; 
bcval (3) =0 ; 
bcdof (4)=12; 
bcval ( 4 ) =0 ; 

ff=zeros (sdof , 1) ; 
kk=zeros (sdof ,  sdof )  ; 
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index=zeros (nnel*ndof ,  1)  ; 
elf orce=zeros (nnel*ndof ,  1)  ; 
eldisp=zeros (nnel*ndof , nnel*ndof ) ; 
stress=zeros (nel, 1) ; 

ff (6) =-600; 
ff  (8)—  600; 

for  iel=l : nel 

nd (1) =nodes (iel,  1)  ; 
nd (2 ) =nodes (iel, 2)  ; 

xl=gcoord (nd ( 1 ) ,1) ; yl=gcoord (nd ( 1 )  ,  2 ) 
x2=gcoord (nd (2 ) ,1) ; y2=gcoord (nd (2 ) f2) 

leng=sqrt ( (x2-xl) A2+ (y2-yl) A2) ; 

if  (x2-xl ) ==0 ; 

beta=2*atan (1) ; 

else 

beta=atan ( (y2-yl) / (x2-xl) ) ; 

end 

el=elprop  (iel, 1) ; 
area=elprop  (iel , 2 ) ; 

[ index] =FEELDOF (nd, nnel , ndof) ; 

k=FETRUSS2 (el, leng, area, 0, beta, 1) ; 

kk=FEASMBLl (kk, k, index) ; 


end 

[ kk, f f ] =FEAPLYC2 ( kk, f f , bcdof , bcval ) ; 

disp=kk\f f ; 

for  iel**!  :  nel 

nd ( 1 ) =nodes ( iel ,  1 )  ; 
nd (2 ) =nodes (iel,  2 )  ; 

xl=gcoord (nd ( 1 ) , 1 ) ; yl=gcoord (nd ( 1 ) , 2 ) 
x2=gcoord (nd (2 ) , 1 ) ; y2=gcoord (nd (2 )  ,  2 ) 

leng=sqrt ( (x2-xl) A2+ (y2-yl) A2) ; 

if  (x2-xl ) ==0 ; 

beta=2*atan (1) ; 

else 

beta=atan ( (y2-yl) / (x2-xl) ) ; 
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end 


el=elprop  (iel, 1) ; 
area=elprop  (iel, 2 ) ; 
index=FEELDOF (nd, nnel , ndof ) ; 

k=FETRUSS2 (el , leng, area, 0 , beta, 1 ) ; 

for  i=l :  (nnel*ndof) 

eldisp (i) =disp (index (i) ) ; 

end 

elf orce=k* eldisp; 

stress (iel) =sqrt (el force (1) A2+elforce (2) A2) /area 

if ( (x2-xl) *el force  (3) )  <0 ; 

stress (iel) ^-stress  (iel) ; 

end 


end 

num=l : 1 : sdof ; 
displ=[num'  disp] 

numm=l : 1 : nel / 
stresses^ [numm 1  stress] 
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FEMD1 


nel=ll ; 
nnel=2 ; 
ndof=2 ; 
nnode=7 ; 

sdof =nnode*ndof ; 

gcoord ( 1 , 1 ) =0 ; gcoord ( 1 , 2 ) =0 ; 
gcoord (2 , 1) =0; gcoord (2 , 2) =  6; 
gcoord (3,1) =0; gcoord (3, 2) =  12; 
gcoord (4,1)=6.365; gcoord (4 , 2 ) =12 ; 
gcoord (5,1)=6.365; gcoord (5, 2) =  6; 
gcoord (6, i ) =  6 . 3  6  5 ; gcoord ( 6, 2 ) =  0 ; 
gcoord (7 , 1) =3 . 1825; gcoord (7 , 2) =  6; 

elprop (1, 1)=10300000; 
elprop ( 1 ,  2 ) = . 25 ; 
elprop (2, 1) =10300000; 
elprop (2^  2)  =  .25; 
elprop (3, 1)=10300000; 
elprop (3  f  2 )  =  . 2  5 ; 
elprop (4, 1) =10300000; 
elprop (4^  2)  =  .25; 
elprop (5, 1) =10300000; 
elprop  (5, 2)  =  . 25; 
elprop {6 , 1)=10300000; 
elprop (6^  2)  =  .25; 
elprop  (7, 1) =10300000; 
elprop (7  f 2)  =  .25; 
elprop (8, 1)=10300000; 
elprop ( 8 , 2 ) = . 25 ; 
elprop (9, 1) =10300000; 
elprop  ( 9 , 2 )  =  . 2  5 ; 
elprop (10, 1) =10300000 ; 
elprop ( 10 , 2 ) = . 25 ; 
elprop (11, 1) =10300000 ; 
elprop (11,2)=. 25; 

nodes (1, 1) =1; nodes (1,2) =2 ; 
nodes (2,1) =2 ; nodes (2, 2) =7; 
nodes (3,1) =2; nodes (3, 2) =3; 
nodes ( 4 , 1 ) =3 ; nodes (4,2) =7 ; 
nodes (5,1) =3; nodes (5,2) =4 ; 
nodes (6,1) =4 ; nodes (6,2) =7 ; 
nodes (7,1) =4; nodes (7, 2) =5; 
nodes (8,1) =5 ; nodes (8, 2) =7; 
nodes (9,1) =5; nodes (9, 2) =6; 
nodes (10, X) =6; nodes (10, 2) =7; 
nodes (11, 1) =7; nodes (11, 2) =1; 


bcdof ( 1 ) =1 ; 
bcval ( 1 ) =0 ; 
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bcdof (2) =2; 
bcval (2 ) =0 ; 
bcdof (3) =11; 
bcval (3 ) =0 ; 
bcdof ( 4 ) =12  ; 
bcval ( 4 ) =0 ; 

ff=zeros (sdof , 1) / 
kk=zeros (sdof ,  sdof)  ; 
index=zeros (nnel*ndof ,  1)  ; 
elf orce=zeros (nnel*ndof ,  1)  ; 
eldisp=zeros (nnel*ndof , nnel*ndof ) ; 
stress=zeros (nel, 1) ; 

ff  (6)=-600; 
ff (8) =-600; 

for  iel=l : nel 

nd (1) =nodes  (iel, 1) ; 
nd (2 ) =nodes (iel, 2 ) ; 

xl=gcoord (nd ( 1 ) ,  1 ) ; yl=gcoord  (nd  ( 1 ) ,  2 ) 
x2=gcoord (nd ( 2 ) , 1 ) ; y2=gcoord (nd ( 2 ) r2) 

leng=sqrt ( (x2-xl) A2+ (y2-yl) A2) ; 

if  (x2-xl ) ==0 ; 

beta=2*atan (1) ; 

else 

beta=atan ( (y2-yl) / (x2-xl) ) ; 

end 

el=elprop  (iel, 1) ; 
area=elprop  (iel, 2 ) ; 

[ index] =FEELDOF (nd, nnel , ndof) ; 

k=FETRUSS2 (el, leng, area, 0 , beta, 1 ) ; 

kk=FEASMBLl (kk, k, index)  ; 


end 

[ kk, f f ] =FEAPLYC2 ( kk, f f , bcdof, bcval ) ; 

disp=kk\f f ; 

for  iel=l : nel 

nd (1) =nodes  (iel, 1) ; 
nd (2 ) =nodes  (iel, 2 ) ; 
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xl=gcoord (nd ( 1 ) ,  1 ) ; yl=gcoord (nd ( 1 )  ,  2 )  ; 
x2=gcoord (nd (2 ) , 1 ) ; y2=gcoord (nd (2 )  ,  2 )  ; 

leng=sqrt ( (x2-xl) A2+ (y2-yl) A2) / 

if  (x2-xl ) ==0 ; 

beta=2*atan (1) ; 

else 

beta=atan ( (y2-yl) / (x2-xl) ) ; 

end 

el=elprop  (iel, 1) ; 
area=elprop  (iel, 2) ; 
index=FEELDOF (nd, nnel, ndof ) ; 

k=FETRUSS2 (el, leng, area, 0 , beta, 1 ) ; 

for  i=l :  (nnel*ndof) 

eldisp (i) =disp (index (i) ) ; 

end 

elf or ce=k* eldisp ; 

stress (iel) =sqrt (elf orce ( 1 ) A2+elf orce (2 ) A2 ) /area 

if ( (x2-xl) *el force  (3) )  <0 ; 

stress (iel) ^-stress  (iel) ; 

end 


end 

num=l : 1 : sdof  ; 
displ=[num'  disp] 

numm=l : 1 : nel / 
stresses^ [numm 1  stress] 
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FEMD5 


nel=ll ; 
nnel=2 ; 
ndof=2 ; 
nnode=7 ; 

sdof =nnode*ndof ; 

gcoord ( 1 , 1 ) =0 ; gcoord ( 1 , 2 ) =0 ; 
gcoord (2 , 1) =0; gcoord (2 , 2) =  6; 
gcoord (3,1) =0; gcoord (3, 2) =  12; 
gcoord (4,1)=6.365; gcoord (4 , 2 ) =12 ; 
gcoord (5,1)=6.365; gcoord (5, 2) =  6; 
gcoord (6, i ) =  6 . 3  6  5 ; gcoord ( 6, 2 ) =  0 ; 
gcoord (7 , 1) =3 . 1825; gcoord (7 , 2) =  6; 

elprop (1, 1)=10300000; 
elprop ( 1 ,  2 ) = . 125 ; 
elprop (2, 1) =10300000; 
elprop (2 , 2 ) = . 125 ; 
elprop (3 , 1)=10300000; 
elprop (3  f  2 )  =  . 12  5 ; 
elprop (4, 1) =10300000; 
elprop ( 4 ,  2 ) = . 125 ; 
elprop (5, 1) =10300000; 
elprop  (5, 2)  =  . 12  5; 
elprop {6 , 1)=10300000; 
elprop (6^  2)  =  . 125; 
elprop  (7, 1) =10300000; 
elprop (7  f 2)  =  . 125; 
elprop (8, 1)=10300000; 
elprop ( 8 , 2 ) = . 125 ; 
elprop (9, 1) =10300000; 
elprop  (9,2)  =  .125; 
elprop (10a1)=10300000; 
elprop (10,2)=. 125; 
elprop (11, 1) =10300000; 
elprop (11,2)=. 125; 

nodes ( 1 , 1 ) =1 ; nodes (1,2) =2; 
nodes (2,1) =2; nodes (2 , 2) =7 ; 
nodes (3,1) =2; nodes (3, 2) =3; 
nodes ( 4 , 1 ) =3 ; nodes ( 4 , 2 ) =7 ; 
nodes (5,1) =3; nodes (5,2) =4 ; 
nodes (6,1) =4 ; nodes (6, 2) =7; 
nodes (7,1) =4; nodes (7, 2) =5; 
nodes (8,1) =5 ; nodes (8, 2) =7; 
nodes (9,1) =5; nodes (9, 2) =6; 
nodes (10, X) =6; nodes (10, 2) =7; 
nodes (11, 1) =7; nodes (11, 2) =1; 


bcdof ( 1 ) =1 ; 
bcval ( 1 ) =0 ; 
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bcdof (2) =2; 
bcval (2 ) =0 ; 
bcdof (3) =11; 
bcval (3 ) =0 ; 
bcdof ( 4 ) =12  ; 
bcval ( 4 ) =0 ; 

ff=zeros (sdof , 1) / 
kk=zeros (sdof ,  sdof)  ; 
index=zeros (nnel*ndof ,  1)  ; 
elf orce=zeros (nnel*ndof ,  1)  ; 
eldisp=zeros (nnel*ndof , nnel*ndof ) ; 
stress=zeros (nel, 1) ; 

ff  (6)=-600; 
ff (8) =-600; 

for  iel=l : nel 

nd (1) =nodes  (iel, 1) ; 
nd (2 ) =nodes (iel, 2 ) ; 

xl=gcoord (nd ( 1 ) ,  1 ) ; yl=gcoord  (nd  ( 1 ) ,  2 ) 
x2=gcoord (nd ( 2 ) , 1 ) ; y2=gcoord (nd ( 2 ) r2) 

leng=sqrt ( (x2-xl) A2+ (y2-yl) A2) ; 

if  (x2-xl ) ==0 ; 

beta=2*atan (1) ; 

else 

beta=atan ( (y2-yl) / (x2-xl) ) ; 

end 

el=elprop  (iel, 1) ; 
area=elprop  (iel, 2 ) ; 

[ index] =FEELDOF (nd, nnel , ndof) ; 

k=FETRUSS2 (el, leng, area, 0 , beta, 1 ) ; 

kk=FEASMBLl (kk, k, index)  ; 


end 

[ kk, f f ] =FEAPLYC2 ( kk, f f , bcdof, bcval ) ; 

disp=kk\f f ; 

for  iel=l : nel 

nd (1) =nodes  (iel, 1) ; 
nd (2 ) =nodes  (iel, 2 ) ; 
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xl=gcoord (nd ( 1 ) ,  1 ) ; yl=gcoord (nd ( 1 )  ,  2 )  ; 
x2=gcoord (nd (2 ) , 1 ) ; y2=gcoord (nd (2 )  ,  2 )  ; 

leng=sqrt ( (x2-xl) A2+ (y2-yl) A2) / 

if  (x2-xl ) ==0 ; 

beta=2*atan (1) ; 

else 

beta=atan ( (y2-yl) / (x2-xl) ) ; 

end 

el=elprop  (iel, 1) ; 
area=elprop  (iel, 2) ; 
index=FEELDOF (nd, nnel, ndof ) ; 

k=FETRUSS2 (el, leng, area, 0 , beta, 1 ) ; 

for  i=l :  (nnel*ndof) 

eldisp (i) =disp (index (i) ) ; 

end 

elf or ce=k* eldisp ; 

stress (iel) =sqrt (elf orce ( 1 ) A2+elf orce (2 ) A2 ) /area 

if ( (x2-xl) *el force  (3) )  <0 ; 

stress (iel) ^-stress  (iel) ; 

end 


end 

num=l : 1 : sdof  ; 
displ=[num'  disp] 

numm=l : 1 : nel / 
stresses^ [numm 1  stress] 
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FEMD25 


nel=ll ; 
nnel=2 ; 
ndof=2 ; 
nnode=7 ; 

sdof =nnode*ndof ; 

gcoord ( 1 , 1 ) =0 ; gcoord ( 1 , 2 ) =0 ; 
gcoord (2 , 1) =0; gcoord (2 , 2) =  6; 
gcoord (3,1) =0; gcoord (3, 2) =  12; 
gcoord (4,1)=6.365; gcoord (4 , 2 ) =12 ; 
gcoord (5,1)=6.365; gcoord (5, 2) =  6; 
gcoord (6, i ) =  6 . 3  6  5 ; gcoord ( 6, 2 ) =  0 ; 
gcoord (7 , 1) =3 . 1825; gcoord (7 , 2) =  6; 

elprop (1, 1)=10300000; 
elprop (1,2)=. 0625; 
elprop (2, 1) =10300000; 
elprop (2,2)=. 0625 ; 
elprop (3, 1)=10300000; 
elprop  (3, 2)  =  . 0625; 
elprop (4, 1) =10300000; 
elprop  (4, 2)  =  . 0625; 
elprop (5, 1) =10300000; 
elprop  (5, 2)  =  . 0625; 
elprop (6, 1)=10300000; 
elprop  (6, 2)  =  . 0625; 
elprop  (7, 1) =10300000; 
elprop (7,2)=. 0625 ; 
elprop (8, 1)=10300000; 
elprop (8, 2) =. 0625; 
elprop (9, 1) =10300000; 
elprop (9, 2) =. 0625; 
elprop (10, 1) =10300000 ; 
elprop (10, 2) =. 0625; 
elprop (11, 1) =10300000; 
elprop  (11, 2)  =  . 0625; 

nodes (1, 1) =1; nodes (1,2) =2 ; 
nodes (2,1) =2 ; nodes (2, 2) =7; 
nodes (3,1) =2; nodes (3, 2) =3; 
nodes ( 4 , 1 ) =3 ; nodes (4,2) =7 ; 
nodes (5,1) =3; nodes (5,2) =4 ; 
nodes (6,1) =4 ; nodes (6, 2) =7; 
nodes (7,1) =4; nodes (7, 2) =5; 
nodes (8,1) =5 ; nodes (8, 2) =7; 
nodes (9,1) =5; nodes (9, 2) =6; 
nodes (10, X) =6; nodes (10, 2) =7; 
nodes (11, 1) =7; nodes (11, 2) =1; 


bcdof ( 1 ) =1 ; 
bcval ( 1 ) =0 ; 
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bcdof (2) =2; 
bcval (2 ) =0 ; 
bcdof (3) =11; 
bcval (3 ) =0 ; 
bcdof ( 4 ) =12  ; 
bcval ( 4 ) =0 ; 

ff=zeros (sdof , 1) / 
kk=zeros (sdof ,  sdof)  ; 
index=zeros (nnel*ndof ,  1)  ; 
elf orce=zeros (nnel*ndof ,  1)  ; 
eldisp=zeros (nnel*ndof , nnel*ndof ) ; 
stress=zeros (nel, 1) ; 

ff  (6)=-600; 
ff (8) =-600; 

for  iel=l : nel 

nd (1) =nodes  (iel, 1) ; 
nd (2 ) =nodes (iel, 2 ) ; 

xl=gcoord (nd ( 1 ) ,  1 ) ; yl=gcoord  (nd  ( 1 ) ,  2 ) 
x2=gcoord (nd ( 2 ) , 1 ) ; y2=gcoord (nd ( 2 ) r2) 

leng=sqrt ( (x2-xl) A2+ (y2-yl) A2) ; 

if  (x2-xl ) ==0 ; 

beta=2*atan (1) ; 

else 

beta=atan ( (y2-yl) / (x2-xl) ) ; 

end 

el=elprop  (iel, 1) ; 
area=elprop  (iel, 2 ) ; 

[ index] =FEELDOF (nd, nnel , ndof) ; 

k=FETRUSS2 (el, leng, area, 0 , beta, 1 ) ; 

kk=FEASMBLl (kk, k, index)  ; 


end 

[ kk, f f ] =FEAPLYC2 ( kk, f f , bcdof, bcval ) ; 

disp=kk\f f ; 

for  iel=l : nel 

nd (1) =nodes  (iel, 1) ; 
nd (2 ) =nodes  (iel, 2 ) ; 
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xl=gcoord (nd ( 1 ) ,  1 ) ; yl=gcoord (nd ( 1 )  ,  2 )  ; 
x2=gcoord (nd (2 ) , 1 ) ; y2=gcoord (nd (2 )  ,  2 )  ; 

leng=sqrt ( (x2-xl) A2+ (y2-yl) A2) / 

if  (x2-xl ) ==0 ; 

beta=2*atan (1) ; 

else 

beta=atan ( (y2-yl) / (x2-xl) ) ; 

end 

el=elprop  (iel, 1) ; 
area=elprop  (iel, 2) ; 
index=FEELDOF (nd, nnel, ndof ) ; 

k=FETRUSS2 (el, leng, area, 0 , beta, 1 ) ; 

for  i=l :  (nnel*ndof) 

eldisp (i) =disp (index (i) ) ; 

end 

elf or ce=k* eldisp ; 

stress (iel) =sqrt (elf orce ( 1 ) A2+elf orce (2 ) A2 ) /area 

if ( (x2-xl) *el force  (3) )  <0 ; 

stress (iel) ^-stress  (iel) ; 

end 


end 

num=l : 1 : sdof  ; 
displ=[num'  disp] 

numm=l : 1 : nel / 
stresses^ [numm 1  stress] 
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FEME1 


nel=12 ; 
nnel=2 ; 
ndof=2 ; 
nnode=8 ; 

sdof =nnode*ndof ; 

gcoord ( 1 , 1 ) =0 ; gcoord (1^2) =0 ; 
gcoord (2,1) =  0 ; gcoord (2, 2) =  6; 
gcoord (3,1) =0; gcoord (3, 2) =  12; 
gcoord (4 , 1 ) =3 . 1825 ; gcoord (4 , 2 ) =12 ; 
gcoord (5,1)=6.365; gcoord ( 5 , 2 ) =12 ; 
gcoord (6, 1) =6 . 3 65; gcoord (6, 2) =6; 
gcoord (7, 1)=6.365; gcoord (7 , 2 ) =0 ; 
gcoord ( 8 , 1 ) =3 . 1825 ; gcoord ( 8 , 2 ) =  6 ; 

elprop (1, 1)=10300000; 
elprop (1, 2) =.25; 
elprop (2, 1) =10300000; 
elprop (2 , 2 ) = . 25; 
elprop  (3,  1)=10300000; 
elprop  (3, 2)  =  .  25,- 
elprop  (4,  1)  =10300000; 
elprop (4, 2) =.25; 
elprop (5, 1) =10300000; 
elprop  (5 , 2 )  =  . 2  5 ; 
elprop (6, 1) =10300000; 
elprop ( 6, 2 ) = . 25; 
elprop (7, 1) =10300000; 
elprop  (7 , 2 )  =  . 25 ; 
elprop (8, 1)=10300000; 
elprop (8, 2) =.25; 
elprop  (9,  1) =10300000; 
elprop ( 9, 2 ) = . 25 ; 
elprop (10, 1) =10300000; 
elprop ( 10 , 2 ) = . 25 ; 
elprop (11, 1) =10300000; 
elprop ( 11 , 2 ) = . 25 ; 
elprop (12, 1) =10300000 ; 
elprop (12,2)=. 25; 

nodes (1, 1) =1; nodes (1,2) =2 ; 
nodes (2,1) =2 ; nodes (2, 2) =8; 
nodes (3,1) =2; nodes (3,2) —4 ; 
nodes (4,1) =2 ; nodes (4,2) =3 ; 
nodes (5, 1) =3; nodes (5, 2) =4; 
nodes (6,1) =4; nodes (6, 2) =6; 
nodes (7, 1) =4; nodes (7, 2) =5; 
nodes (8,1) =5 ; nodes (8, 2) =6; 
nodes (9,1) =6; nodes (9, 2) =8; 
nodes (10, 1) =6; nodes (10, 2) =7; 
nodes ( 11 , 1 ) =7 ; nodes (11, 2) =8; 
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nodes ( 12 , 1 ) =8 ; nodes ( 12 , 2 ) =1 ; 


bcdof ( 1 ) =1 ; 
bcval ( 1 ) =0 ; 
bcdof (2 ) =2 ; 
bcval (2 ) =0 ; 
bcdof (3)=13; 
bcval ( 3 ) =0 ; 
bcdof (4) =14; 
bcval ( 4 ) =0 ; 

ff=zeros (sdof , 1) / 
kk=zeros (sdof, sdof) ; 
index=zeros (nnel*ndof ,  1)  ; 
elf orce=zeros (nnel*ndof ,  1)  ; 
eldisp=zeros (nnel*ndof , nnel*ndof ) ; 
stress=zeros (nel, 1) ; 

ff  (6)  ^400,- 
ff  (8)  =-400; 
ff  (10)— 400; 

for  iel=l : nel 

nd ( 1 ) =nodes ( iel , 1 ) ; 
nd (2 ) =nodes (iel, 2) ; 

xl=gcoord (nd ( 1 ) ,  1 ) ; yl=gcoord (nd  ( 1 ) ,2) 
x2=gcoord (nd (2 ) , 1 ) ; y2=gcoord (nd (2 ) r2) 

leng=sqrt ( (x2-xl) A2+ (y2-yl) A2) ; 

if  (x2-xl ) ==0 ; 

beta=2*atan (1) ; 

else 

beta=atan ( (y2-yl) / (x2-xl) ) ; 

end 

el=elprop  (iel, 1) ; 
area=elprop  (iel, 2 ) ; 

[ index] =FEELDOF (nd, nnel , ndof) ; 

k=FETRUSS2 (el , leng, area, 0 , beta, 1 ) ; 

kk=FEASMBLl (kk, k, index)  ; 


end 

[kk, f f ] =FEAPLYC2 (kk, f f , bcdof, bcval) ; 
disp=kk\f f ; 
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for  iel=l:nel 


nd (1) =nodes (iel, 1)  ; 
nd (2 ) =nodes (iel, 2 )  ; 

xl=gcoord (nd ( 1 ) , 1 ) ; yl=gcoord (nd ( 1 )  ,  2 )  ; 
x2=gcoord (nd (2 ) , 1 ) ; y2=gcoord (nd (2 )  ,  2 )  ; 

leng=sqrt ( (x2-xl) A2+ (y2-yl) A2) / 

if  (x2-xl ) ==0 ; 

beta=2*atan (1) ; 

else 

beta=atan ( (y2-yl) / (x2-xl) ) ; 

end 

el=elprop  (iel, 1) ; 
area=elprop  (iel, 2) ; 
index=FEELDOF (nd, nnel , ndof ) ; 

k=FETRUSS2 (el, leng, area, 0 , beta, 1 ) ; 

for  i=l :  (nnel*ndof) 

eldisp (i) =disp (index (i)  )  ; 

end 

elf or ce=k* eldisp ; 

stress (iel) =sqrt (elforce (1) A2+elforce (2) A2) /area 

if ( (x2-xl) *el force  (3) )  <0 ; 

stress (iel) ^-stress  (iel) ; 

end 


end 

num=l : 1 : sdof ; 
displ=[num'  disp] 

numm=l : 1 : nel / 
stresses^ [numm '  stress] 
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FEME  5 


nel=12 ; 
nnel=2 ; 
ndof=2 ; 
nnode=8 ; 

sdof =nnode*ndof ; 

gcoord ( 1 , 1 ) =0 ; gcoord (1^2) =0 ; 
gcoord (2,1) =  0 ; gcoord (2, 2) =  6; 
gcoord (3,1) =0; gcoord (3, 2) =  12; 
gcoord (4 , 1 ) =3 . 1825 ; gcoord (4 , 2 ) =12 ; 
gcoord (5,1)=6.365; gcoord ( 5 , 2 ) =12 ; 
gcoord (6, 1) =6 . 3 65; gcoord (6, 2) =6; 
gcoord (7, 1)=6.365; gcoord (7 , 2 ) =0 ; 
gcoord ( 8 , 1 ) =3 . 1825 ; gcoord ( 8 , 2 ) =  6 ; 

elprop (1, 1)=10300000; 
elprop (1,2)=. 125; 
elprop (2, 1) =10300000; 
elprop  (2 , 2 )  =  . 12  5 ; 
elprop  (3,  1)=10300000; 
elprop (3,2)=. 125; 
elprop (4, 1) =10300000; 
elprop (4,2)=. 125; 
elprop (5, 1) =10300000; 
elprop (5 , 2 ) = . 125 ; 
elprop (6, 1) =10300000; 
elprop ( 6, 2 ) = . 125 ; 
elprop (7, 1) =10300000; 
elprop  (7 , 2 )  =  . 12  5 ; 
elprop  (8, 1)=10300000; 
elprop ( 8 , 2 ) = . 125 ; 
elprop (9, 1) =10300000; 
elprop  ( 9, 2 )  =  . 12  5 ; 
elprop (10, 1) =10300000 ; 
elprop (10,2)=. 125; 
elprop (11,1)=10300000; 
elprop (11, 2) =. 125; 
elprop (12, 1) =10300000 ; 
elprop (12,2)=. 125; 

nodes (1, 1) =1; nodes (1,2) =2 ; 
nodes (2,1) =2 ; nodes (2, 2) =8; 
nodes (3,1) =2; nodes (3,2) —4 ; 
nodes (4,1) =2 ; nodes (4,2) =3 ; 
nodes (5, 1) =3; nodes (5, 2) =4; 
nodes (6,1) =4; nodes (6, 2) =6; 
nodes (7, 1) =4; nodes (7, 2) =5; 
nodes (8,1) =5 ; nodes (8, 2) =6; 
nodes (9,1) =6; nodes (9, 2) =8; 
nodes (10, 1) =6; nodes (10, 2) =7; 
nodes ( 11 , 1 ) =7 ; nodes (11, 2) =8; 
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nodes ( 12 , 1 ) =8 ; nodes ( 12 , 2 ) =1 ; 


bcdof ( 1 ) =1 ; 
bcval ( 1 ) =0 ; 
bcdof (2 ) =2 ; 
bcval (2 ) =0 ; 
bcdof (3)=13; 
bcval ( 3 ) =0 ; 
bcdof (4) =14; 
bcval ( 4 ) =0 ; 

ff=zeros (sdof , 1) / 
kk=zeros (sdof, sdof) ; 
index=zeros (nnel*ndof ,  1)  ; 
elf orce=zeros (nnel*ndof ,  1)  ; 
eldisp=zeros (nnel*ndof , nnel*ndof ) ; 
stress=zeros (nel, 1) ; 

ff  (6)  ^400,- 
ff  (8)  =-400; 
ff  (10)— 400; 

for  iel=l : nel 

nd ( 1 ) =nodes ( iel , 1 ) ; 
nd (2 ) =nodes (iel, 2) ; 

xl=gcoord (nd ( 1 ) ,  1 ) ; yl=gcoord (nd  ( 1 ) ,2) 
x2=gcoord (nd (2 ) , 1 ) ; y2=gcoord (nd (2 ) r2) 

leng=sqrt ( (x2-xl) A2+ (y2-yl) A2) ; 

if  (x2-xl ) ==0 ; 

beta=2*atan (1) ; 

else 

beta=atan ( (y2-yl) / (x2-xl) ) ; 

end 

el=elprop  (iel, 1) ; 
area=elprop  (iel, 2 ) ; 

[ index] =FEELDOF (nd, nnel , ndof) ; 

k=FETRUSS2 (el , leng, area, 0 , beta, 1 ) ; 

kk=FEASMBLl (kk, k, index)  ; 


end 

[kk, f f ] =FEAPLYC2 (kk, f f , bcdof, bcval) ; 
disp=kk\f f ; 
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for  iel=l:nel 


nd (1) =nodes (iel, 1)  ; 
nd (2 ) =nodes (iel, 2 )  ; 

xl=gcoord (nd ( 1 ) , 1 ) ; yl=gcoord (nd ( 1 )  ,  2 )  ; 
x2=gcoord (nd (2 ) , 1 ) ; y2=gcoord (nd (2 )  ,  2 )  ; 

leng=sqrt ( (x2-xl) A2+ (y2-yl) A2) / 

if  (x2-xl ) ==0 ; 

beta=2*atan (1) ; 

else 

beta=atan ( (y2-yl) / (x2-xl) ) ; 

end 

el=elprop  (iel, 1) ; 
area=elprop  (iel, 2) ; 
index=FEELDOF (nd, nnel , ndof ) ; 

k=FETRUSS2 (el, leng, area, 0 , beta, 1 ) ; 

for  i=l :  (nnel*ndof) 

eldisp (i) =disp (index (i)  )  ; 

end 

elf or ce=k* eldisp ; 

stress (iel) =sqrt (elforce (1) A2+elforce (2) A2) /area 

if ( (x2-xl) *el force  (3) )  <0 ; 

stress (iel) ^-stress  (iel) ; 

end 


end 

num=l : 1 : sdof ; 
displ=[num'  disp] 

numm=l : 1 : nel / 
stresses^ [numm '  stress] 
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FEME25 


nel=12 ; 
nnel=2 ; 
ndof=2 ; 
nnode=8 ; 

sdof =nnode*ndof ; 

gcoord ( 1 , 1 ) =0 ; gcoord (1^2) =0 ; 
gcoord (2,1) =  0 ; gcoord (2, 2) =  6; 
gcoord (3,1) =0; gcoord (3, 2) =  12; 
gcoord (4 , 1 ) =3 . 1825 ; gcoord (4 , 2 ) =12 ; 
gcoord (5,1)=6.365; gcoord ( 5 , 2 ) =12 ; 
gcoord (6, 1) =6 . 3 65; gcoord (6, 2) =6; 
gcoord (7, 1)=6.365; gcoord (7 , 2 ) =0 ; 
gcoord ( 8 , 1 ) =3 . 1825 ; gcoord ( 8 , 2 ) =  6 ; 

elprop (1, 1)=10300000; 
elprop (1,2)=. 0625; 
elprop (2, 1) =10300000; 
elprop  (2,2)=. 0625; 
elprop  (3,  1)=10300000; 
elprop (3,2)=.0625; 
elprop (4, 1) =10300000; 
elprop (4,2)=. 0625 ; 
elprop (5, 1) =10300000; 
elprop  (5, 2)  =  . 0625; 
elprop (6, 1) =10300000; 
elprop (6, 2) =. 0625; 
elprop (7, 1) =10300000; 
elprop  (7,2)=. 0625; 
elprop  (8, 1)=10300000; 
elprop  (8, 2)  =  . 0625; 
elprop (9, 1) =10300000; 
elprop  (9,2)  =  .0625; 
elprop (10, 1) =10300000 ; 
elprop (10, 2) =. 0625; 
elprop (11,1)=10300000; 
elprop (11, 2) =. 0625; 
elprop (12, 1) =10300000 ; 
elprop (12, 2) =. 0625; 

nodes (1, 1) =1; nodes (1,2) =2 ; 
nodes (2,1) =2 ; nodes (2, 2) =8; 
nodes (3,1) =2; nodes (3,2) —4 ; 
nodes (4,1) =2 ; nodes (4,2) =3 ; 
nodes (5, 1) =3; nodes (5, 2) =4; 
nodes (6,1) =4; nodes (6, 2) =6; 
nodes (7, 1) =4; nodes (7, 2) =5; 
nodes (8,1) =5 ; nodes (8, 2) =6; 
nodes (9,1) =6; nodes (9, 2) =8; 
nodes (10, 1) =6; nodes (10, 2) =7; 
nodes ( 11 , 1 ) =7 ; nodes (11, 2) =8; 
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nodes ( 12 , 1 ) =8 ; nodes ( 12 , 2 ) =1 ; 


bcdof ( 1 ) =1 ; 
bcval ( 1 ) =0 ; 
bcdof (2 ) =2 ; 
bcval (2 ) =0 ; 
bcdof (3)=13; 
bcval ( 3 ) =0 ; 
bcdof (4) =14; 
bcval ( 4 ) =0 ; 

ff=zeros (sdof , 1) / 
kk=zeros (sdof, sdof) ; 
index=zeros (nnel*ndof ,  1)  ; 
elf orce=zeros (nnel*ndof ,  1)  ; 
eldisp=zeros (nnel*ndof , nnel*ndof ) ; 
stress=zeros (nel, 1) ; 

ff  (6)  ^400,- 
ff  (8)  =-400; 
ff  (10)— 400; 

for  iel=l : nel 

nd ( 1 ) =nodes ( iel , 1 ) ; 
nd (2 ) =nodes (iel, 2) ; 

xl=gcoord (nd ( 1 ) ,  1 ) ; yl=gcoord (nd  ( 1 ) ,2) 
x2=gcoord (nd (2 ) , 1 ) ; y2=gcoord (nd (2 ) r2) 

leng=sqrt ( (x2-xl) A2+ (y2-yl) A2) ; 

if  (x2-xl ) ==0 ; 

beta=2*atan (1) ; 

else 

beta=atan ( (y2-yl) / (x2-xl) ) ; 

end 

el=elprop  (iel, 1) ; 
area=elprop  (iel, 2 ) ; 

[ index] =FEELDOF (nd, nnel , ndof) ; 

k=FETRUSS2 (el , leng, area, 0 , beta, 1 ) ; 

kk=FEASMBLl (kk, k, index)  ; 


end 

[kk, f f ] =FEAPLYC2 (kk, f f , bcdof, bcval) ; 
disp=kk\f f ; 


94 


for  iel=l:nel 


nd (1) =nodes (iel, 1)  ; 
nd (2 ) =nodes (iel, 2 )  ; 

xl=gcoord (nd ( 1 ) , 1 ) ; yl=gcoord (nd ( 1 )  ,  2 )  ; 
x2=gcoord (nd (2 ) , 1 ) ; y2=gcoord (nd (2 )  ,  2 )  ; 

leng=sqrt ( (x2-xl) A2+ (y2-yl) A2) / 

if  (x2-xl ) ==0 ; 

beta=2*atan (1) ; 

else 

beta=atan ( (y2-yl) / (x2-xl) ) ; 

end 

el=elprop  (iel, 1) ; 
area=elprop  (iel, 2) ; 
index=FEELDOF (nd, nnel , ndof ) ; 

k=FETRUSS2 (el, leng, area, 0 , beta, 1 ) ; 

for  i=l :  (nnel*ndof) 

eldisp (i) =disp (index (i)  )  ; 

end 

elf or ce=k* eldisp ; 

stress (iel) =sqrt (elforce (1) A2+elforce (2) A2) /area 

if ( (x2-xl) *el force  (3) )  <0 ; 

stress (iel) ^-stress  (iel) ; 

end 


end 

num=l : 1 : sdof ; 
displ=[num'  disp] 

numm=l : 1 : nel / 
stresses^ [numm '  stress] 
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FEMF1 


nel=7 ; 
nnel=2 ; 
ndof=2 ; 
nnode=5 ; 

sdof =nnode*ndof ; 

gcoord ( 1 , 1 ) =0 ; gcoord ( 1 , 2 ) =0 ; 
gcoord (2 , 1) =0 ; gcoord (2 , 2 ) =12  ; 
gcoord (3, 1) =6 . 3 65; gcoord (3, 2) =  12; 
gcoord (4,1)=6.365; gcoord (4 , 2 ) =0 ; 
gcoord  (5 , 1 ) =3 . 1825 ;  gcoord  (5 , 2 ) =6; 

elprop (1, 1) =10300000; 
elprop ( 1 , 2 ) = . 25 ; 
elprop (2 ,  1)=10300000; 
elprop  (2  ,  2 )  =  . 2  5 ; 
elprop (3,1) =10300000; 
elprop  (3, 2)  =  . 25; 
elprop (4, 1)=10300000; 
elprop  ( 4  ,  2 )  =  . 2  5 ; 
elprop (5, 1) =10300000; 
elprop ( 5 ,  2)  =  .25; 
elprop (6, 1) =10300000; 
elprop ( 6 ,  2)  =  .25; 
elprop (7 , 1)=10300000; 
elprop  ( 7 , 2 )  =  . 2  5 ; 

nodes ( 1 , 1 ) =1 ; nodes (1,2) =2; 
nodes (2,1) =2 ; nodes (2 , 2) =5; 
nodes (3,1) =2; nodes (3, 2) =3; 
nodes (4,1) =3 ; nodes (4,2) =5 ; 
nodes (5,1) =3; nodes (5, 2) =4; 
nodes (6, 1) =4; nodes (6, 2) =5; 
nodes (7,1) =5; nodes (7, 2) =1; 

bcdof ( 1 ) =1 ; 
bcval ( 1 ) =0 ; 
bcdof (2) =2; 
bcval (2 ) =0 ; 
bcdof (3) =7; 
bcval (3 ) =0 ; 
bcdof (4) =8; 
bcval ( 4 ) =0 ; 

ff=zeros (sdof , 1) ; 
kk=zeros (sdof,  sdof)  ; 
index=zeros (nnel*ndof ,  1)  ; 
elf orce=zeros (nnel*ndof ,  1)  ; 
eldisp=zeros (nnel*ndof , nnel*ndof ) 
stress=zeros (nel, 1) ; 


ff (4)=-600; 
ff (6)=-600; 


for  iel=l : nel 

nd (1) =nodes  (iel, 1) ; 
nd (2 ) =nodes (iel, 2 ) ; 

xl=gcoord (nd ( 1 ) ,  1 ) ; yl=gcoord (nd ( 1 ) ,  2 ) 
x2=gcoord (nd (2 ) ,  1 ) ; y2=gcoord (nd (2 ) ,2) 

leng=sqrt ( (x2-xl) A2+ (y2-yl) A2) ; 

if  (x2-xl ) ==0 ; 

beta=2*atan (1) ; 

else 

beta=atan ( (y2-yl) / (x2-xl) ) ; 

end 

el=elprop  (iel, 1) ; 
area=elprop  (ielf  2) ; 

[ index] =FEELDOF (ndf nnel A ndof) ; 

k=FETRUSS2 (el,leng, area, 0 , beta, 1 ) ; 

kk=FEASMBLl (kk, k, index)  ; 


end 

[ kk, f f ] =FEAPLYC2 ( kk, f f , bcdof , bcval ) ; 

disp=kk\f f ; 

for  iel=l : nel 

nd (1) =nodes  (iel, 1) ; 
nd (2 ) =nodes  (iel, 2 ) ; 

xl=gcoord (nd ( 1 ) , 1 ) ; yl=gcoord (nd  ( 1 ) , 2 ) 
x2=gcoord (nd (2 ) , 1 ) ; y2=gcoord (nd (2 )  ,  2 ) 

leng=sqrt ( (x2-xl) A2+ (y2-yl) A2) ; 

if  (x2-xl ) ==0 ; 

beta=2^atan (1) ; 

else 

beta=atan ( (y2-yl) / (x2-xl) ) / 

end 

el=elprop  (iel, 1) ; 
area=elprop  (iel, 2 ) ; 
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index=FEELDOF (nd, nnel , ndof ) ; 


k=FETRUSS2 (el, leng, area, 0 , beta,  1 ) ; 

for  i=l :  (nnel*ndof) 

eldisp (i) =disp (index (i)  )  ; 

end 

elf orce=k*eldisp; 

stress (iel) =sqrt (elforce (1) A2+elforce (2) A2) /area 

if ( (x2-xl ) *el force  ( 3 ) )  <0 ; 

stress (iel) ^-stress  (iel) ; 

end 


end 

num=l : 1 : sdof  ; 
displ=[num'  disp] 

numm=l : 1 : nel / 
stresses^ [numm 1  stress] 
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FEMF5 


nel=7 ; 
nnel=2 ; 
ndof=2 ; 
nnode=5 ; 

sdof =nnode*ndof ; 

gcoord ( 1 , 1 ) =0 ; gcoord ( 1 , 2 ) =0 ; 
gcoord (2 , 1) =0 ; gcoord (2 , 2 ) =12  ; 
gcoord (3, 1) =6 . 3 65; gcoord (3, 2) =  12; 
gcoord (4,1)=6.365; gcoord (4 , 2 ) =0 ; 
gcoord  (5 , 1 ) =3 . 1825 ;  gcoord  (5 , 2 ) =6; 

elprop (1, 1) =10300000; 
elprop ( 1 , 2 ) = . 125 ; 
elprop (2 ,  1)=10300000; 
elprop  (2  ,  2 )  =  . 12  5 ; 
elprop (3,1) =10300000; 
elprop  (3, 2)  =  . 12  5; 
elprop (4, 1)=10300000; 
elprop  ( 4  ,  2 )  =  . 12  5 ; 
elprop (5, 1) =10300000; 
elprop ( 5 ,  2)  =  . 125; 
elprop (6, 1) =10300000; 
elprop ( 6 , 2 ) = . 125 ; 
elprop (7 , 1)=10300000; 
elprop  (7 , 2)  =  . 125; 

nodes ( 1 , 1 ) =1 ; nodes ( 1 , 2 ) =2 ; 
nodes (2,1) =2 ; nodes (2 ,  2 ) =5 ; 
nodes (3,1) =2; nodes (3, 2) =3; 
nodes (4,1) =3 ; nodes (4,2) =5 ; 
nodes (5,1) =3; nodes (5, 2) =4; 
nodes (6, 1) =4; nodes (6, 2) =5; 
nodes (7,1) =5; nodes (7, 2) =1; 

bcdof ( 1 ) =1 ; 
bcval ( 1 ) =0 ; 
bcdof (2) =2; 
bcval (2 ) =0 ; 
bcdof (3) =7; 
bcval (3 ) =0 ; 
bcdof (4) =8; 
bcval ( 4 ) =0 ; 

ff=zeros (sdof , 1) ; 
kk=zeros (sdof,  sdof)  ; 
index=zeros (nnel*ndof ,  1)  ; 
elf orce=zeros (nnel*ndof ,  1)  ; 
eldisp=zeros (nnel*ndof , nnel*ndof ) 
stress=zeros (nel, 1) ; 


ff (4)=-600; 
ff (6)=-600; 


for  iel=l : nel 

nd (1) =nodes  (iel, 1) ; 
nd (2 ) =nodes (iel, 2 ) ; 

xl=gcoord (nd ( 1 ) ,  1 ) ; yl=gcoord (nd ( 1 ) ,  2 ) 
x2=gcoord (nd (2 ) ,  1 ) ; y2=gcoord (nd (2 ) ,2) 

leng=sqrt ( (x2-xl) A2+ (y2-yl) A2) ; 

if  (x2-xl ) ==0 ; 

beta=2*atan (1) ; 

else 

beta=atan ( (y2-yl) / (x2-xl) ) ; 

end 

el=elprop  (iel, 1) ; 
area=elprop  (ielf  2) ; 

[ index] =FEELDOF (ndf nnel A ndof) ; 

k=FETRUSS2 (el,leng, area, 0 , beta, 1 ) ; 

kk=FEASMBLl (kk, k, index)  ; 


end 

[ kk, f f ] =FEAPLYC2 ( kk, f f , bcdof , bcval ) ; 

disp=kk\f f ; 

for  iel=l : nel 

nd (1) =nodes  (iel, 1) ; 
nd (2 ) =nodes  (iel, 2 ) ; 

xl=gcoord (nd ( 1 ) , 1 ) ; yl=gcoord (nd  ( 1 ) , 2 ) 
x2=gcoord (nd (2 ) , 1 ) ; y2=gcoord (nd (2 )  ,  2 ) 

leng=sqrt ( (x2-xl) A2+ (y2-yl) A2) ; 

if  (x2-xl ) ==0 ; 

beta=2^atan (1) ; 

else 

beta=atan ( (y2-yl) / (x2-xl) ) / 

end 

el=elprop  (iel, 1) ; 
area=elprop  (iel, 2 ) ; 
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index=FEELDOF (nd, nnel , ndof ) ; 


k=FETRUSS2 (el, leng, area, 0 , beta,  1 ) ; 

for  i=l :  (nnel*ndof) 

eldisp (i) =disp (index (i)  )  ; 

end 

elf orce=k*eldisp; 

stress (iel) =sqrt (elforce (1) A2+elforce (2) A2) /area 

if ( (x2-xl ) *el force  ( 3 ) )  <0 ; 

stress (iel) ^-stress  (iel) ; 

end 


end 

num=l : 1 : sdof  ; 
displ=[num'  disp] 

numm=l : 1 : nel ; 
stresses^ [numm 1  stress] 
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FEMF25 


nel=7 ; 
nnel=2 ; 
ndof=2 ; 
nnode=5 ; 

sdof =nnode*ndof ; 

gcoord ( 1 , 1 ) =0 ; gcoord ( 1 , 2 ) =0 ; 
gcoord (2 , 1) =0 ; gcoord (2 , 2 ) =12  ; 
gcoord (3, 1) =6 . 3 65; gcoord (3, 2) =  12; 
gcoord (4,1)=6.365; gcoord (4 , 2 ) =0 ; 
gcoord  (5 , 1 ) =3 . 1825 ;  gcoord  (5 , 2 ) =6; 

elprop (1, 1) =10300000; 
elprop (1,2)=. 0625 ; 
elprop (2, 1)=10300000; 
elprop  (2,2)  =  .0625; 
elprop (3, 1) =10300000 ; 
elprop  (3,2)  =  . 0625 ; 
elprop (4, 1)=10300000; 
elprop  (4,2)  =  .0625; 
elprop (5, 1) =10300000; 
elprop  (5, 2)  =  . 0625; 
elprop (6, 1) =10300000; 
elprop (6,2)=.0625; 
elprop (7, 1)=10300000; 
elprop  (7, 2)  =  . 0625; 

nodes (1, 1) =1; nodes (1,2) =2 ; 
nodes (2,1) =2 ; nodes (2,2) =5 ; 
nodes (3,1) =2; nodes (3, 2) =3; 
nodes (4,1) =3 ; nodes (4, 2) =5; 
nodes (5,1) =3; nodes (5, 2) =4; 
nodes (6, 1) =4; nodes (6, 2) =5; 
nodes (7,1) =5; nodes (7, 2) =1; 

bcdof ( 1 ) =1 ; 
bcval ( 1 ) =0 ; 
bcdof (2) =2; 
bcval (2 ) =0 ; 
bcdof (3) =7; 
bcval  (3 ) =0 ; 
bcdof (4) =8; 
bcval ( 4 ) =0 ; 

ff=zeros (sdof , 1) ; 
kk=zeros (sdof,  sdof)  ; 
index=zeros (nnel*ndof ,  1)  ; 
elf orce=zeros (nnel*ndof ,  1)  ; 
eldisp=zeros (nnel*ndof , nnel*ndof ) 
stress=zeros (nel, 1) ; 


ff (4)=-600; 
ff (6)=-600; 


for  iel=l : nel 

nd (1) =nodes  (iel, 1) ; 
nd (2 ) =nodes (iel, 2 ) ; 

xl=gcoord (nd ( 1 ) ,  1 ) ; yl=gcoord (nd ( 1 ) ,  2 ) 
x2=gcoord (nd (2 ) ,  1 ) ; y2=gcoord (nd (2 ) ,2) 

leng=sqrt ( (x2-xl) A2+ (y2-yl) A2) ; 

if  (x2-xl ) ==0 ; 

beta=2*atan (1) ; 

else 

beta=atan ( (y2-yl) / (x2-xl) ) ; 

end 

el=elprop  (iel, 1) ; 
area=elprop  (ielf  2) ; 

[ index] =FEELDOF (ndf nnel A ndof) ; 

k=FETRUSS2 (el,leng, area, 0 , beta, 1 ) ; 

kk=FEASMBLl (kk, k, index)  ; 


end 

[ kk, f f ] =FEAPLYC2 ( kk, f f , bcdof , bcval ) ; 

disp=kk\f f ; 

for  iel=l : nel 

nd (1) =nodes  (iel, 1) ; 
nd (2 ) =nodes  (iel, 2 ) ; 

xl=gcoord (nd ( 1 ) , 1 ) ; yl=gcoord (nd  ( 1 ) , 2 ) 
x2=gcoord (nd (2 ) , 1 ) ; y2=gcoord (nd (2 )  ,  2 ) 

leng=sqrt ( (x2-xl) A2+ (y2-yl) A2) ; 

if  (x2-xl ) ==0 ; 

beta=2^atan (1) ; 

else 

beta=atan ( (y2-yl) / (x2-xl) ) / 

end 

el=elprop  (iel, 1) ; 
area=elprop  (iel, 2 ) ; 
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index=FEELDOF (nd, nnel , ndof ) ; 


k=FETRUSS2 (el, leng, area, 0 , beta,  1 ) ; 

for  i=l :  (nnel*ndof) 

eldisp (i) =disp (index (i)  )  ; 

end 

elf orce=k*eldisp; 

stress (iel) =sqrt (elforce (1) A2+elforce (2) A2) /area 

if ( (x2-xl ) *el force  ( 3 ) )  <0 ; 

stress (iel) ^-stress  (iel) ; 

end 


end 

num=l : 1 : sdof  ; 
displ=[num'  disp] 

numm=l : 1 : nel / 
stresses^ [numm 1  stress] 
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FEMG1 


nel=13 ; 
nnel=2 ; 
ndof=2 ; 
nnode=8 ; 

sdof =nnode*ndof ; 

gcoord ( 1 , 1 ) =0 ; gcoord (1^2) =0 ; 
gcoord (2,1) =  0 ; gcoord (2, 2) =  6; 
gcoord (3,1) =0; gcoord (3, 2) =  12; 
gcoord (4,1)=6.365; gcoord (4 , 2 ) =12 ; 
gcoord (5,1)=6.365; gcoord (5, 2) =6; 
gcoord (6, I) =6 . 365; gcoord ( 6, 2 ) =0 ; 
gcoord (7 , 1 ) =3 . 1825 ; gcoord (7 , 2 ) =3; 
gcoord ( 8 , 1 ) =3 . 1825 ; gcoord ( 8 , 2 ) =9 ; 

elprop (1, 1)=10300000; 
elprop (1, 2) =.25; 
elprop (2, 1) =10300000; 
elprop (2 , 2 ) = . 25; 
elprop  (3,  1)=10300000; 
elprop  (3, 2)  =  .  25,- 
elprop  (4,  1)  =10300000; 
elprop (4, 2) =.25; 
elprop (5, 1) =10300000; 
elprop  (5 , 2 )  =  . 2  5 ; 
elprop (6, 1) =10300000; 
elprop ( 6, 2 ) = . 25; 
elprop (7, 1) =10300000; 
elprop  (7 , 2 )  =  . 25 ; 
elprop (8, 1)=10300000; 
elprop (8, 2) =.25; 
elprop  (9,  1) =10300000; 
elprop ( 9, 2 ) = . 25 ; 
elprop (10, 1) =10300000; 
elprop ( 10 , 2 ) = . 25 ; 
elprop (11, 1) =10300000; 
elprop ( 11 , 2 ) = . 25 ; 
elprop (12, 1) =10300000 ; 
elprop ( 12 , 2 ) = . 25 ; 
elprop (13, 1) =10300000 ; 
elprop (13,2)=. 25; 

nodes  (1,1)  =71 ;  nodes  (1,2)  =2 ; 
nodes (2,1) =2 ; nodes (2,2) =7 ; 
nodes (3,1) =2; nodes (3, 2) =8; 
nodes (4,1) =2 ; nodes (4,2) =3 ; 
nodes (5,1) =3; nodes (5, 2) =8; 
nodes (6, 1) =3; nodes (6,2) =4 ; 
nodes (7,1) =4; nodes (7, 2) =8; 
nodes (8, 1) =4; nodes (8,2) =5 ; 
nodes (9, 1) =5; nodes (9, 2) =8; 
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nodes ( 10 , 1 ) =5 ; nodes ( 10 , 2 ) =7 ; 
nodes ( 11 , 1 ) =5 ; nodes ( 11 , 2 ) =  6; 
nodes ( 12 , 1 ) =6; nodes ( 12 , 2 ) =7 ; 
nodes ( 13 , 1) =7 ; nodes (13, 2) =1/ 

bcdof (1) =1; 
bcval ( 1 ) =0 ; 
bcdof (2) =2; 
bcval (2 ) =0 ; 
bcdof (3) =11; 
bcval (3 ) =0 ; 
bcdof (4)=12; 
bcval ( 4 ) =0 ; 

ff=zeros (sdof, 1) ; 
kk=zeros (sdof ,  sdof ) ; 
index=zeros (nnel*ndof ,  1)  ; 
elf orce=zeros (nnel*ndof ,  1)  ; 
eldisp=zeros (nnel*ndof , nnel*ndof ) ; 
stress=zeros (nel, 1) ; 

ff  (6)=-600; 
ff (8) =-600; 

for  iel=l : nel 

nd (1) =nodes  (iel, 1) ; 
nd (2 ) =nodes (iel, 2 ) ; 

xl=gcoord (nd ( 1 ) ,  1 ) ; yl=gcoord  (nd  ( 1 ) ,  2 ) 
x2=gcoord (nd ( 2 ) , 1 ) ; y2=gcoord (nd ( 2 ) ,2) 

leng=sqrt ( (x2-xl) A2+ (y2-yl) A2) ; 

if  (x2-xl ) ==0 ; 

beta=2*atan (1) ; 

else 

beta=atan ( (y2-yl) / (x2-xl) ) ; 

end 

el=elprop  (iel, 1) ; 
area=elprop  ( iel , 2 ) ; 

[ index] =FEELDOF (nd, nnel , ndof) ; 

k=FETRUSS2 (el , leng, area, 0 , beta, 1 ) ; 

kk=FEASMBLl (kk, k, index)  ; 


end 

[ kk, f f ] =FEAPLYC2 ( kk, f f , bcdof, bcval ) ; 


disp=kk\f f ; 
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for  iel=l : nel 


nd (1) =nodes (iel, 1)  ; 
nd (2 ) =nodes (iel,  2 )  ; 

xl=gcoord (nd ( 1 ) ,  1 ) ; yl=gcoord (nd ( 1 )  ,  2 )  ; 
x2=gcoord (nd (2 ) ,  1 ) ; y2=gcoord (nd (2 )  ,2)  ; 

leng=sqrt ( (x2-xl) A2+ (y2-yl) A2) ; 

if  (x2-xl ) ==0 ; 

beta=2*atan (1) ; 

else 

beta=atan ( (y2-yl) / (x2-xl) ) ; 

end 

el=elprop  (iel, 1) ; 
area=elprop  (iel, 2 ) ; 
index=FEELDOF (nd, nnel , ndof ) ; 

k=FETRUSS2 (el, leng, area, 0 , beta, 1 ) ; 

for  i=l :  (nnel*ndof) 

eldisp (i) =disp (index (i) ) ; 

end 

elf orce=k*eldisp; 

stress (iel) =sqrt (el force (1) A2+elforce (2) A2) /area 

if ( (x2-xl) *elforce (3) )<0; 

stress (iel) ^-stress  (iel) ; 

end 


end 

num=l : 1 : sdof ; 
displ=[num'  disp] 

numm=l : 1 : nel ; 
stresses^ [numm 1  stress] 
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FEMG5 


nel=13 ; 
nnel=2 ; 
ndof=2 ; 
nnode=8 ; 

sdof =nnode*ndof ; 

gcoord ( 1 , 1 ) =0 ; gcoord (1^2) =0 ; 
gcoord (2,1) =  0 ; gcoord (2, 2) =  6; 
gcoord (3,1) =0; gcoord (3, 2) =  12; 
gcoord (4,1)=6.365; gcoord (4 , 2 ) =12 ; 
gcoord (5,1)=6.365; gcoord (5, 2) =6; 
gcoord (6, I) =6 . 365; gcoord ( 6, 2 ) =0 ; 
gcoord (7 , 1 ) =3 . 1825 ; gcoord (7 , 2 ) =3; 
gcoord ( 8 , 1 ) =3 . 1825 ; gcoord ( 8 , 2 ) =9 ; 

elprop (1, 1)=10300000; 
elprop (1,2)=. 125; 
elprop (2, 1) =10300000; 
elprop  (2 , 2 )  =  . 12  5 ; 
elprop  (3,  1)=10300000; 
elprop (3,2)=. 125; 
elprop (4, 1) =10300000; 
elprop (4,2)=. 125; 
elprop (5, 1) =10300000; 
elprop (5 , 2 ) = . 125 ; 
elprop (6, 1) =10300000; 
elprop ( 6, 2 ) = . 125 ; 
elprop (7, 1) =10300000; 
elprop  (7 , 2 )  =  . 12  5 ; 
elprop  (8, 1)=10300000; 
elprop ( 8 , 2 ) = . 125 ; 
elprop (9, 1) =10300000; 
elprop  ( 9, 2 )  =  . 12  5 ; 
elprop (10, 1) =10300000 ; 
elprop (10,2)=. 125; 
elprop (11,1)=10300000; 
elprop (11, 2) =. 125; 
elprop (12, 1) =10300000 ; 
elprop (12,2)=. 125; 
elprop (13, 1) =10300000; 
elprop (13, 2) =. 125; 

nodes  (1,1)  =71 ;  nodes  (1,2)  =2 ; 
nodes (2,1) =2 ; nodes (2,2) =7 ; 
nodes (3,1) =2; nodes (3, 2) =8; 
nodes (4,1) =2 ; nodes (4,2) =3 ; 
nodes (5,1) =3; nodes (5, 2) =8; 
nodes (6, 1) =3; nodes (6,2) =4 ; 
nodes (7,1) =4; nodes (7, 2) =8; 
nodes (8, 1) =4; nodes (8,2) =5 ; 
nodes (9, 1) =5; nodes (9, 2) =8; 
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nodes ( 10 , 1 ) =5 ; nodes ( 10 , 2 ) =7 ; 
nodes ( 11 , 1 ) =5 ; nodes ( 11 , 2 ) =  6; 
nodes ( 12 , 1 ) =6; nodes ( 12 , 2 ) =7 ; 
nodes ( 13 , 1) =7 ; nodes (13, 2) =1/ 

bcdof (1) =1; 
bcval ( 1 ) =0 ; 
bcdof (2) =2; 
bcval (2 ) =0 ; 
bcdof (3) =11; 
bcval (3 ) =0 ; 
bcdof (4)=12; 
bcval ( 4 ) =0 ; 

ff=zeros (sdof, 1) ; 
kk=zeros (sdof ,  sdof ) ; 
index=zeros (nnel*ndof ,  1)  ; 
elf orce=zeros (nnel*ndof ,  1)  ; 
eldisp=zeros (nnel*ndof , nnel*ndof ) ; 
stress=zeros (nel, 1) ; 

ff  (6)=-600; 
ff (8) =-600; 

for  iel=l : nel 

nd (1) =nodes  (iel, 1) ; 
nd (2 ) =nodes (iel, 2 ) ; 

xl=gcoord (nd ( 1 ) ,  1 ) ; yl=gcoord  (nd  ( 1 ) ,  2 ) 
x2=gcoord (nd ( 2 ) , 1 ) ; y2=gcoord (nd ( 2 ) ,2) 

leng=sqrt ( (x2-xl) A2+ (y2-yl) A2) ; 

if  (x2-xl ) ==0 ; 

beta=2*atan (1) ; 

else 

beta=atan ( (y2-yl) / (x2-xl) ) ; 

end 

el=elprop  (iel, 1) ; 
area=elprop  ( iel , 2 ) ; 

[ index] =FEELDOF (nd, nnel , ndof) ; 

k=FETRUSS2 (el , leng, area, 0 , beta, 1 ) ; 

kk=FEASMBLl (kk, k, index)  ; 


end 

[ kk, f f ] =FEAPLYC2 ( kk, f f , bcdof, bcval ) ; 


disp=kk\f f ; 
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for  iel=l : nel 


nd (1) =nodes (iel, 1)  ; 
nd (2 ) =nodes (iel,  2 )  ; 

xl=gcoord (nd ( 1 ) ,  1 ) ; yl=gcoord (nd ( 1 )  ,  2 )  ; 
x2=gcoord (nd (2 ) ,  1 ) ; y2=gcoord (nd (2 )  ,2)  ; 

leng=sqrt ( (x2-xl) A2+ (y2-yl) A2) ; 

if  (x2-xl ) ==0 ; 

beta=2*atan (1) ; 

else 

beta=atan ( (y2-yl) / (x2-xl) ) ; 

end 

el=elprop  (iel, 1) ; 
area=elprop  (iel, 2 ) ; 
index=FEELDOF (nd, nnel , ndof ) ; 

k=FETRUSS2 (el, leng, area, 0 , beta, 1 ) ; 

for  i=l :  (nnel*ndof) 

eldisp (i) =disp (index (i) ) ; 

end 

elf orce=k*eldisp; 

stress (iel) =sqrt (el force (1) A2+elforce (2) A2) /area 

if ( (x2-xl) *elforce (3) )<0; 

stress (iel) ^-stress  (iel) ; 

end 


end 

num=l : 1 : sdof ; 
displ=[num'  disp] 

numm=l : 1 : nel ; 
stresses^ [numm 1  stress] 
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FEMG25 


nel=13 ; 
nnel=2 ; 
ndof=2 ; 
nnode=8 ; 

sdof =nnode*ndof ; 

gcoord ( 1 , 1 ) =0 ; gcoord (1^2) =0 ; 
gcoord (2,1) =  0 ; gcoord (2, 2) =  6; 
gcoord (3,1) =0; gcoord (3, 2) =  12; 
gcoord (4,1)=6.365; gcoord (4 , 2 ) =12 ; 
gcoord (5,1)=6.365; gcoord (5, 2) =6; 
gcoord (6, I) =6 . 365; gcoord ( 6, 2 ) =0 ; 
gcoord (7 , 1 ) =3 . 1825 ; gcoord (7 , 2 ) =3; 
gcoord ( 8 , 1 ) =3 . 1825 ; gcoord ( 8 , 2 ) =9 ; 

elprop (1, 1)=10300000; 
elprop (1,2)=. 0625; 
elprop (2, 1) =10300000; 
elprop  (2,2)=. 0625; 
elprop  (3,  1)=10300000; 
elprop (3,2)=.0625; 
elprop (4, 1) =10300000; 
elprop (4,2)=. 0625 ; 
elprop (5, 1) =10300000; 
elprop  (5, 2)  =  . 0625; 
elprop (6, 1) =10300000; 
elprop (6, 2) =. 0625; 
elprop (7, 1) =10300000; 
elprop  (7,2)=. 0625; 
elprop  (8, 1)=10300000; 
elprop  (8, 2)  =  . 0625; 
elprop (9, 1) =10300000; 
elprop  (9,2)  =  .0625; 
elprop (10, 1) =10300000 ; 
elprop (10, 2) =. 0625; 
elprop (11,1)=10300000; 
elprop (11, 2) =. 0625; 
elprop (12, 1) =10300000 ; 
elprop (12, 2) =. 0625; 
elprop (13, 1) =10300000 ; 
elprop (13, 2) =. 0625; 

nodes  (1,1)  =71 ;  nodes  (1,2)  =2 ; 
nodes (2,1) =2 ; nodes (2,2) =7 ; 
nodes (3,1) =2; nodes (3, 2) =8; 
nodes (4,1) =2 ; nodes (4,2) =3 ; 
nodes (5,1) =3; nodes (5, 2) =8; 
nodes (6, 1) =3; nodes (6,2) =4 ; 
nodes (7,1) =4; nodes (7, 2) =8; 
nodes (8, 1) =4; nodes (8,2) =5 ; 
nodes (9, 1) =5; nodes (9, 2) =8; 
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nodes ( 10 , 1 ) =5 ; nodes ( 10 , 2 ) =7 ; 
nodes ( 11 , 1 ) =5 ; nodes ( 11 , 2 ) =  6; 
nodes ( 12 , 1 ) =6; nodes ( 12 , 2 ) =7 ; 
nodes ( 13 , 1) =7 ; nodes (13, 2) =1/ 

bcdof (1) =1; 
bcval ( 1 ) =0 ; 
bcdof (2) =2; 
bcval (2 ) =0 ; 
bcdof (3) =11; 
bcval (3 ) =0 ; 
bcdof (4)=12; 
bcval ( 4 ) =0 ; 

ff=zeros (sdof, 1) ; 
kk=zeros (sdof ,  sdof ) ; 
index=zeros (nnel*ndof ,  1)  ; 
elf orce=zeros (nnel*ndof ,  1)  ; 
eldisp=zeros (nnel*ndof , nnel*ndof ) ; 
stress=zeros (nel, 1) ; 

ff  (6)=-600; 
ff (8) =-600; 

for  iel=l : nel 

nd (1) =nodes  (iel, 1) ; 
nd (2 ) =nodes (iel, 2 ) ; 

xl=gcoord (nd ( 1 ) ,  1 ) ; yl=gcoord  (nd  ( 1 ) ,  2 ) 
x2=gcoord (nd ( 2 ) , 1 ) ; y2=gcoord (nd ( 2 ) ,2) 

leng=sqrt ( (x2-xl) A2+ (y2-yl) A2) ; 

if  (x2-xl ) ==0 ; 

beta=2*atan (1) ; 

else 

beta=atan ( (y2-yl) / (x2-xl) ) ; 

end 

el=elprop  (iel, 1) ; 
area=elprop  ( iel , 2 ) ; 

[ index] =FEELDOF (nd, nnel , ndof) ; 

k=FETRUSS2 (el , leng, area, 0 , beta, 1 ) ; 

kk=FEASMBLl (kk, k, index)  ; 


end 

[ kk, f f ] =FEAPLYC2 ( kk, f f , bcdof, bcval ) ; 


disp=kk\f f ; 
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for  iel=l : nel 


nd (1) =nodes (iel, 1)  ; 
nd (2 ) =nodes (iel,  2 )  ; 

xl=gcoord (nd ( 1 ) ,  1 ) ; yl=gcoord (nd ( 1 )  ,  2 )  ; 
x2=gcoord (nd (2 ) ,  1 ) ; y2=gcoord (nd (2 )  ,2)  ; 

leng=sqrt ( (x2-xl) A2+ (y2-yl) A2) ; 

if  (x2-xl ) ==0 ; 

beta=2*atan (1) ; 

else 

beta=atan ( (y2-yl) / (x2-xl) ) ; 

end 

el=elprop  (iel, 1) ; 
area=elprop  (iel, 2 ) ; 
index=FEELDOF (nd, nnel , ndof ) ; 

k=FETRUSS2 (el, leng, area, 0 , beta, 1 ) ; 

for  i=l :  (nnel*ndof) 

eldisp (i) =disp (index (i) ) ; 

end 

elf orce=k*eldisp; 

stress (iel) =sqrt (el force (1) A2+elforce (2) A2) /area 

if ( (x2-xl) *elforce (3) )<0; 

stress (iel) ^-stress  (iel) ; 

end 


end 

num=l : 1 : sdof ; 
displ=[num'  disp] 

numm=l : 1 : nel ; 
stresses^ [numm 1  stress] 
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FEMH1 


nel=14 ; 
nnel=2 ; 
ndof=2 ; 
nnode=8 ; 

sdof =nnode*ndof ; 

gcoord ( 1 , 1 ) =0 ; gcoord ( 1 , 2 ) =0 ; 
gcoord (2,1) =  0 ; gcoord (2, 2) =6; 
gcoord (3,1) =0; gcoord (3, 2) =  12; 
gcoord (4,1)=6.365; gcoord (4 , 2 ) =12 ; 
gcoord (5,1)=6.365; gcoord (5, 2) =6; 
gcoord (6, I) =6 . 365; gcoord ( 6, 2 ) =0 ; 
gcoord (7 , 1 ) =3 . 1825 ; gcoord (7 , 2 ) =3; 
gcoord ( 8 , 1 ) =3 . 1825 ; gcoord ( 8 , 2 ) =9 ; 

elprop (1, 1)=10300000; 
elprop (1, 2) =.25; 
elprop  (2, 1) =10300000; 
elprop (2 , 2 ) = . 25; 
elprop (3, 1)=10300000; 
elprop  (3, 2)  =  .  25,- 
elprop  (4,  1)  =10300000; 
elprop (4, 2) =.25; 
elprop (5, 1) =10300000; 
elprop  (5 , 2 )  =  . 2  5 ; 
elprop (6, 1) =10300000; 
elprop ( 6, 2 ) = . 25; 
elprop (7, 1) =10300000; 
elprop  (7 , 2 )  =  . 25 ; 
elprop (8, 1)=10300000; 
elprop (8, 2) =.25; 
elprop  (9,  1) =10300000; 
elprop ( 9, 2 ) = . 25 ; 
elprop (10, 1) =10300000; 
elprop ( 10 , 2 ) = . 25 ; 
elprop (11, 1) =10300000; 
elprop ( 11 , 2 ) = . 25 ; 
elprop (12, 1) =10300000 ; 
elprop ( 12 , 2 ) = . 25 ; 
elprop (13, 1) =10300000 ; 
elprop ( 13 , 2 ) = . 25 ; 
elprop (14, 1) =10300000 ; 
elprop (14,2)=. 25; 

nodes ( 1 , 1 ) =1 ; nodes (1,2) =2 ; 
nodes (2,1) =2 ; nodes (2,2) =7 ; 
nodes (3,1) =2 ; nodes (3, 2) =5; 
nodes (4,1) =2 ; nodes (4, 2) =8; 
nodes (5,1) =2 ; nodes (5,2) =3 ; 
nodes (6, 1) =3; nodes (6, 2) =8; 
nodes (7,1) =3; nodes (7,2) =4 ; 
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nodes ( 8 , 1 ) =4 ; nodes ( 8 , 2 ) =8 ; 
nodes ( 9, 1 ) =4 ; nodes ( 9, 2 ) =5/ 
nodes ( 10 , 1 ) =5 ; nodes (10r  2) =8; 
nodes ( 11 , 1 ) =5 ; nodes ( 11 , 2 ) =7 ; 
nodes ( 12 , 1 ) =5 ; nodes ( 12 , 2 ) =6; 
nodes ( 13 , 1 ) =  6; nodes ( 13 , 2 ) =7 ; 
nodes ( 14 , 1 ) =7 ; nodes ( 14 , 2 ) =1 ; 

bcdof (1) =1; 
bcval ( 1 ) =0 ; 
bcdof (2) =2; 
bcval (2 ) =0 ; 
bcdof (3) =11; 
bcval ( 3 ) =0 ; 
bcdof (4) =12; 
bcval ( 4 ) =0 ; 

ff=zeros (sdof, 1) ; 
kk=zeros (sdof , sdof ) ; 
index=zeros (nnel^ndof ,  1)  ; 
elf orce=zeros (nnel*ndof ,  1)  ; 
eldisp=zeros (nnel*ndof , nnel^ndof )  ; 
stress=zeros (nel,  1)  ; 

ff (6) =-600; 
ff (8)=-600; 

for  iel=l : nel 

nd ( 1 ) =nodes ( iel ,  1 )  ; 
nd (2 ) =nodes (iel,  2 )  ; 

xl=gcoord (nd ( 1 ) , 1 ) ; yl=gcoord (nd ( 1 )  ,2) 
x2=gcoord (nd (2 ) , 1 ) ; y2=gcoord (nd (2 ) ,2) 

leng=sqrt ( (x2-xl) A2+ (y2-yl) A2) ; 

if  (x2-xl ) ==0 ; 

beta=2^atan (1) ; 

else 

beta=atan ( (y2-yl) / (x2-xl) ) ; 

end 

el=elprop  (iel, 1) ; 
area=elprop  ( iel , 2 ) ; 

[ index] =FEELDOF (nd, nnel , ndof) ; 

k=FETRUSS2 (el , leng, area, 0 , beta, 1 ) ; 

kk=FEASMBLl (kk, k, index)  ; 


end 


115 


[kk, f f ] =FEAPLYC2 (kk,  f  f , b c do  f  , bcval) ; 


disp=kk\f f ; 

for  iel=l : nel 

nd (1) =nodes  (iel, 1) ; 
nd (2 ) =nodes (iel, 2 ) ; 

xl=gcoord (nd ( 1 ) ,  1 ) ; yl=gcoord (nd  ( 1 ) ,  2 )  ; 
x2=gcoord (nd (2 ) ,  1 ) ; y2=gcoord (nd (2 ) ,2) ; 

leng=sqrt ( (x2-xl) A2+ (y2-yl) A2) ; 

if  (x2-xl ) ==0 ; 

beta=2*atan (1) ; 

else 

beta=atan ( (y2-yl) / (x2-xl) ) / 

end 

el=elprop  (iel, 1) ; 
area=elprop  (iel, 2 ) ; 
index=FEELDOF (nd, nnel , ndof ) ; 

k=FETRUSS2 (el, leng, area, 0 , beta, 1 ) ; 

for  i=l :  (nnel*ndof) 

eldisp (i) =disp (index (i) ) ; 

end 

elf orce=k*eldisp; 

stress (iel) =sqrt (elforce (1) A2+elforce (2) A2) /area 

if ( (x2-xl ) *elforce(3) )<0; 

stress (iel) ^-stress  (iel) ; 

end 


end 

num=l : 1 : sdof ; 
displ=[num'  disp] 

numm=l : 1 : nel / 
stresses^ [numm 1  stress] 
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FEMH5 


nel=14 ; 
nnel=2 ; 
ndof=2 ; 
nnode=8 ; 

sdof =nnode*ndof ; 

gcoord ( 1 , 1 ) =0 ; gcoord ( 1 , 2 ) =0 ; 
gcoord (2,1) =  0 ; gcoord (2, 2) =6; 
gcoord (3,1) =0; gcoord (3, 2) =  12; 
gcoord (4,1)=6.365; gcoord (4 , 2 ) =12 ; 
gcoord (5,1)=6.365; gcoord (5, 2) =6; 
gcoord (6, I) =6 . 365; gcoord ( 6, 2 ) =0 ; 
gcoord (7 , 1 ) =3 . 1825 ; gcoord (7 , 2 ) =3; 
gcoord ( 8 , 1 ) =3 . 1825 ; gcoord ( 8 , 2 ) =9 ; 

elprop (1, 1)=10300000; 
elprop (1,2)=. 125; 
elprop  (2, 1) =10300000; 
elprop (2 , 2 ) = . 125 ; 
elprop (3, 1)=10300000; 
elprop  (3,2)=. 125; 
elprop  (4, 1) =10300000; 
elprop (4,2)=. 125; 
elprop (5, 1) =10300000; 
elprop (5 , 2 ) = . 125 ; 
elprop (6, 1) =10300000; 
elprop ( 6, 2 ) = . 125 ; 
elprop (7, 1) =10300000; 
elprop (7 , 2 ) = . 125 ; 
elprop (8, 1)=10300000; 
elprop  ( 8 , 2 )  =  . 12  5 ; 
elprop (9, 1) =10300000; 
elprop  ( 9, 2 )  =  . 12  5 ; 
elprop (10, 1) =10300000 ; 
elprop (10,2)=. 125; 
elprop (11,1)=10300000; 
elprop (11, 2) =. 125; 
elprop (12, 1) =10300000 ; 
elprop (12,2)=. 125; 
elprop (13,1)=10300000; 
elprop (13,2)=. 125; 
elprop (14, 1) =10300000; 
elprop (14, 2) =. 125; 

nodes ( 1 , 1 ) =1 ; nodes (1,2) =2 ; 
nodes (2,1) =2 ; nodes (2,2) =7 ; 
nodes (3,1) =2 ; nodes (3, 2) =5; 
nodes (4,1) =2 ; nodes (4, 2) =8; 
nodes (5,1) =2 ; nodes (5,2) =3 ; 
nodes (6, 1) =3; nodes (6, 2) =8; 
nodes (7,1) =3; nodes (7,2) =4 ; 
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nodes ( 8 , 1 ) =4 ; nodes ( 8 , 2 ) =8 ; 
nodes ( 9, 1 ) =4 ; nodes ( 9, 2 ) =5/ 
nodes ( 10 , 1 ) =5 ; nodes (10r  2) =8; 
nodes ( 11 , 1 ) =5 ; nodes ( 11 , 2 ) =7 ; 
nodes ( 12 , 1 ) =5 ; nodes ( 12 , 2 ) =6; 
nodes ( 13 , 1 ) =  6; nodes ( 13 , 2 ) =7 ; 
nodes ( 14 , 1 ) =7 ; nodes ( 14 , 2 ) =1 ; 

bcdof (1) =1; 
bcval ( 1 ) =0 ; 
bcdof (2) =2; 
bcval (2 ) =0 ; 
bcdof (3) =11; 
bcval ( 3 ) =0 ; 
bcdof (4) =12; 
bcval ( 4 ) =0 ; 

ff=zeros (sdof, 1) ; 
kk=zeros (sdof , sdof ) ; 
index=zeros (nnel^ndof ,  1)  ; 
elf orce=zeros (nnel*ndof ,  1)  ; 
eldisp=zeros (nnel*ndof , nnel^ndof )  ; 
stress=zeros (nel,  1)  ; 

ff (6) =-600; 
ff (8)=-600; 

for  iel=l : nel 

nd ( 1 ) =nodes ( iel ,  1 )  ; 
nd (2 ) =nodes (iel,  2 )  ; 

xl=gcoord (nd ( 1 ) , 1 ) ; yl=gcoord (nd ( 1 )  ,2) 
x2=gcoord (nd (2 ) , 1 ) ; y2=gcoord (nd (2 ) ,2) 

leng=sqrt ( (x2-xl) A2+ (y2-yl) A2) ; 

if  (x2-xl ) ==0 ; 

beta=2^atan (1) ; 

else 

beta=atan ( (y2-yl) / (x2-xl) ) ; 

end 

el=elprop  (iel, 1) ; 
area=elprop  ( iel , 2 ) ; 

[ index] =FEELDOF (nd, nnel , ndof) ; 

k=FETRUSS2 (el , leng, area, 0 , beta, 1 ) ; 

kk=FEASMBLl (kk, k, index)  ; 


end 
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[kk, f f ] =FEAPLYC2 (kk,  f  f , b c do  f  , bcval) ; 


disp=kk\f f ; 

for  iel=l : nel 

nd (1) =nodes  (iel, 1) ; 
nd (2 ) =nodes (iel, 2 ) ; 

xl=gcoord (nd ( 1 ) ,  1 ) ; yl=gcoord (nd  ( 1 ) ,  2 )  ; 
x2=gcoord (nd (2 ) ,  1 ) ; y2=gcoord (nd (2 ) ,2) ; 

leng=sqrt ( (x2-xl) A2+ (y2-yl) A2) ; 

if  (x2-xl ) ==0 ; 

beta=2*atan (1) ; 

else 

beta=atan ( (y2-yl) / (x2-xl) ) / 

end 

el=elprop  (iel, 1) ; 
area=elprop  (iel, 2 ) ; 
index=FEELDOF (nd, nnel , ndof ) ; 

k=FETRUSS2 (el, leng, area, 0 , beta, 1 ) ; 

for  i=l :  (nnel*ndof) 

eldisp (i) =disp (index (i) ) ; 

end 

elf orce=k*eldisp; 

stress (iel) =sqrt (elforce (1) A2+elforce (2) A2) /area 

if ( (x2-xl ) *elforce(3) )<0; 

stress (iel) ^-stress  (iel) ; 

end 


end 

num=l : 1 : sdof ; 
displ=[num'  disp] 

numm=l : 1 : nel / 
stresses^ [numm 1  stress] 
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FEMH25 


nel=14 ; 
nnel=2 ; 
ndof=2 ; 
nnode=8 ; 

sdof =nnode*ndof ; 

gcoord ( 1 , 1 ) =0 ; gcoord ( 1 , 2 ) =0 ; 
gcoord (2,1) =  0 ; gcoord (2, 2) =6; 
gcoord (3,1) =0; gcoord (3, 2) =  12; 
gcoord (4,1)=6.365; gcoord (4 , 2 ) =12 ; 
gcoord (5,1)=6.365; gcoord (5, 2) =6; 
gcoord (6, I) =6 . 365; gcoord ( 6, 2 ) =0 ; 
gcoord (7 , 1 ) =3 . 1825 ; gcoord (7 , 2 ) =3; 
gcoord ( 8 , 1 ) =3 . 1825 ; gcoord ( 8 , 2 ) =9 ; 

elprop (1, 1)=10300000; 
elprop (1,2)=. 0625; 
elprop  (2, 1) =10300000; 
elprop (2,2)=. 0625; 
elprop (3, 1)=10300000; 
elprop  (3,2)=.0625; 
elprop  (4, 1) =10300000; 
elprop  (4,2)  =  . 0625 ; 
elprop (5, 1) =10300000; 
elprop (5, 2) =. 0625; 
elprop (6, 1) =10300000; 
elprop  (6, 2)  =  . 0625; 
elprop (7, 1) =10300000; 
elprop (7,2)=.0625; 
elprop (8, 1)=10300000; 
elprop  (8, 2)  =  . 0625; 
elprop (9, 1) =10300000; 
elprop  (9,2)  =  .0625; 
elprop (10, 1) =10300000 ; 
elprop (10, 2) =. 0625; 
elprop (11,1)=10300000; 
elprop (11, 2) =. 0625; 
elprop (12, 1) =10300000 ; 
elprop  (12, 2)  =  . 0625; 
elprop (13, 1) =10300000 ; 
elprop (13, 2) =. 0625; 
elprop (14, 1) =10300000; 
elprop (14, 2) =. 0625; 

nodes ( 1 , 1 ) =1 ; nodes (1,2) =2 ; 
nodes (2,1) =2 ; nodes (2,2) =7 ; 
nodes (3,1) =2 ; nodes (3, 2) =5; 
nodes (4,1) =2 ; nodes (4, 2) =8; 
nodes (5,1) =2 ; nodes (5,2) =3 ; 
nodes (6, 1) =3; nodes (6, 2) =8; 
nodes (7,1) =3; nodes (7,2) =4 ; 
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nodes ( 8 , 1 ) =4 ; nodes ( 8 , 2 ) =8 ; 
nodes ( 9, 1 ) =4 ; nodes ( 9, 2 ) =5/ 
nodes ( 10 , 1 ) =5 ; nodes (10r  2) =8; 
nodes ( 11 , 1 ) =5 ; nodes ( 11 , 2 ) =7 ; 
nodes ( 12 , 1 ) =5 ; nodes ( 12 , 2 ) =6; 
nodes ( 13 , 1 ) =  6; nodes ( 13 , 2 ) =7 ; 
nodes ( 14 , 1 ) =7 ; nodes ( 14 , 2 ) =1 ; 

bcdof (1) =1; 
bcval ( 1 ) =0 ; 
bcdof (2) =2; 
bcval (2 ) =0 ; 
bcdof (3) =11; 
bcval ( 3 ) =0 ; 
bcdof (4) =12; 
bcval ( 4 ) =0 ; 

ff=zeros (sdof, 1) ; 
kk=zeros (sdof , sdof ) ; 
index=zeros (nnel^ndof ,  1)  ; 
elf orce=zeros (nnel*ndof ,  1)  ; 
eldisp=zeros (nnel*ndof , nnel^ndof )  ; 
stress=zeros (nel,  1)  ; 

ff (6) =-600; 
ff (8)=-600; 

for  iel=l : nel 

nd ( 1 ) =nodes ( iel ,  1 )  ; 
nd (2 ) =nodes (iel,  2 )  ; 

xl=gcoord (nd ( 1 ) , 1 ) ; yl=gcoord (nd ( 1 )  ,2) 
x2=gcoord (nd (2 ) , 1 ) ; y2=gcoord (nd (2 ) ,2) 

leng=sqrt ( (x2-xl) A2+ (y2-yl) A2) ; 

if  (x2-xl ) ==0 ; 

beta=2^atan (1) ; 

else 

beta=atan ( (y2-yl) / (x2-xl) ) ; 

end 

el=elprop  (iel, 1) ; 
area=elprop  ( iel , 2 ) ; 

[ index] =FEELDOF (nd, nnel , ndof) ; 

k=FETRUSS2 (el , leng, area, 0 , beta, 1 ) ; 

kk=FEASMBLl (kk, k, index)  ; 


end 
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[kk, f f ] =FEAPLYC2 (kk,  f  f , b c do  f  , bcval) ; 


disp=kk\f f ; 

for  iel=l : nel 

nd (1) =nodes  (iel, 1) ; 
nd (2 ) =nodes (iel, 2 ) ; 

xl=gcoord (nd ( 1 ) ,  1 ) ; yl=gcoord (nd  ( 1 ) ,  2 )  ; 
x2=gcoord (nd (2 ) ,  1 ) ; y2=gcoord (nd (2 ) ,2) ; 

leng=sqrt ( (x2-xl) A2+ (y2-yl) A2) ; 

if  (x2-xl ) ==0 ; 

beta=2*atan (1) ; 

else 

beta=atan ( (y2-yl) / (x2-xl) ) / 

end 

el=elprop  (iel, 1) ; 
area=elprop  (iel, 2 ) ; 
index=FEELDOF (nd, nnel , ndof ) ; 

k=FETRUSS2 (el, leng, area, 0 , beta, 1 ) ; 

for  i=l :  (nnel*ndof) 

eldisp (i) =disp (index (i) ) ; 

end 

elf orce=k*eldisp; 

stress (iel) =sqrt (elforce (1) A2+elforce (2) A2) /area 

if ( (x2-xl ) *elforce(3) )<0; 

stress (iel) ^-stress  (iel) ; 

end 


end 

num=l : 1 : sdof ; 
displ=[num'  disp] 

numm=l : 1 : nel / 
stresses^ [numm 1  stress] 
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APPENDIX  B 


H0L2F 

H  05" 

HI" 
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64000  -0.0047 
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7  otxx)  aoooo 
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7.0000  OLOOOO 

84000  -0.0092 

84000  -0.0046 
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84000  -0.0056 
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84000  440B6 
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94000  04026 
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94000  OLOOOO 

104000  44045 
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304000  40014 

300000  -0.0093 

10L0000  -0.0047 

laoooo  -0.0023 

114000  44000 

114000  44000 

114000  44000 

114000  010000 

110000  04000 

110000  04000 

11 0000  04005 

110000  04003 

110000  a  00m 

124000  04000 

124000  04000 

124000  04000 

124000  40000 

12.0000  40000 

124000  40000 

12.0000  -00056 

124000  44028 

124000  44014 

134000  44000 

134000  44000 

134000  44000 

134000  aoooo 

130000  04000 

134000  04000 

13LOOOO  04000 

134000  04000 

134000  nmm 

144000  44017 

144000  4400B 

144000  40004 

144000  40000 

14.0000  40000 

144000  40000 

14.0000  04000 

144000  04000 

144000  aoooo 

nam  44000 

iqnmn  44000 

mmm  -nmm 

154000  0 

15l0000  0 

154000  0 

15.0000  nmm 

wmm  nmm 

154000  nmm 

160000  44072 

164000  44036 

164000  44012 

juGuuG  -0.0112 

l&SHB  -0.0056 

JULUUUU  -0.0028 

laoooo  oloooo 

16  nmn  OLOOOO 

164000  aoooo 

stresses  = 

stresses = 

stresses  = 

stresses  = 

stresses = 

stresses  = 

streses= 

stresses  = 

stresses  = 

14fe-40B* 

10*409* 

10*409* 

14fe-403" 

10*409* 

LOb-rfXB" 

LOe+OOB* 

10*409* 

10*409* 

04010  7.7597 

04010  1879B 

0L0Q1D  19999 

04010  9.6000 

04030  4.8000 

04010  2.4000 

04010  9.6000 

04010  4.8C00 

04010  2.4000 

00020  -7.GR30 

04070  -13415 

a  ram  46*i7 

nmm  omm 

mnom  omm 

omm  (uixi) 

nmm  14973 

nmm  08487 

nmm  04743 

04090  16564 

04030  L8282 

010090  09141 

04090  40000 

04090  -aoooo 

04090  44000 

04030  36223 

04090  -18111 

00090  49056 

04040  2J419 

04040  11709 

00040  45855 

04040  9.6000 

04040  4.8000 

04040  2.4000 

04040  64000 

04040  12000 

04040  16000 

04050  7.9936 

04050  3.9968 

a 0050  1.998^1 

04050  0 

04050  0 

04050  0 

04050  44000 

04050  40000 

nonm  40000 

00060  23419 

04060  -LOT 

a 0060  45855 

04060  04000 

04060  04000 

04060  04000 

00060  36223 

04060  -18111 

04060  49066 

04070  L7D41 

04070  08521 

a  0070  a42G0 

04070  40000 

04070  -aoooo 

04070  44000 

04070  44000 

04070  40000 

04070  -OLOOOO 

04080  -23419 

04080  -L1OT 

a 0080  45855 

04080  9.6000 

04080  4.8000 

04080  2.4000 

04080  64000 

04080  12000 

04080  16000 

04090  7.9936 

04090  3.9968 

a 0090  1.9984 

04090  04000 

04090  04000 

04090  aoooo 

04090  16973 

04090  08487 

04090  04243 

04100  -23419 

04100  -11309 

aoioo  45855 

00100  40000 

aoioo  -aoooo 

04100  44000 

CLOUD  9.6000 

04100  4.8C00 

04300  2.4000 

nmm  Z6830 

04110  -13415 

a  QUO  46707 

nmiD  9.6000 

04110  4.8000 

04110  2.4000 

00130  44000 

nmm  40000 

nmm  44000 

00120  7.7597 

00120  3L879B 

a012D  L9999 

nmm  04000 

04120  04000 

nmm  aoooo 

nmm  n nmn 

nmm  oloooo 

nmm  aoooo 

00130  -Z6B30 

00130  -13415 

a 0130  46707 

04130  40000 

00130  40000 

04330  44000 

00140  -2j6B9D 

00140  -13415 

a Q140  46707 
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COJS* 


D  0L2P 

DOJ 

D  1* 

cfhpl- 

dhpl- 

dtapl- 

LOOM  -0.0000 

L0000  -0X000 

LOOOO  -0X000 

2-0000  0-0000 

2_0000  0.0000 

2-0000  OlOOOO 

3  0000  0.0000 

3,0000  0X000 

9-0000  OLOOOO 

4-0000  -00094 

4X000  -0  0Q17 

4X000  -0X008 

■umo  -n  firms 

5-0000  -0X009 

5.0000  -0  0002 

umo  -0.0068 

6X000  -0.0034 

6X000  -0.0017 

7X1000  nf¥WK 

7X000  0,0009 

7X000  0X002 

LOOOO  -0.0068 

8X000  -0.0034 

8X000  -0.0017 

9.0000  0-0000 

9  0000  0.0000 

910000  OLOOOO 

10.0000  -OXOM 

10.0000  -0-0017 

jnnnnp  -0  QPQB 

11.0000  0X000 

ii  nrmn  QjDOQQ 

ii  nrmn  (mOQQ 

12.0000  0 

12X000  0 

12X000  0 

13.0000  0X000 

13X000  0X000 

13X000  OLOOOO 

14.0000  -0L0CH2 

14X000  -0X016 

140000  -0X008 

stresses  = 

stresses = 

stresses  = 

lXe+C03* 

IXchOCB* 

lXetOOB  * 

O.OOLO  5.7981 

0X010  2.8990 

OlOOlD  1.4495 

0X020  0 

0X020  0 

0X020  0 

0X090  5.7981 

0.0030  2.8990 

00(00  1.4495 

0X040  -43096 

0X040  -2.1518 

0X040  -1X759 

0X050  2X166 

0X050  IXOfB 

0X050  05041 

0X060  -43096 

0X060  -2.1518 

00060  -1X759 

0X070  5.7981 

0X070  2.8990 

0X070  1.4495 

0X080  0 

0X080  0 

0X080  0 

0X090  5.7981 

0.0090  2.8990 

0X090  1.4495 

0X100  -43036 

0X100  -2.1518 

0X100  -1X759 

nmin  -43036 

0X110  -2.1518 

0X110  -1X759 

P  P  h  6  P  P  9  •  11  iPbioAiPii 


COJ 

CP 

BOJP 

BQJ 

Blv 

cflspl- 

dhpl- 

dspl- 

dhpl- 

dbpl- 

1X000  -0-0000 

1X000  -QJQOOO 

1X000  0 

1X000  0 

1X000  0 

2_0000  -0,0000 

2-0000  -0X000 

2X000  -OlOOOO 

i  nrmn  -Q  QPOO 

2-0000  -0X000 

3.0000  010026 

50000  0XQ13 

3-0000  -0X053 

3J00Q  -0XQ26 

90000  -0X013 

4X000  -0X028 

4X000  -00014 

4-0000  -0X056 

4X000  -0X028 

4X000  -00014 

5.0000  -0-0000 

50000  -0-0000 

5-0000  -0X105 

5X000  -0X053 

5X000  -0X036 

6X000  -0.0056 

6X000  -0.0028 

6X000  -0.0112 

6X000  -0.0056 

6X000  -0.0028 

7.0000  -0-0000 

7-0000  -0X000 

7X000  -0X105 

7X000  -0X053 

7X000  -00026 

8X000  -0.0056 

8X000  -0.0028 

8X000  -0.0112 

8X000  -0.0056 

8X000  -0.0028 

9X000  00026 

4X000  nnma 

9  0000  -0X053 

9-0000  -0X026 

9X000  -00013 

10X000  -0X028 

-0  0014 

10-0000  -0.0056 

IPOPPO  _n tmm 

IQPfmp  -QX014 

11X000  0 

11X000  0 

11X000  0 

11X000  0 

11X000  0 

12X000  0 

12X000  0 

12X000  0 

12X000  0 

12X000  0 

stresses  = 

stresses = 

stresses 

stresses  = 

stresses  = 

lXe*003* 

LOe+OCB* 

UOe+OQ3* 

LOe-KXB  * 

lXe+XCB* 

0X010  0X000 

00010  0.0000 

0X010  9.6000 

0X010  4.8000 

0X010  2.4000 

0X020  4.8000 

0X020  2.4000 

0X020  0X000 

OX02D  0X000 

0X020  0X000 

0X030  0 

0X090  0 

0  0050  -OLOOOO 

0X030  -OXOOO 

0L0090  -OXOOO 

0X040  4.8000 

0X040  2.4000 

0X040  9.6000 

0X040  4.8000 

0X040  2.4000 

o  0050  -OXOOO 

OlOOSO  -0X000 

nooso  0-0000 

PQQ«^p  OlOOOO 

OLOOSO  0-0000 

0-0060  -OXOOO 

Q.0060  -OXOOO 

0X060  0 

0X060  0 

0X060  0 

0X070  4.8000 

0X070 

2.4000 

0X070  9.6000 

0X070  4.8000 

0X070  2.4000 

0X080  4.8000 

0X080 

2.4000 

0X080  9.6000 

0X080  4.8000 

0X080  2.4000 
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F  025T 

FOJ 

FI" 

A  027 

A  05* 

Al" 

dapl. 

dhpl- 

dipl- 

dbpl- 

dtapl- 

dtapl- 

i  fwwm  -ooqoq 

I  pppp  -ppppo 

I  fifing  -Qpppp 

1X000  0 

1X000  0 

1X000  0 

i  rmn  Q^popp 

2,0000  QJOOPO 

2.000Q  ojoooo 

2000Q  0X000 

20000  0X000 

20000  Qoooo 

5JPPP0  -QQPPff 

5.0000  -0X009 

5-0000  -00002 

9X000  -0.0211 

9X000  -0.0105 

9X000  -0.0053 

4JOOOO  --0.0068 

4J0D0D  -0.0034 

4X000  -aooi7 

4X000  -0X112 

4^0000  -QjQQSff 

4^0000  -nnnaa 

jpppp  QLOOOtr 

5X000  nnnfn 

5X000  00007 

5X000  -0.0211 

5X000  -0.0105 

5X000  -0.0053 

6JD000  --0.0068 

6X000  -0.0034 

6X000  -0L0O17 

6X000  -0X112 

OlQOPO  -QjQQSff 

6.0000  -OOfW 

7JDOOO  0 

7X000  0 

7X000  0 

7X000  0 

7X000  0 

7X000  0 

8JOOOO  0 

8X000  0 

8X000  0 

8X000  0 

8X000  0 

8X000  0 

ppppp  QlOQPP 

90000  OJQOOO 

90000  QjQOOO 

innnnn  .nnngn 

IDjOOPP  -0X016 

1PJPPPP  -0X008 

itjonea  ■ 

1mm 

llBUi  ■ 

Dm- 

1IMXB* 

ucmxb* 

KMX »• 

U0a*<xa* 

LOetOOS* 

1XM0B* 

GLOOUO  9.6000 

nnmn  4.8000 

OOOLO  2.400C 

finnan  0[0Q0Q 

nnnan  OLQpQp 

nnmn  n nmn 

nnmn  5.7981 

0X010  2.8990 

0.0010  14495 

ppppO  qloqqp 

00090  ojdoo 

0^0090  ooooo 

nnmn  -43CB6 

0X000  -2.1518 

0X020  -1X759 

00040  9.6000 

0X040  4.8000 

00040  2.400C 

nnnan  2X166 

nnnan  l  nan 

0X090  05041 

OJQQHQ  4  909b 

0X040  2.1518 

0X040  1X759 

OJOOBO  5.7981 

0X050  2.8990 

ouoon  14495 

nnnan  -A  anas 

0X060  -2.1518 

o  oom  -1X759 

0X070  -4309b 

0X070  -2.1518 

0X070  -1X759 
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APPENDIX  C 


rjWS. ) swjwtte  ‘rot*  ,'Jf 


Malarial 

Coat  (0/hgl 

Relative 

Cost 

[Ekra/bi*S] 

■  l«P»l 
[10*9  pal] 

PolaaonTa 
retie  H 

«Y[Mpa)[kil] 

*t  (Mpa)  [leal] 

K 

elongation 

itsal  Alley  ASS,  plst% '  <  rolled 

0.SMM 

1 

0.3 

'  I 

23 

V-jjr  ateel  Alky  i  y  ✓ 

1.1M.3Q 

4 

Oaf 

40 

Call  Iran  Gl  BOO 

1.90 

4.7 

SB 

Call  tom  6ri  id*  SO-4^11 

I.M-IjOO 

1.9 

0.19 

39SHHI 

t§ 

5ilr,  :  Alley  7373,  /  y  Ti 

om.to 

13.4 

OafS 

S ' 

11 

Nun  Inurn  Alley  3S6,  eutteni  tm  TO 

11.00 

m 

SSRSHi 

9.33 

EltBMMi 

iiUtDB 

1.3 

f  .flip:  01720%  s*  Mr  S30C  age 

Hi 

■ 

1,29  [0.2f8] 

i2i  ii  / 

u 

063-1209  [i  > 
17*1 

1140-1110  [IS  - 
1*0] 

Mugiseeiussi  AIieyA£310,  iheet 

11.00 

E£QM 

1.77 102)199] 

SRH 

ui 

220  [32] 

9Rnn 

u 

Titanium  Alley,  Commercial  purity  A3TM 
gnm  i 

mmmm 

■ 

4.31  fl}.;  S] 

109  [14.9] 

0.34 

170 [231 

140  [913 

90 

Geld  [cold  worked) 

14100 

19.12  [0>  19?’] 

77  [11.29 

M2 

ifIBWi 

lSTSHMI 

4 

Plarinm  [cold  worked) 

11,400-1*400 

11.43  [0  w 

171124.9] 

0.39 

SE33EJB 

1-3 

Site  [cold  worked) 

»V'>'  / 

271 

10.49  [0.3791 

0.37 

lf8Ml _ 

19 

NkkeliOO 

19-29.00 

sl.4 

9.39  [0.321] 

0.31 

143  [21.9] 

47 

Monel  400 

LS.90-LS.90 

22.9 

B.90E0.319] 

E 

0.32 

'40  >5’ 

KHI 

40 

Tip,  jtB'jw'ji  pirltf  {* 

5.BS-B.B3 

11.2 

m% 

IBH 

17 

L20 

la? 

104.9  [If  .2] 

o.ai 

13+199  [1949231 

SMI 
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APPENDIX  D 


M_Bastidty 
m  (E*6  psl) 

Density  (lbs/lnA3) 

Truss 

Width 

(In) 

Mass  of 
Full  Wall 
(ln*3) 

A  (lnA2) 

Volume 
of  Voids 
(ln*3) 

VPR(X) 

Mass 

One  Wall 
(lbs) 

•  max 

(psl) 

dmax 

(In) 

F=a*A 

(Ibsf) 

k=F/d 

(Ibs/in) 

32.25 

038045 

025 

37.88679 

0.0625 

17.40453 

0933471 

26.61 

5798.1 

00022 

362.3813 

164718.8 

32.25 

038045 

0625 

37.88679 

015625 

14.98602 

0784814 

28.18 

2319.2 

0000663 

362.375 

419901.5 

47,56104161 

0568671756 

0847976 

56.6307 

0213994 

13.61462 

0712994 

43.45 

1709.4 

0000431 

362.3827 

840210.3 

16.93895839 

0392228244 

0402024 

39.14289 

0300506 

1680712 

0859237 

13.77 

3605.6 

P 

1 

362.3841 

139378.5 

47  56104161 

0568671756 

0402024 

56.6307 

0.100506 

1680712 

0859237 

40.75 

3605.6 

fMMMEEL 

362.3841 

398355.6 

32.25 

038045 

0625 

37.88679 

015625 

34.98602 

0784834 

28.18 

2319.2 

0000663 

362.375 

419901.5 

16.93895839 

0.192228244 

0847976 

19.14289 

0233994 

13.61462 

0712994 

14.69 

1709.4 

00012 

362.3827 

301985.6 

32.25 

QJ8D45 

1 

37.88679 

025 

32.70833 

0665521 

29.66 

1449.5 

0000539 

362.375 

671811.3 

32.25 

0697 

0625 

6941016 

015625 

34.98602 

0784814 

51.63 

2319.2 

PQ00B63 

362.375 

419901.5 

47  56104161 

0392228244 

0847976 

19.14289 

0213994 

33.61462 

0712994 

14.69 

1709.4 

0000431 

362.3827 

840210.3 

16.93895839 

0568671756 

0847976 

56.6307 

0213994 

33.61462 

0712994 

43.45 

1709.4 

00012 

362.3827 

301985.6 

32.25 

038045 

0625 

37.88679 

015625 

34.98602 

0784814 

28.18 

2319.2 

0000663 

362.375 

419901.5 

lfi.93895839 

0568671756 

0402024 

56.6307 

0.100506 

16.40712 

0859237 

40.75 

3605.6 

P 

1 

362.3641 

139378.5 

32.25 

00639 

0.625 

6363428 

0.15625 

1438602 

0.784814 

4.73 

2319.2 

0000063 

362.375 

419901.5 

6.5 

038045 

0.625 

3738679 

0.15625 

14.98602 

0.784814 

28.18 

2319.2 

00043 

362.375 

84273.26 

32.25 

038045 

0.625 

3738679 

0.15625 

14.98602 

0.784814 

28.18 

2319.2 

0.000063 

362.375 

419901.5 

47  56104161 

0.192228244 

0402024 

19.14289 

0.100506 

16.40712 

0859237 

13.77 

3605.6 

0.00091 

362.3841 

398355.6 

32.25 

038045 

0625 

37,88679 

015625 

34.90602 

0784814 

28.18 

2319.2 

0000863 

362.375 

419901.5 

32.25 

038045 

0625 

3788679 

015625 

34.98602 

0784814 

28.18 

2319.2 

0000863 

362.375 

419901.5 

58 

038045 

0625 

3788679 

015625 

14.98602 

0784814 

28.18 

2319.2 

000048 

362.375 

755262.6 
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APPENDIX  E 


FEMFla 

nel=7 ; 
nnel=2 ; 
ndof =2 ; 
nnode=5 ; 

sdof=nnode*ndof ; 

gcoord ( 1 , 1 ) =0 ; gcoord ( 1 , 2 ) =0 ; 
gcoord (2 , 1 ) =0 ; gcoord (2 , 2 ) =12 ; 
gcoord (3,1) =  6. 3 65; gcoord (3, 2) =12; 
gcoord (4,1)=6.365; gcoord (4, 2) =0; 
gcoord (5, 1) =3 . 1825; gcoord (5, 2) =  6; 

elprop (1, 1)=10300000; 
elprop ( 1 , 2 ) = . 25 ; 
elprop (2, 1) =10300000; 
elprop (2 , 2 ) = . 25 ; 
elprop (3, 1)=10300000; 
elprop  (3 , 2 )  =  . 2  5 ; 
elprop (4, 1) =10300000; 
elprop (4 , 2 ) = . 25 ; 
elprop (5, 1) =10300000; 
elprop  (5, 2)  =  . 25; 
elprop (6, 1)=10300000; 
elprop ( 6, 2 ) = . 25; 
elprop  (7, 1) =10300000; 
elprop (7 , 2 ) = . 25 ; 

nodes ( 1 , 1 ) =1 ; nodes ( 1 , 2 ) =2 ; 
nodes (2 , 1 ) =2 ; nodes (2 , 2 ) =5; 
nodes (3,1) =2; nodes (3, 2) =3; 
nodes ( 4 , 1 ) =3 ; nodes ( 4 , 2 ) =5 ; 
nodes (5,1) =3; nodes (5, 2) =4; 
nodes ( 6, 1 ) =4 ; nodes ( 6, 2 ) =5; 
nodes ( 7 , 1 ) =5 ; nodes ( 7 , 2 ) =1 ; 

bcdof (1) =1; 
bcval ( 1 ) =0 ; 
bcdof (2 ) =2; 
bcval (2) =0 ; 
bcdof (3) =7; 
bcval ( 3 ) =0 ; 
bcdof (4) =8; 
bcval ( 4 ) =0 ; 

ff=zeros (sdof , 1) ; 
kk=zeros (sdof,  sdof)  ; 
index=zeros (nnel^ndof ,  1)  ; 
elf orce=zeros (nnel^ndof ,  1)  ; 
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eldisp=zeros (nnel*ndof , nnel*ndof ) ; 
stress=zeros (nel, 1) ; 

ff (3) =600; 
ff (4) =-600; 
ff (5) =600; 
ff  (6)— 600; 

for  iel=l : nel 

nd (1) =nodes (iel, 1)  ; 
nd (2 ) =nodes (iel, 2)  ; 

xl=gcoord (nd ( 1 ) ,1) ; yl=gcoord (nd ( 1 )  ,  2 ) 
x2=gcoord (nd (2 ) ,1) ; y2=gcoord (nd (2 ) f2) 

leng=sqrt ( (x2-xl) A2+ (y2-yl) A2) ; 

if  (x2-xl ) ==0 ; 

beta=2*atan (1) ; 

else 

beta=atan ( (y2-yl) / (x2-xl) ) ; 

end 

el=elprop  (iel, 1) ; 
area=elprop  (iel , 2 ) ; 

[ index] =FEELDOF (nd, nnel , ndof) ; 

k=FETRUSS2 (el, leng, area, 0, beta, 1) ; 

kk=FEASMBLl (kk, k, index) ; 


end 

[ kk, f f ] =FEAPLYC2 ( kk, f f , bcdof , bcval ) ; 

disp=kk\f f ; 

for  iel=l : nel 

nd ( 1 ) =nodes  ( iel , 1 ) ; 
nd (2 ) =nodes  (iel, 2 ) ; 

xl=gcoord (nd ( 1 ) , 1 ) ; yl=gcoord (nd  ( 1 ) , 2 ) 
x2=gcoord (nd (2 ) , 1 ) ; y2=gcoord (nd (2 )  ,  2 ) 

leng=sqrt ( (x2-xl) A2+ (y2-yl) A2) ; 

if  (x2-xl ) ==0 ; 

beta=2*atan (1) ; 

else 

beta=atan ( (y2-yl) / (x2-xl) ) ; 
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end 


el=elprop  (iel, 1) ; 
area=elprop  (iel, 2 ) ; 
index=FEELDOF (nd, nnel , ndof ) ; 

k=FETRUSS2 (el , leng, area, 0 , beta, 1 ) ; 

for  i=l :  (nnel*ndof) 

eldisp (i) =disp (index (i) ) ; 

end 

elf orce=k* eldisp; 

stress (iel) =sqrt (el force (1) A2+elforce (2) A2) /area 

if ( (x2-xl) *el force  (3) )  <0 ; 

stress (iel) ^-stress  (iel) ; 

end 


end 

num=l : 1 : sdof ; 
displ=[num'  disp] 

numm=l : 1 : nel / 
stresses^ [numm 1  stress] 
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FEMFlb 


nel=7 ; 
nnel=2 ; 
ndof=2 ; 
nnode=5 ; 

sdof =nnode*ndof ; 

gcoord ( 1 , 1 ) =0 ; gcoord ( 1 , 2 ) =0 ; 
gcoord (2 , 1) =0 ; gcoord (2 , 2 ) =12  ; 
gcoord (3, 1) =6 . 3 65; gcoord (3, 2) =  12; 
gcoord (4,1)=6.365; gcoord (4 , 2 ) =0 ; 
gcoord  (5 , 1 ) =3 . 1825 ;  gcoord  (5 , 2 ) =6; 

elprop (1, 1) =10300000; 
elprop ( 1 , 2 ) = . 25 ; 
elprop (2 ,  1)=10300000; 
elprop  (2  ,  2 )  =  . 2  5 ; 
elprop (3,1) =10300000; 
elprop  (3, 2)  =  . 25; 
elprop (4, 1)=10300000; 
elprop  ( 4  ,  2 )  =  . 2  5 ; 
elprop (5, 1) =10300000; 
elprop ( 5 ,  2)  =  .25; 
elprop (6, 1) =10300000; 
elprop ( 6 ,  2)  =  .25; 
elprop (7 , 1)=10300000; 
elprop  ( 7 , 2 )  =  . 2  5 ; 

nodes ( 1 , 1 ) =1 ; nodes (1,2) =2; 
nodes (2,1) =2 ; nodes (2 , 2) =5; 
nodes (3,1) =2; nodes (3, 2) =3; 
nodes (4,1) =3 ; nodes (4,2) =5 ; 
nodes (5,1) =3; nodes (5, 2) =4; 
nodes (6, 1) =4; nodes (6, 2) =5; 
nodes (7,1) =5; nodes (7, 2) =1; 

bcdof ( 1 ) =1 ; 
bcval ( 1 ) =0 ; 
bcdof (2) =2; 
bcval (2 ) =0 ; 
bcdof (3) =7; 
bcval (3 ) =0 ; 
bcdof (4) =8; 
bcval ( 4 ) =0 ; 

ff=zeros (sdof , 1) ; 
kk=zeros (sdof,  sdof)  ; 
index=zeros (nnel*ndof ,  1)  ; 
elf orce=zeros (nnel*ndof ,  1)  ; 
eldisp=zeros (nnel*ndof , nnel*ndof ) 
stress=zeros (nel, 1) ; 


ff (3)=-600; 
ff (4)«-600; 
ff (5) =-600; 
ff (6) =-600; 

for  iel=l:nel 

nd ( 1 ) =nodes ( iel , 1 ) ; 
nd (2 ) =nodes (iel, 2 ) ; 

xl=gcoord (nd ( 1 ) , 1 ) ; yl=gcoord (nd ( 1 ) ,2) 
x2=gcoord (nd ( 2 ) ,  1 ) ; y2=gcoord (nd ( 2 ) ,2) 

leng=sqrt ( (x2-xl) A2+ (y2-yl) A2) ; 

if  (x2-xl ) ==0 ; 

beta=2*atan (1) ; 

else 

beta=atan ( (y2-yl) / (x2-xl) ) ; 

end 

el=elprop  (iel, 1) ; 
area=elprop  ( iel , 2 ) ; 

[ index] =FEELDOF (nd, nnel , ndof) ; 

k=FETRUSS2 (el , leng, area, 0 , beta, 1 ) ; 

kk=FEASMBLl (kk, k, index)  ; 


end 

[kk, f f ] =FEAPLYC2 (kk, f f , bcdof , bcval) ; 

disp=kk\f f ; 

for  iel=l:nel 

nd (1) =nodes  (iel, 1) ; 
nd (2 ) =nodes (iel, 2 ) ; 

xl=gcoord (nd ( 1 ) , 1 ) ; yl=gcoord  (nd  ( 1 ) , 2 ) 
x2=gcoord (nd ( 2 )  ,  1 ) ; y2=gcoord (nd ( 2 )  ,  2 ) 

leng=sqrt ( (x2-xl) A2+ (y2-yl) A2) ; 

if  (x2-xl ) ==0 ; 

beta=2*atan (1) ; 

else 

beta=atan ( (y2-yl) / (x2-xl) ) ; 

end 
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el=elprop  (iel, 1) / 
area=elprop  (iel , 2 ) ; 
index=FEELDOF (nd, nnel , ndof ) ; 

k=FETRUSS2 (el , leng, area, 0 , beta, 1 ) ; 

for  i=l :  (nnel*ndof) 

eldisp (i) =disp (index (i) ) ; 

end 

elf orce=k*eldisp; 

stress (iel) =sqrt (el force (1) A2+elforce (2) A2) /area 

if ( (x2-xl) *elforce (3) )<0; 

stress (iel) ^-stress  (iel) ; 

end 


end 

num=l : 1 : sdof ; 
displ=[num'  disp] 

numm=l : 1 : nel / 
stresses^ [numm 1  stress] 
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FEMFlc 


nel=7 ; 
nnel=2 ; 
ndof=2 ; 
nnode=5 ; 

sdof =nnode*ndof ; 

gcoord ( 1 , 1 ) =0 ; gcoord ( 1 , 2 ) =0 ; 
gcoord (2 , 1) =0 ; gcoord (2 , 2 ) =12  ; 
gcoord (3, 1) =6 . 3 65; gcoord (3, 2) =  12; 
gcoord (4,1)=6.365; gcoord (4 , 2 ) =0 ; 
gcoord  (5 , 1 ) =3 . 1825 ;  gcoord  (5 , 2 ) =6; 

elprop (1, 1) =10300000; 
elprop ( 1 , 2 ) = . 25 ; 
elprop (2 ,  1)=10300000; 
elprop  (2  ,  2 )  =  . 2  5 ; 
elprop (3,1) =10300000; 
elprop  (3, 2)  =  . 25; 
elprop (4, 1)=10300000; 
elprop  ( 4  ,  2 )  =  . 2  5 ; 
elprop (5, 1) =10300000; 
elprop ( 5 ,  2)  =  .25; 
elprop (6, 1) =10300000; 
elprop ( 6 ,  2)  =  .25; 
elprop (7 , 1)=10300000; 
elprop  ( 7 , 2 )  =  . 2  5 ; 

nodes ( 1 , 1 ) =1 ; nodes (1,2) =2; 
nodes (2,1) =2 ; nodes (2 , 2) =5; 
nodes (3,1) =2; nodes (3, 2) =3; 
nodes (4,1) =3 ; nodes (4,2) =5 ; 
nodes (5,1) =3; nodes (5, 2) =4; 
nodes (6, 1) =4; nodes (6, 2) =5; 
nodes (7,1) =5; nodes (7, 2) =1; 

bcdof ( 1 ) =1 ; 
bcval ( 1 ) =0 ; 
bcdof (2) =2; 
bcval (2 ) =0 ; 
bcdof (3) =7; 
bcval (3 ) =0 ; 
bcdof (4) =8; 
bcval ( 4 ) =0 ; 

ff=zeros (sdof , 1) ; 
kk=zeros (sdof,  sdof)  ; 
index=zeros (nnel*ndof ,  1)  ; 
elf orce=zeros (nnel*ndof ,  1)  ; 
eldisp=zeros (nnel*ndof , nnel*ndof ) 
stress=zeros (nel, 1) ; 


ff (4)=600; 
ff (6)=600; 


for  iel=l : nel 

nd (1) =nodes  (iel, 1) ; 
nd (2 ) =nodes (iel, 2 ) ; 

xl=gcoord (nd ( 1 ) ,  1 ) ; yl=gcoord (nd ( 1 ) ,  2 ) 
x2=gcoord (nd (2 ) ,  1 ) ; y2=gcoord (nd (2 ) ,2) 

leng=sqrt ( (x2-xl) A2+ (y2-yl) A2) ; 

if  (x2-xl ) ==0 ; 

beta=2*atan (1) ; 

else 

beta=atan ( (y2-yl) / (x2-xl) ) ; 

end 

el=elprop  (iel, 1) ; 
area=elprop  (ielf  2) ; 

[ index] =FEELDOF (ndf nnel A ndof) ; 

k=FETRUSS2 (el,leng, area, 0 , beta, 1 ) ; 

kk=FEASMBLl (kk, k, index)  ; 


end 

[ kk, f f ] =FEAPLYC2 ( kk, f f , bcdof , bcval ) ; 

disp=kk\f f ; 

for  iel=l : nel 

nd (1) =nodes  (iel, 1) ; 
nd (2 ) =nodes  (iel, 2 ) ; 

xl=gcoord (nd ( 1 ) , 1 ) ; yl=gcoord (nd  ( 1 ) , 2 ) 
x2=gcoord (nd (2 ) , 1 ) ; y2=gcoord (nd (2 )  ,  2 ) 

leng=sqrt ( (x2-xl) A2+ (y2-yl) A2) ; 

if  (x2-xl ) ==0 ; 

beta=2^atan (1) ; 

else 

beta=atan ( (y2-yl) / (x2-xl) ) / 

end 

el=elprop  (iel, 1) ; 
area=elprop  (iel, 2 ) ; 
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index=FEELDOF (nd, nnel , ndof ) ; 


k=FETRUSS2 (el, leng, area, 0 , beta,  1 ) ; 

for  i=l :  (nnel*ndof) 

eldisp (i) =disp (index (i)  )  ; 

end 

elf orce=k*eldisp; 

stress (iel) =sqrt (elforce (1) A2+elforce (2) A2) /area 

if ( (x2-xl ) *el force  ( 3 ) )  <0 ; 

stress (iel) ^-stress  (iel) ; 

end 


end 

num=l : 1 : sdof  ; 
displ=[num'  disp] 

numm=l : 1 : nel / 
stresses^ [numm 1  stress] 


139 


FEMFld 


nel=7 ; 
nnel=2 ; 
ndof=2 ; 
nnode=5 ; 

sdof =nnode*ndof ; 

gcoord ( 1 , 1 ) =0 ; gcoord ( 1 , 2 ) =0 ; 
gcoord (2 , 1) =0 ; gcoord (2 , 2 ) =12  ; 
gcoord (3, 1) =6 . 3 65; gcoord (3, 2) =  12; 
gcoord (4,1)=6.365; gcoord (4 , 2 ) =0 ; 
gcoord  (5 , 1 ) =3 . 1825 ;  gcoord  (5 , 2 ) =6; 

elprop (1, 1) =10300000; 
elprop ( 1 , 2 ) = . 25 ; 
elprop (2 ,  1)=10300000; 
elprop  (2  ,  2 )  =  . 2  5 ; 
elprop (3,1) =10300000; 
elprop  (3, 2)  =  . 25; 
elprop (4, 1)=10300000; 
elprop  ( 4  ,  2 )  =  . 2  5 ; 
elprop (5, 1) =10300000; 
elprop ( 5 ,  2)  =  .25; 
elprop (6, 1) =10300000; 
elprop ( 6 ,  2)  =  .25; 
elprop (7 , 1)=10300000; 
elprop  ( 7 , 2 )  =  . 2  5 ; 

nodes ( 1 , 1 ) =1 ; nodes (1,2) =2; 
nodes (2,1) =2 ; nodes (2 , 2) =5; 
nodes (3,1) =2; nodes (3, 2) =3; 
nodes (4,1) =3 ; nodes (4,2) =5 ; 
nodes (5,1) =3; nodes (5, 2) =4; 
nodes (6, 1) =4; nodes (6, 2) =5; 
nodes (7,1) =5; nodes (7, 2) =1; 

bcdof ( 1 ) =1 ; 
bcval ( 1 ) =0 ; 
bcdof (2) =2; 
bcval (2 ) =0 ; 
bcdof (3) =7; 
bcval (3 ) =0 ; 
bcdof (4) =8; 
bcval ( 4 ) =0 ; 

ff=zeros (sdof , 1) ; 
kk=zeros (sdof,  sdof)  ; 
index=zeros (nnel*ndof ,  1)  ; 
elf orce=zeros (nnel*ndof ,  1)  ; 
eldisp=zeros (nnel*ndof , nnel*ndof ) 
stress=zeros (nel, 1) ; 


f f (3) =600; 
ff (4)=600; 
ff (5) =600; 
ff (6) =600; 

for  iel=l:nel 

nd ( 1 ) =nodes ( iel , 1 ) ; 
nd (2 ) =nodes  (iel, 2 ) ; 

xl=gcoord (nd ( 1 ) , 1 ) ; yl=gcoord  (nd  ( 1 ) ,  2 ) 
x2=gcoord (nd ( 2 ) , 1 ) ; y2=gcoord (nd ( 2 ) , 2 ) 

leng=sqrt ( (x2-xl) A2+ (y2-yl) A2) ; 

if  (x2-xl ) ==0 ; 

beta=2*atan (1) ; 

else 

beta=atan ( (y2-yl) / (x2-xl) ) ; 

end 

el=elprop  (iel, 1) ; 
area=elprop  ( iel ,  2 ) ; 

[ index] =FEELDOF (nd, nnel , ndof) ; 

k=FETRUSS2 (el, lengf area, 0 , beta, 1 ) ; 

kk=FEASMBLl (kk, k, index)  ; 


end 

[kk, f f ] =FEAPLYC2 (kk, f f , bcdof , bcval) ; 

disp=kk\f f ; 

for  iel=l:nel 

nd (1) =nodes  (iel, 1) ; 
nd (2 ) =nodes (iel, 2 ) ; 

xl=gcoord (nd ( 1 ) , 1 ) ; yl=gcoord  (nd  ( 1 ) , 2 ) 
x2=gcoord (nd (2)  ,  1 ) ; y2=gcoord (nd ( 2 )  ,  2 ) 

leng=sqrt ( (x2-xl) A2+ (y2-yl) A2) ; 

if  (x2-xl ) ==0 ; 

beta=2^atan (1) ; 

else 

beta=atan ( (y2-yl) / (x2-xl) ) ; 

end 
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el=elprop  (iel, 1) / 
area=elprop  (iel , 2 ) ; 
index=FEELDOF (nd, nnel , ndof ) ; 

k=FETRUSS2 (el , leng, area, 0 , beta, 1 ) ; 

for  i=l :  (nnel*ndof) 

eldisp (i) =disp (index (i) ) ; 

end 

elf orce=k*eldisp; 

stress (iel) =sqrt (el force (1) A2+elforce (2) A2) /area 

if ( (x2-xl) *elforce (3) )<0; 

stress (iel) ^-stress  (iel) ; 

end 


end 

num=l : 1 : sdof ; 
displ=[num'  disp] 

numm=l : 1 : nel / 
stresses^ [numm 1  stress] 
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FEMFle 


nel=7 ; 
nnel=2 ; 
ndof=2 ; 
nnode=5 ; 

sdof =nnode*ndof ; 

gcoord ( 1 , 1 ) =0 ; gcoord ( 1 , 2 ) =0 ; 
gcoord (2 , 1) =0 ; gcoord (2 , 2 ) =12  ; 
gcoord (3, 1) =6 . 3 65; gcoord (3, 2) =  12; 
gcoord (4,1)=6.365; gcoord (4 , 2 ) =0 ; 
gcoord  (5 , 1 ) =3 . 1825 ;  gcoord  (5 , 2 ) =6; 

elprop (1, 1) =10300000; 
elprop ( 1 , 2 ) = . 25 ; 
elprop (2 ,  1)=10300000; 
elprop  (2  ,  2 )  =  . 2  5 ; 
elprop (3,1) =10300000; 
elprop  (3, 2)  =  . 25; 
elprop (4, 1)=10300000; 
elprop  ( 4  ,  2 )  =  . 2  5 ; 
elprop (5, 1) =10300000; 
elprop ( 5 ,  2)  =  .25; 
elprop (6, 1) =10300000; 
elprop ( 6 ,  2)  =  .25; 
elprop (7 , 1)=10300000; 
elprop  ( 7 , 2 )  =  . 2  5 ; 

nodes ( 1 , 1 ) =1 ; nodes (1,2) =2; 
nodes (2,1) =2 ; nodes (2 , 2) =5; 
nodes (3,1) =2; nodes (3, 2) =3; 
nodes (4,1) =3 ; nodes (4,2) =5 ; 
nodes (5,1) =3; nodes (5, 2) =4; 
nodes (6, 1) =4; nodes (6, 2) =5; 
nodes (7,1) =5; nodes (7, 2) =1; 

bcdof ( 1 ) =1 ; 
bcval ( 1 ) =0 ; 
bcdof (2) =2; 
bcval (2 ) =0 ; 
bcdof (3) =7; 
bcval (3 ) =0 ; 
bcdof (4) =8; 
bcval ( 4 ) =0 ; 

ff=zeros (sdof , 1) ; 
kk=zeros (sdof,  sdof)  ; 
index=zeros (nnel*ndof ,  1)  ; 
elf orce=zeros (nnel*ndof ,  1)  ; 
eldisp=zeros (nnel*ndof , nnel*ndof ) 
stress=zeros (nel, 1) ; 


ff (3)=-600; 
ff (4)=600; 
ff (5) =-600; 
ff (6) =600; 

for  iel=l:nel 

nd ( 1 ) =nodes ( iel , 1 ) ; 
nd (2 ) =nodes (iel, 2 ) ; 

xl=gcoord (nd ( 1 ) , 1 ) ; yl=gcoord (nd ( 1 ) ,2) 
x2=gcoord (nd ( 2 ) ,  1 ) ; y2=gcoord (nd ( 2 ) ,2) 

leng=sqrt ( (x2-xl) A2+ (y2-yl) A2) ; 

if  (x2-xl ) ==0 ; 

beta=2*atan (1) ; 

else 

beta=atan ( (y2-yl) / (x2-xl) ) ; 

end 

el=elprop  (iel, 1) ; 
area=elprop  ( iel , 2 ) ; 

[ index] =FEELDOF (nd, nnel , ndof) ; 

k=FETRUSS2 (el , leng, area, 0 , beta, 1 ) ; 

kk=FEASMBLl (kk, k, index)  ; 


end 

[kk, f f ] =FEAPLYC2 (kk, f f , bcdof , bcval) ; 

disp=kk\f f ; 

for  iel=l:nel 

nd (1) =nodes  (iel, 1) ; 
nd (2 ) =nodes (iel, 2 ) ; 

xl=gcoord (nd ( 1 ) , 1 ) ; yl=gcoord  (nd  ( 1 ) , 2 ) 
x2=gcoord (nd ( 2 )  ,  1 ) ; y2=gcoord (nd ( 2 )  ,  2 ) 

leng=sqrt ( (x2-xl) A2+ (y2-yl) A2) ; 

if  (x2-xl ) ==0 ; 

beta=2*atan (1) ; 

else 

beta=atan ( (y2-yl) / (x2-xl) ) ; 

end 
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el=elprop  (iel, 1) / 
area=elprop  (iel , 2 ) ; 
index=FEELDOF (nd, nnel , ndof ) ; 

k=FETRUSS2 (el , leng, area, 0 , beta, 1 ) ; 

for  i=l :  (nnel*ndof) 

eldisp (i) =disp (index (i) ) ; 

end 

elf orce=k*eldisp; 

stress (iel) =sqrt (el force (1) A2+elforce (2) A2) /area 

if ( (x2-xl) *elforce (3) )<0; 

stress (iel) ^-stress  (iel) ; 

end 


end 

num=l : 1 : sdof ; 
displ=[num'  disp] 

numm=l : 1 : nel / 
stresses^ [numm 1  stress] 
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FEMFlf 


nel=7 ; 
nnel=2 ; 
ndof=2 ; 
nnode=5 ; 

sdof =nnode*ndof ; 

gcoord ( 1 , 1 ) =0 ; gcoord ( 1 , 2 ) =0 ; 
gcoord (2 , 1) =0 ; gcoord (2 , 2 ) =12  ; 
gcoord (3, 1) =6 . 3 65; gcoord (3, 2) =  12; 
gcoord (4,1)=6.365; gcoord (4 , 2 ) =0 ; 
gcoord  (5 , 1 ) =3 . 1825 ;  gcoord  (5 , 2 ) =6; 

elprop (1, 1) =10300000; 
elprop ( 1 , 2 ) = . 25 ; 
elprop (2 ,  1)=10300000; 
elprop  (2  ,  2 )  =  . 2  5 ; 
elprop (3,1) =10300000; 
elprop  (3, 2)  =  . 25; 
elprop (4, 1)=10300000; 
elprop  ( 4  ,  2 )  =  . 2  5 ; 
elprop (5, 1) =10300000; 
elprop ( 5 ,  2)  =  .25; 
elprop (6, 1) =10300000; 
elprop ( 6 ,  2)  =  .25; 
elprop (7 , 1)=10300000; 
elprop  ( 7 , 2 )  =  . 2  5 ; 

nodes ( 1 , 1 ) =1 ; nodes (1,2) =2; 
nodes (2,1) =2 ; nodes (2 , 2) =5; 
nodes (3,1) =2; nodes (3, 2) =3; 
nodes (4,1) =3 ; nodes (4,2) =5 ; 
nodes (5,1) =4; nodes (5, 2) =5; 
nodes (6,1) =4 ; nodes (6, 2) =1; 
nodes (7,1) =5; nodes (7, 2) =1; 

bcdof ( 1 ) =5 ; 
bcval ( 1 ) =0 ; 
bcdof (2) =6; 
bcval (2 ) =0 ; 
bcdof (3) =7; 
bcval (3 ) =0 ; 
bcdof (4) =8; 
bcval ( 4 ) =0 ; 

ff=zeros (sdof , 1) ; 
kk=zeros (sdof, sdof)  ; 
index=zeros (nnel*ndof ,  1)  ; 
elf orce=zeros (nnel*ndof ,  1)  ; 
eldisp=zeros (nnel*ndof , nnel*ndof ) 
stress=zeros (nel, 1)  ; 


ff (2)=-600; 
ff (4)=-600; 


for  iel=l : nel 

nd (1) =nodes  (iel, 1) ; 
nd (2 ) =nodes (iel, 2 ) ; 

xl=gcoord (nd ( 1 ) ,  1 ) ; yl=gcoord (nd ( 1 ) ,  2 ) 
x2=gcoord (nd (2 ) ,  1 ) ; y2=gcoord (nd (2 ) ,2) 

leng=sqrt ( (x2-xl) A2+ (y2-yl) A2) ; 

if  (x2-xl ) ==0 ; 

beta=2*atan (1) ; 

else 

beta=atan ( (y2-yl) / (x2-xl) ) ; 

end 

el=elprop  (iel, 1) ; 
area=elprop  (ielf  2) ; 

[ index] =FEELDOF (ndf nnel A ndof) ; 

k=FETRUSS2 (el,leng, area, 0 , beta, 1 ) ; 

kk=FEASMBLl (kk, k, index)  ; 


end 

[ kk, f f ] =FEAPLYC2 ( kk, f f , bcdof , bcval ) ; 

disp=kk\f f ; 

for  iel=l : nel 

nd (1) =nodes  (iel, 1) ; 
nd (2 ) =nodes  (iel, 2 ) ; 

xl=gcoord (nd ( 1 ) , 1 ) ; yl=gcoord (nd  ( 1 ) , 2 ) 
x2=gcoord (nd (2 ) , 1 ) ; y2=gcoord (nd (2 )  ,  2 ) 

leng=sqrt ( (x2-xl) A2+ (y2-yl) A2) ; 

if  (x2-xl ) ==0 ; 

beta=2^atan (1) ; 

else 

beta=atan ( (y2-yl) / (x2-xl) ) / 

end 

el=elprop  (iel, 1) ; 
area=elprop  (iel, 2 ) ; 
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index=FEELDOF (nd, nnel , ndof ) ; 


k=FETRUSS2 (el, leng, area, 0 , beta,  1 ) ; 

for  i=l :  (nnel*ndof) 

eldisp (i) =disp (index (i)  )  ; 

end 

elf orce=k*eldisp; 

stress (iel) =sqrt (elforce (1) A2+elforce (2) A2) /area 

if ( (x2-xl ) *el force  ( 3 ) )  <0 ; 

stress (iel) ^-stress  (iel) ; 

end 


end 

num=l : 1 : sdof  ; 
displ=[num'  disp] 

numm=l : 1 : nel / 
stresses^ [numm 1  stress] 
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FEMFlg 


nel=7 ; 
nnel=2 ; 
ndof=2 ; 
nnode=5 ; 

sdof =nnode*ndof ; 

gcoord ( 1 , 1 ) =0 ; gcoord ( 1 , 2 ) =0 ; 
gcoord (2 , 1) =0 ; gcoord (2 , 2 ) =12  ; 
gcoord (3, 1) =6 . 3 65; gcoord (3, 2) =  12; 
gcoord (4,1)=6.365; gcoord (4 , 2 ) =0 ; 
gcoord  (5 , 1 ) =3 . 1825 ;  gcoord  (5 , 2 ) =6; 

elprop (1, 1) =10300000; 
elprop ( 1 , 2 ) = . 25 ; 
elprop (2 ,  1)=10300000; 
elprop  (2  ,  2 )  =  . 2  5 ; 
elprop (3,1) =10300000; 
elprop  (3, 2)  =  . 25; 
elprop (4, 1)=10300000; 
elprop  ( 4  ,  2 )  =  . 2  5 ; 
elprop (5, 1) =10300000; 
elprop ( 5 ,  2)  =  .25; 
elprop (6, 1) =10300000; 
elprop ( 6 ,  2)  =  .25; 
elprop (7 , 1)=10300000; 
elprop  ( 7 , 2 )  =  . 2  5 ; 

nodes ( 1 , 1 ) =1 ; nodes (1,2) =2; 
nodes (2,1) =2 ; nodes (2 , 2) =5; 
nodes (3,1) =2; nodes (3, 2) =3; 
nodes (4,1) =3 ; nodes (4,2) =5 ; 
nodes (5,1) =4; nodes (5, 2) =5; 
nodes (6,1) =4 ; nodes (6, 2) =1; 
nodes (7,1) =5; nodes (7, 2) =1; 

bcdof ( 1 ) =5 ; 
bcval ( 1 ) =0 ; 
bcdof (2) =6; 
bcval (2 ) =0 ; 
bcdof (3) =7; 
bcval (3 ) =0 ; 
bcdof (4) =8; 
bcval ( 4 ) =0 ; 

ff=zeros (sdof , 1) ; 
kk=zeros (sdof, sdof)  ; 
index=zeros (nnel*ndof ,  1)  ; 
elf orce=zeros (nnel*ndof ,  1)  ; 
eldisp=zeros (nnel*ndof , nnel*ndof )  ; 
stress=zeros (nel, 1)  ; 
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ff (1) =600; 
ff (2)«-600; 
ff (3) =600; 
ff (4) =-600; 

for  iel=l:nel 

nd ( 1 ) =nodes ( iel , 1 ) ; 
nd (2 ) =nodes  (iel, 2 ) ; 

xl=gcoord (nd ( 1 ) , 1 ) ; yl=gcoord  (nd  ( 1 ) ,  2 ) 
x2=gcoord (nd ( 2 ) , 1 ) ; y2=gcoord (nd ( 2 )  ,  2 ) 

leng=sqrt ( (x2-xl) A2+ (y2-yl) A2) ; 

if  (x2-xl ) ==0 ; 

beta=2*atan (1) ; 

else 

beta=atan ( (y2-yl) / (x2-xl) ) ; 

end 

el=elprop  (iel, 1) ; 
area=elprop  ( iel ,  2 ) ; 

[ index] =FEELDOF (nd, nnel , ndof) ; 

k=FETRUSS2 (el, lengf area, 0 , beta, 1 ) ; 

kk=FEASMBLl (kk, k, index)  ; 


end 

[kk, f f ] =FEAPLYC2 (kk, f f , bcdof , bcval) ; 

disp=kk\f f ; 

for  iel=l:nel 

nd (1) =nodes  (iel, 1) ; 
nd (2 ) =nodes (iel, 2 ) ; 

xl=gcoord (nd ( 1 ) , 1 ) ; yl=gcoord  (nd  ( 1 ) , 2 ) 
x2=gcoord (nd (2)  ,  1 ) ; y2=gcoord (nd ( 2 )  ,  2 ) 

leng=sqrt ( (x2-xl) A2+ (y2-yl) A2) ; 

if  (x2-xl ) ==0 ; 

beta=2^atan (1) ; 

else 

beta=atan ( (y2-yl) / (x2-xl) ) ; 

end 
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el=elprop  (iel, 1) / 
area=elprop  (iel , 2 ) ; 
index=FEELDOF (nd, nnel , ndof ) ; 

k=FETRUSS2 (el , leng, area, 0 , beta, 1 ) ; 

for  i=l :  (nnel*ndof) 

eldisp (i) =disp (index (i) ) ; 

end 

elf orce=k*eldisp; 

stress (iel) =sqrt (el force (1) A2+elforce (2) A2) /area 

if ( (x2-xl) *elforce (3) )<0; 

stress (iel) ^-stress  (iel) ; 

end 


end 

num=l : 1 : sdof ; 
displ=[num'  disp] 

numm=l : 1 : nel / 
stresses^ [numm 1  stress] 
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FEMFlh 


nel=7 ; 
nnel=2 ; 
ndof=2 ; 
nnode=5 ; 

sdof =nnode*ndof ; 

gcoord ( 1 , 1 ) =0 ; gcoord ( 1 , 2 ) =0 ; 
gcoord (2 , 1) =0 ; gcoord (2 , 2 ) =12  ; 
gcoord (3, 1) =6 . 3 65; gcoord (3, 2) =  12; 
gcoord (4,1)=6.365; gcoord (4 , 2 ) =0 ; 
gcoord  (5 , 1 ) =3 . 1825 ;  gcoord  (5 , 2 ) =6; 

elprop (1, 1) =10300000; 
elprop ( 1 , 2 ) = . 25 ; 
elprop (2 ,  1)=10300000; 
elprop  (2  ,  2 )  =  . 2  5 ; 
elprop (3,1) =10300000; 
elprop  (3, 2)  =  . 25; 
elprop (4, 1)=10300000; 
elprop  ( 4  ,  2 )  =  . 2  5 ; 
elprop (5, 1) =10300000; 
elprop ( 5 ,  2)  =  .25; 
elprop (6, 1) =10300000; 
elprop ( 6 ,  2)  =  .25; 
elprop (7 , 1)=10300000; 
elprop  ( 7 , 2 )  =  . 2  5 ; 

nodes ( 1 , 1 ) =1 ; nodes (1,2) =2; 
nodes (2,1) =2 ; nodes (2 , 2) =5; 
nodes (3,1) =2; nodes (3, 2) =3; 
nodes (4,1) =3 ; nodes (4,2) =5 ; 
nodes (5,1) =4; nodes (5, 2) =5; 
nodes (6,1) =4 ; nodes (6, 2) =1; 
nodes (7,1) =5; nodes (7, 2) =1; 

bcdof ( 1 ) =5 ; 
bcval ( 1 ) =0 ; 
bcdof (2) =6; 
bcval (2 ) =0 ; 
bcdof (3) =7; 
bcval (3 ) =0 ; 
bcdof (4) =8; 
bcval ( 4 ) =0 ; 

ff=zeros (sdof , 1) ; 
kk=zeros (sdof, sdof)  ; 
index=zeros (nnel*ndof ,  1)  ; 
elf orce=zeros (nnel*ndof ,  1)  ; 
eldisp=zeros (nnel*ndof , nnel*ndof ) 
stress=zeros (nel, 1)  ; 


ff (l)=-600; 
ff  (2)—  600; 
ff (3) =-600; 
ff (4) =-600; 

for  iel=l:nel 

nd ( 1 ) =nodes ( iel , 1 ) ; 
nd (2 ) =nodes (iel, 2 ) ; 

xl=gcoord (nd ( 1 ) , 1 ) ; yl=gcoord  (nd  ( 1 ) ,  2 ) 
x2=gcoord (nd ( 2 )  ,  1 ) ; y2=gcoord (nd ( 2 )  ,  2 ) 

leng=sqrt ( (x2-xl) A2+ (y2-yl) A2) ; 

if  (x2-xl ) ==0 ; 

beta=2*atan (1) ; 

else 

beta=atan ( (y2-yl) / (x2-xl) ) ; 

end 

el=elprop  (iel, 1) ; 
area=elprop  ( iel ,  2 ) ; 

[ index] =FEELDOF (nd, nnel , ndof) ; 

k=FETRUSS2 (el, lengf area, 0 , beta, 1 ) ; 

kk=FEASMBLl (kk, k, index)  ; 


end 

[kk, f f ] =FEAPLYC2 (kk, f f , bcdof , bcval) ; 

disp=kk\f f ; 

for  iel=l:nel 

nd (1) =nodes  (iel, 1) ; 
nd (2 ) =nodes (iel, 2 ) ; 

xl=gcoord (nd ( 1 ) , 1 ) ; yl=gcoord  (nd  ( 1 ) , 2 ) 
x2=gcoord (nd (2)  ,  1 ) ; y2=gcoord (nd ( 2 )  ,  2 ) 

leng=sqrt ( (x2-xl) A2+ (y2-yl) A2) ; 

if  (x2-xl ) ==0 ; 

beta=2^atan (1) ; 

else 

beta=atan ( (y2-yl) / (x2-xl) ) ; 

end 
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el=elprop  (iel, 1) / 
area=elprop  (iel , 2 ) ; 
index=FEELDOF (nd, nnel , ndof ) ; 

k=FETRUSS2 (el , leng, area, 0 , beta, 1 ) ; 

for  i=l :  (nnel*ndof) 

eldisp (i) =disp (index (i) ) ; 

end 

elf orce=k*eldisp; 

stress (iel) =sqrt (el force (1) A2+elforce (2) A2) /area 

if ( (x2-xl) *elforce (3) )<0; 

stress (iel) ^-stress  (iel) ; 

end 


end 

num=l : 1 : sdof ; 
displ=[num'  disp] 

numm=l : 1 : nel / 
stresses^ [numm 1  stress] 
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FEMFli 


nel=7 ; 
nnel=2 ; 
ndof=2 ; 
nnode=5 ; 

sdof =nnode*ndof ; 

gcoord ( 1 , 1 ) =0 ; gcoord ( 1 , 2 ) =0 ; 
gcoord (2 , 1) =0 ; gcoord (2 , 2 ) =12  ; 
gcoord (3, 1) =6 . 3 65; gcoord (3, 2) =  12; 
gcoord (4,1)=6.365; gcoord (4 , 2 ) =0 ; 
gcoord  (5 , 1 ) =3 . 1825 ;  gcoord  (5 , 2 ) =6; 

elprop (1, 1) =10300000; 
elprop ( 1 , 2 ) = . 25 ; 
elprop (2 ,  1)=10300000; 
elprop  (2  ,  2 )  =  . 2  5 ; 
elprop (3,1) =10300000; 
elprop  (3, 2)  =  . 25; 
elprop (4, 1)=10300000; 
elprop  ( 4  ,  2 )  =  . 2  5 ; 
elprop (5, 1) =10300000; 
elprop ( 5 ,  2)  =  .25; 
elprop (6, 1) =10300000; 
elprop ( 6 ,  2)  =  .25; 
elprop (7 , 1)=10300000; 
elprop  ( 7 , 2 )  =  . 2  5 ; 

nodes ( 1 , 1 ) =1 ; nodes (1,2) =2; 
nodes (2,1) =2 ; nodes (2 , 2) =5; 
nodes (3,1) =2; nodes (3, 2) =3; 
nodes (4,1) =3 ; nodes (4,2) =5 ; 
nodes (5,1) =4; nodes (5, 2) =5; 
nodes (6,1) =4 ; nodes (6, 2) =1; 
nodes (7,1) =5; nodes (7, 2) =1; 

bcdof ( 1 ) =5 ; 
bcval ( 1 ) =0 ; 
bcdof (2) =6; 
bcval (2 ) =0 ; 
bcdof (3) =7; 
bcval (3 ) =0 ; 
bcdof (4) =8; 
bcval ( 4 ) =0 ; 

ff=zeros (sdof , 1) ; 
kk=zeros (sdof, sdof)  ; 
index=zeros (nnel*ndof ,  1)  ; 
elf orce=zeros (nnel*ndof ,  1)  ; 
eldisp=zeros (nnel*ndof , nnel*ndof ) 
stress=zeros (nel, 1)  ; 


ff (2)=600; 
ff (4)=600; 


for  iel=l : nel 

nd (1) =nodes  (iel, 1) ; 
nd (2 ) =nodes (iel, 2 ) ; 

xl=gcoord (nd ( 1 ) ,  1 ) ; yl=gcoord (nd ( 1 ) ,  2 ) 
x2=gcoord (nd (2 ) ,  1 ) ; y2=gcoord (nd (2 ) ,2) 

leng=sqrt ( (x2-xl) A2+ (y2-yl) A2) ; 

if  (x2-xl ) ==0 ; 

beta=2*atan (1) ; 

else 

beta=atan ( (y2-yl) / (x2-xl) ) ; 

end 

el=elprop  (iel, 1) ; 
area=elprop  (ielf  2) ; 

[ index] =FEELDOF (ndf nnel A ndof) ; 

k=FETRUSS2 (el,leng, area, 0 , beta, 1 ) ; 

kk=FEASMBLl (kk, k, index)  ; 


end 

[ kk, f f ] =FEAPLYC2 ( kk, f f , bcdof , bcval ) ; 

disp=kk\f f ; 

for  iel=l : nel 

nd (1) =nodes  (iel, 1) ; 
nd (2 ) =nodes  (iel, 2 ) ; 

xl=gcoord (nd ( 1 ) , 1 ) ; yl=gcoord (nd  ( 1 ) , 2 ) 
x2=gcoord (nd (2 ) , 1 ) ; y2=gcoord (nd (2 )  ,  2 ) 

leng=sqrt ( (x2-xl) A2+ (y2-yl) A2) ; 

if  (x2-xl ) ==0 ; 

beta=2^atan (1) ; 

else 

beta=atan ( (y2-yl) / (x2-xl) ) / 

end 

el=elprop  (iel, 1) ; 
area=elprop  (iel, 2 ) ; 
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index=FEELDOF (nd, nnel , ndof ) ; 


k=FETRUSS2 (el, leng, area, 0 , beta,  1 ) ; 

for  i=l :  (nnel*ndof) 

eldisp (i) =disp (index (i)  )  ; 

end 

elf orce=k*eldisp; 

stress (iel) =sqrt (elforce (1) A2+elforce (2) A2) /area 

if ( (x2-xl ) *el force  ( 3 ) )  <0 ; 

stress (iel) ^-stress  (iel) ; 

end 


end 

num=l : 1 : sdof  ; 
displ=[num'  disp] 

numm=l : 1 : nel / 
stresses^ [numm 1  stress] 
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FEMFlj 


nel=7 ; 
nnel=2 ; 
ndof=2 ; 
nnode=5 ; 

sdof =nnode*ndof ; 

gcoord ( 1 , 1 ) =0 ; gcoord ( 1 , 2 ) =0 ; 
gcoord (2 , 1) =0 ; gcoord (2 , 2 ) =12  ; 
gcoord (3, 1) =6 . 3 65; gcoord (3, 2) =  12; 
gcoord (4,1)=6.365; gcoord (4 , 2 ) =0 ; 
gcoord  (5 , 1 ) =3 . 1825 ;  gcoord  (5 , 2 ) =6; 

elprop (1, 1) =10300000; 
elprop ( 1 , 2 ) = . 25 ; 
elprop (2 ,  1)=10300000; 
elprop  (2  ,  2 )  =  . 2  5 ; 
elprop (3,1) =10300000; 
elprop  (3, 2)  =  . 25; 
elprop (4, 1)=10300000; 
elprop  ( 4  ,  2 )  =  . 2  5 ; 
elprop (5, 1) =10300000; 
elprop ( 5 ,  2)  =  .25; 
elprop (6, 1) =10300000; 
elprop ( 6 ,  2)  =  .25; 
elprop (7 , 1)=10300000; 
elprop  ( 7 , 2 )  =  . 2  5 ; 

nodes ( 1 , 1 ) =1 ; nodes (1,2) =2; 
nodes (2,1) =2 ; nodes (2 , 2) =5; 
nodes (3,1) =2; nodes (3, 2) =3; 
nodes (4,1) =3 ; nodes (4,2) =5 ; 
nodes (5,1) =4; nodes (5, 2) =5; 
nodes (6,1) =4 ; nodes (6, 2) =1; 
nodes (7,1) =5; nodes (7, 2) =1; 

bcdof ( 1 ) =5 ; 
bcval ( 1 ) =0 ; 
bcdof (2) =6; 
bcval (2 ) =0 ; 
bcdof (3) =7; 
bcval (3 ) =0 ; 
bcdof (4) =8; 
bcval ( 4 ) =0 ; 

ff=zeros (sdof , 1) ; 
kk=zeros (sdof, sdof)  ; 
index=zeros (nnel*ndof ,  1)  ; 
elf orce=zeros (nnel*ndof ,  1)  ; 
eldisp=zeros (nnel*ndof , nnel*ndof )  ; 
stress=zeros (nel, 1)  ; 
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f f (1) =600; 
ff (2)=600; 
ff (3) =600; 
ff (4) =600; 

for  iel=l:nel 

nd ( 1 ) =nodes ( iel , 1 ) ; 
nd (2 ) =nodes  (iel, 2 ) ; 

xl=gcoord (nd ( 1 ) , 1 ) ; yl=gcoord  (nd  ( 1 ) ,  2 ) 
x2=gcoord (nd ( 2 ) , 1 ) ; y2=gcoord (nd ( 2 ) , 2 ) 

leng=sqrt ( (x2-xl) A2+ (y2-yl) A2) ; 

if  (x2-xl ) ==0 ; 

beta=2*atan (1) ; 

else 

beta=atan ( (y2-yl) / (x2-xl) ) ; 

end 

el=elprop  (iel, 1) ; 
area=elprop  ( iel ,  2 ) ; 

[ index] =FEELDOF (nd, nnel , ndof) ; 

k=FETRUSS2 (el, lengf area, 0 , beta, 1 ) ; 

kk=FEASMBLl (kk, k, index)  ; 


end 

[kk, f f ] =FEAPLYC2 (kk, f f , bcdof , bcval) ; 

disp=kk\f f ; 

for  iel=l:nel 

nd (1) =nodes  (iel, 1) ; 
nd (2 ) =nodes (iel, 2 ) ; 

xl=gcoord (nd ( 1 ) , 1 ) ; yl=gcoord  (nd  ( 1 ) , 2 ) 
x2=gcoord (nd (2)  ,  1 ) ; y2=gcoord (nd ( 2 )  ,  2 ) 

leng=sqrt ( (x2-xl) A2+ (y2-yl) A2) ; 

if  (x2-xl ) ==0 ; 

beta=2^atan (1) ; 

else 

beta=atan ( (y2-yl) / (x2-xl) ) ; 

end 
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el=elprop  (iel, 1) / 
area=elprop  (iel , 2 ) ; 
index=FEELDOF (nd, nnel , ndof ) ; 

k=FETRUSS2 (el , leng, area, 0 , beta, 1 ) ; 

for  i=l :  (nnel*ndof) 

eldisp (i) =disp (index (i) ) ; 

end 

elf orce=k*eldisp; 

stress (iel) =sqrt (el force (1) A2+elforce (2) A2) /area 

if ( (x2-xl) *elforce (3) )<0; 

stress (iel) ^-stress  (iel) ; 

end 


end 

num=l : 1 : sdof ; 
displ=[num'  disp] 

numm=l : 1 : nel / 
stresses^ [numm 1  stress] 
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FEMFlk 


nel=7 ; 
nnel=2 ; 
ndof=2 ; 
nnode=5 ; 

sdof =nnode*ndof ; 

gcoord ( 1 , 1 ) =0 ; gcoord ( 1 , 2 ) =0 ; 
gcoord (2 , 1) =0 ; gcoord (2 , 2 ) =12  ; 
gcoord (3, 1) =6 . 3 65; gcoord (3, 2) =  12; 
gcoord (4,1)=6.365; gcoord (4 , 2 ) =0 ; 
gcoord  (5 , 1 ) =3 . 1825 ;  gcoord  (5 , 2 ) =6; 

elprop (1, 1) =10300000; 
elprop ( 1 , 2 ) = . 25 ; 
elprop (2 ,  1)=10300000; 
elprop  (2  ,  2 )  =  . 2  5 ; 
elprop (3,1) =10300000; 
elprop  (3, 2)  =  . 25; 
elprop (4, 1)=10300000; 
elprop  ( 4  ,  2 )  =  . 2  5 ; 
elprop (5, 1) =10300000; 
elprop ( 5 ,  2)  =  .25; 
elprop (6, 1) =10300000; 
elprop ( 6 ,  2)  =  .25; 
elprop (7 , 1)=10300000; 
elprop  ( 7 , 2 )  =  . 2  5 ; 

nodes ( 1 , 1 ) =1 ; nodes (1,2) =2; 
nodes (2,1) =2 ; nodes (2 , 2) =5; 
nodes (3,1) =2; nodes (3, 2) =3; 
nodes (4,1) =3 ; nodes (4,2) =5 ; 
nodes (5,1) =4; nodes (5, 2) =5; 
nodes (6,1) =4 ; nodes (6, 2) =1; 
nodes (7,1) =5; nodes (7, 2) =1; 

bcdof ( 1 ) =5 ; 
bcval ( 1 ) =0 ; 
bcdof (2) =6; 
bcval (2 ) =0 ; 
bcdof (3) =7; 
bcval (3 ) =0 ; 
bcdof (4) =8; 
bcval ( 4 ) =0 ; 

ff=zeros (sdof , 1) ; 
kk=zeros (sdof, sdof)  ; 
index=zeros (nnel*ndof ,  1)  ; 
elf orce=zeros (nnel*ndof ,  1)  ; 
eldisp=zeros (nnel*ndof , nnel*ndof ) 
stress=zeros (nel, 1)  ; 


f  f  (1)  —600; 
ff (2)=600; 
ff (3) =-600; 
ff (4) =600; 

for  iel=l:nel 

nd ( 1 ) =nodes ( iel , 1 ) ; 
nd (2 ) =nodes  (iel, 2 ) ; 

xl=gcoord (nd ( 1 ) , 1 ) ; yl=gcoord (nd  ( 1 ) ,  2 ) 
x2=gcoord (nd ( 2 ) , 1 ) ; y2=gcoord (nd ( 2 ) , 2 ) 

leng=sqrt ( (x2-xl) A2+ (y2-yl) A2) ; 

if  (x2-xl ) ==0 ; 

beta=2*atan (1) ; 

else 

beta=atan ( (y2-yl) / (x2-xl) ) ; 

end 

el=elprop  (iel, 1) ; 
area=elprop  ( iel ,  2 ) ; 

[ index] =FEELDOF (nd, nnel , ndof) ; 

k=FETRUSS2 (el, lengf area, 0 , beta, 1 ) ; 

kk=FEASMBLl (kk, k, index)  ; 


end 

[kk, f f ] =FEAPLYC2 (kk, f f , bcdof , bcval) ; 

disp=kk\f f ; 

for  iel=l:nel 

nd (1) =nodes  (iel, 1) ; 
nd (2 ) =nodes (iel, 2 ) ; 

xl=gcoord (nd ( 1 ) , 1 ) ; yl=gcoord  (nd  ( 1 ) , 2 ) 
x2=gcoord (nd (2)  ,  1 ) ; y2=gcoord (nd ( 2 )  ,  2 ) 

leng=sqrt ( (x2-xl) A2+ (y2-yl) A2) ; 

if  (x2-xl ) ==0 ; 

beta=2^atan (1) ; 

else 

beta=atan ( (y2-yl) / (x2-xl) ) ; 

end 
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el=elprop  (iel, 1) / 
area=elprop  (iel , 2 ) ; 
index=FEELDOF (nd, nnel , ndof ) ; 

k=FETRUSS2 (el , leng, area, 0 , beta, 1 ) ; 

for  i=l :  (nnel*ndof) 

eldisp (i) =disp (index (i) ) ; 

end 

elf orce=k*eldisp; 

stress (iel) =sqrt (el force (1) A2+elforce (2) A2) /area 

if ( (x2-xl) *elforce (3) )<0; 

stress (iel) ^-stress  (iel) ; 

end 


end 

num=l : 1 : sdof ; 
displ=[num'  disp] 

numm=l : 1 : nel / 
stresses^ [numm 1  stress] 
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APPENDIX  F 


stress  and  deflections  from  compound  loading  of  configuration  f 


FI 

Fla 

Fib 

Flc 

Fid 

Fie 

Flf 

Fig 

displ  = 

displ  s 

displ  = 

displ  s 

displ = 

displ  - 

displ  - 

displ  - 

1,0000  -0.0000 

2.0000  0,0000 

1.0000 

2.0000 

0,0000 

0,0000 

1.0000  -0,0000 

2.0000  0.0000 

1,0000 

2,0000 

0.0000 

-0.0000 

1.0000 

2,0000 

0.0000 

,0.0000 

1.0000 

2.0000 

■0.0000 

■0.0000 

1.0000 

2.0000 

0.0003 

■0,0045 

1,0000 

2,0000 

0.0022 

-0.0047 

3-0000  -0.0003 

3.0000 

0.0242 

3.0000  -0.0245 

3,0000 

0.0002 

3.0000 

0.Q24S 

3.0000 

-O.02d2 

3.0000 

-0.0008 

3.CCCC 

0-0006 

4.0000  -0.0017 

4.0000 

O.OO30 

4.0000  -0,0070 

4.0000 

0.0017 

4,0000 

0-0070 

4.0000 

-0-0036 

4.0000 

0,0045 

4.0000 

-00043 

5,0000  0  0002 

5.0000 

0,0245 

5.0000  -0,0242 

5.0000 

-0,0002 

5,0000 

0.0242 

5.0000 

-0,0245 

5-0000 

0.0000 

5,0000 

-0,0000 

6X000  -0,0017 

6,0000 

-0.0070 

6,0000  0,0036 

6,0000 

0,0017 

6,0000 

0,0036 

6.0000 

0,0070 

6-0000 

-0,0000 

6,0000 

0,0000 

7.0000  0 

7,0000 

-0.0000 

7,0000  0.0000 

7,0000 

0 

7,0000  ■ 

■0,0000 

7,0000 

0,0000 

7,0000 

0,0000 

7,0000 

-0,0000 

£.0000  O 

3,0000 

0,0000 

3.0000  0.0000 

6,0000 

O 

6,0000  ■ 

■0,0000 

3,0000 

-0,0000 

8,0000 

0,0000 

8.0000 

-0,0000 

9.OO00  0.0000 

9,0000 

0,0072 

9.0000  -0,0072 

9,0000 

-0.OO0O 

9.0000 

0,0072 

9,0000 

-0.0072 

9.0000 

-0.0000 

9,0000 

0.0005 

10.0000  -o.ooos 

10.0000 

-0,0008 

10,0000  -0,0003 

10.0000 

0,0008 

10.0000 

0.0003 

10.0000 

0.0008 

10.0000 

-0.0020 

10,0000 

0.0020 

stresses = 

stresses  = 

stresses  - 

stresses  = 

stresses  - 

stresses  = 

stresses  = 

stresses  = 

l.Oe+003  m 

0.0010  l-ms 

l.De+003 

0.0010 

3.0752 

l,Qe+003  ■ 

0.0010  5.9743 

1.0S+OQ3 

a.o0io 

1.4495 

1,064003 

0.13010 

5.9743 

l.Oe+003 

0.0010 

3.0752 

l.Oe+003 

0.0010 

0 

1.0e4003 

0.0010 

0.3216 

0.0020  -1.0759 

0,0020 

6,197* 

0.0020  4.0459 

0.OO20 

1.0759 

0.0020  -4.0459 

0.0020 

6,1978 

0.0020 

-2,7167 

0  0020 

-3.0807 

00030  0  5041 

0.0030 

05041 

0.0030  0.5041 

0.0030 

-0  5041 

0.0030 

-0-5041 

0.0030 

-0.5041 

0.0030 

1-2730 

0  0030 

-0.9564 

0.0040  -1,0759 

0.0040 

4,0459 

0.0040  -6,1976 

0,0040 

3.0759 

0,0040 

6.1976 

0,0040  -4,0459 

0.0040 

2,7167 

0,0040 

2,3527 

0.0O50  1.4435 

0.0050 

5,9743 

0,0050  3 ,0752 

0,0050 

1.4495 

0,0050 

3,0752 

0,0050 

5,9743 

0,0050 

-2.7167 

0,0650 

-3,0807 

0.0O6O  -1,0759 

0.0060 

-6,1978 

0.0060  4.0459 

0,0060 

1,0759 

0.OO60 

■4.0459 

0.0060 

6,1978 

0.0060 

-:.275G 

0,0060 

-3,5024 

0.0070  -1,0759 

0,0070 

4,0459 

0.0070  -6.1978 

0,0070 

1.0759 

0.0070 

6. 1978 

0.0070  4.0459 

0.0070 

2,7167 

0,0070 

2.3527 

Fill 

Fli 

Flj 

Flk 

displ  - 

displ  = 

displ  - 

displ  s 

1.0000  -0.0006 

1.0000 

"O.OOOE 

1.0000  0.0006 

1,0000 

-0.0022 

2.0000  0.0043 

2.0000 

0,0045 

2.0000  O.0O43 

2.0000 

0.0047 

3.0000  -0.0032 

4  0000  -0.0047 

3.0000 

4.0000 

O.0O0& 

0X045 

3.0000  0.0052 

4.0000  00047 

3.0000 

4.0000 

-0,0006 

0.0043 

5,0000  0.0000 

5,0000 

-0.0000 

5.0000  -0,0000 

5,0000 

0.0000 

0,0000  -0,0000 

6,0000 

0,0000 

6,0000  0,0000 

6,0000 

-0.0O0O 

7,0000  -0,0000 

7,0000 

-0.0000 

7,0000  0  OOOO 

7,0000 

0,0000 

£.0000  -0,0000 

8,0000 

-0,0000 

3.0000  O.0OQO 

S,OOCO 

0,0000 

9,0000  -0,0005 

9,0000 

0,0000 

9.0000  O.O0O5 

9,0000 

-0.0005 

10.0000  -  0.0020 

10.0000 

0.0020 

10,0000  0.0020 

10.0000 

0.0020 

stresses = 

stresses  = 

stresses  = 

stresses  = 

1 Xe+003  * 

l.Oe+003 

* 

l,Ge+003  * 

l.Oe+003 

* 

0,0010  0  3216 

0.0020 

0 

0.0010  0.3216 

0.0010 

0.3216 

O.0O2O  -2.3527 

0,0020 

2.7167 

0.0020  2,3527 

0.0020 

3.0S07 

O-003O  3-5024 

0.0030 

-1.2730 

0.0030  -3.5024 

0.OO30 

0.9564 

0,0040  3.0807 

0.0040 

-2.7167 

0.0040  -3,0607 

0,0040 

-2,3527 

0.0050  -2,3527 

0.0050 

2,7167 

0.0050  2.3527 

0,0050 

3,0807 

O.006O  0.9564 

0,0060 

:.2730 

O.OO60  -0,9564 

□,0060 

3.5024 

O.0O70  3.0S07 

0,0070 

-2.7167 

0.0070  -3.0807 

0,0070 

-2.3527 
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APPENDIX  G 


Absolute 

Max 

deflection 

Absolute 

Max 

Stress 

a(psi) 

A 

(in  A2) 

F=ct*A 

(lbs) 

k=F/5 

(Ibs/in) 

%  of 
full  k 

Per  (lbs) 

Status 

FI 

0.0017 

1.4495 

1449.5 

0.25 

362.375 

213,162 

15.60 

919.2069 

Safe 

Flal 

0.0245 

6.1978 

6197.8 

0.25 

1549.45 

63,243 

4.63 

919.2069 

Flbl 

0.0245 

6.1978 

6197.8 

0.25 

1549.45 

63,243 

4.63 

919.2069 

Flcl 

0.0017 

1.4495 

1449.5 

0.25 

362.375 

213,162 

15.60 

919.2069 

Safe  | 

Fldl 

0.0245 

6.1978 

6197.8 

0.25 

1549.45 

63,243 

4.63 

919.2069 

Flel 

0.0245 

6.1978 

6197.8 

0.25 

1549.45 

63,243 

4.63 

919.2069 

Flfl 

0.0045 

2.7167 

2716.7 

0.25 

679.175 

150,928 

11.05 

919.2069 

Safe 

Flgl 

0.0047 

3.5024 

3502.4 

0.25 

875.6 

186,298 

13.64 

919.2069 

Safe 

Flhl 

0.0047 

3.5024 

3502.4 

0.25 

875.6 

186,298 

13.64 

919.2069 

Safe 

Flil 

0.0045 

2.7167 

2716.7 

0.25 

679.175 

150,928 

11.05 

919.2069 

Safe 

Fiji 

0.0047 

3.5024 

3502.4 

0.25 

875.6 

186,298 

13.64 

919.2069 

Safe 

Flkl 

0.0047 

3.5024 

3502.4 

0.25 

875.6 

186,298 

13.64 

919.2069 

Safe 
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